3688: 折线统计
Description
二维平面上有n个点(xi, yi),现在这些点中取若干点构成一个集合S,对它们按照x坐标排序,顺次连接,将会构成一些连续上升、下降的折线,设其数量为f(S)。如下图中,1->2,2->3,3->5,5->6(数字为下图中从左到右的点编号),将折线分为了4部分,每部分连续上升、下降。
现给定k,求满足f(S) = k的S集合个数。
Input
第一行两个整数n和k,以下n行每行两个数(xi, yi)表示第i个点的坐标。所有点的坐标值都在[1, 100000]内,且不存在两个点,x坐标值相等或y坐标值相等
Output
输出满足要求的方案总数 mod 100007的结果
Sample Input
5 1
5 5
3 2
4 4
2 3
1 1
Sample Output
19
【解题报告】
首先我们将点按照x轴排序。
令dp[i][j][0]和dp[i][j][1]分别表示前i个点,选择j段,最后一段是下降或者上升的方案数。
状态转移方程如下:
dp[i][j][0]=∑(dp[k][j][0]+dp[k][j-1][1]) (k
/**************************************************************
Problem: 3688
User: onepointo
Language: C++
Result: Accepted
Time:1240 ms
Memory:25044 kb
****************************************************************/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 100010
#define mod 100007
#define lowbit(x) ((x)&(-x))
struct Point
{
int x,y;
bool friend operator < (Point a,Point b)
{return a.x<b.x;}
}p[N];
int n,m,k;
int dp[N][15][2],bit[N][15][2];
void add(int p,int a,int b,int v)
{
for(int i=p;i<N;i+=lowbit(i))
bit[i][a][b]=(bit[i][a][b]+v)%mod;
}
int query(int p,int a,int b)
{
int sum=0;
for(int i=p;i;i-=lowbit(i))
sum=(sum+bit[i][a][b])%mod;
return sum;
}
int main()
{
scanf("%d%d",&n,&k);
for(int i=1;i<=n;++i)
scanf("%d%d",&p[i].x,&p[i].y);
sort(p+1,p+n+1);
for(int i=1;i<=n;++i)
{
dp[i][0][0]=dp[i][0][1]=1;
add(p[i].y,0,0,1);add(p[i].y,0,1,1);
for(int j=1;j<=k;++j)
{
dp[i][j][0]=(dp[i][j][0]+query(p[i].y-1,j,0)+query(p[i].y-1,j-1,1))%mod;
dp[i][j][1]=(dp[i][j][1]+query(N-1,j,1)-query(p[i].y,j,1)+query(N-1,j-1,0)-query(p[i].y,j-1,0))%mod;
if(dp[i][j][1]<0) dp[i][j][1]+=mod;
add(p[i].y,j,0,dp[i][j][0]);
add(p[i].y,j,1,dp[i][j][1]);
}
}
int ans=0;
for(int i=1;i<=n;i++)
ans=(ans+dp[i][k][0]+dp[i][k][1])%mod;
printf("%d\n",ans);
return 0;
}