元月初七hu测 T1.折线统计（dp+树状数组）

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

Hint

对于 10%的数据，n <= 10, k = 1
对于 30%的数据，n <= 1000, k = 1
对于 50%的数据，n <= 1000, k <= 10
对于 100%的数据，n <= 50000，0 < k <= 10

分析：
原题：bzoj3688

简单的，对于30%的数据
$f[i]$表示$i$作为子集中的最后一个点，且序列上升时的方案数
$g[i]$表示$i$作为子集中的最后一个点，且序列下降时的方案数
$f[i]=Σf[j]，j＜i，a[j]＜a[i]$
$g[i]=Σg[j]，j＜i，a[j]＞a[i]$

奇怪的是，我写的这一部分的代码，竟然只有10分？？？
我发现如果直接用暴力枚举30%很稳
但是如果用上述的dp，初始化是什么呢？问题就出在初始化上，没法破。。。所以最终得分只有10

算了，我们并不关心部分分，重点是正解：
$f[i][s]$表示已经有$s$段，最后一个点是$i$，最后一段是上升的方案数
$g[i][s]$表示已经有$s$段，最后一个点是$i$，最后一段是下降的方案数
$f[i][s]= Σf[j][s](j<i，a[j]>a[i])+Σg[j][s-1](j<i，a[j]<a[i])$
$g[i][s]=Σg[j][s](j<i，a[j]>a[i])+Σf[j][s-1](j<i，a[j]>a[i])$

（此时我奇迹般的发现，我其实想出了正解，只是没有完善出来。。。）

在转移的时候，因为具有前缀和的形式，可以用树状数组维护一下
（顺序转移，因此在插入树状数组的时候，就已经满足$j<i$这个条件了）

tip

注意树状数组的大小是由坐标范围确定的
我惊奇的发现，自己已经想到的了dp的基本思路以及用树状数组优化
看来还是要充分的相信自己，fight！

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

const int p=100007;
const int N=100003;
int n,m,c[N][12][2],mxy=0,f[N][12][2];
struct node{
    int x,y;
    bool operator <(const node &a) const
    {
        return x<a.x;
    }
};
node a[N];

void add(int x,int y,int z,int v)
{
    for (int i=x;i<=mxy;i+=(i&(-i)))
        c[i][y][z]=(c[i][y][z]+v)%p;
}

int ask(int x,int y,int z)
{
    int ans=0;
    for (int i=x;i>0;i-=(i&(-i)))
        ans=(ans+c[i][y][z])%p;
    return ans;
}

void solve()
{
    for (int i=1;i<=n;i++)
    {
        f[i][0][0]=f[i][0][1]=1;
        add(a[i].y,0,0,1); add(a[i].y,0,1,1);
        for (int j=1;j<=m;j++)
        {
            f[i][j][0]+=ask(a[i].y-1,j,0)+ask(a[i].y-1,j-1,1);
            f[i][j][1]+=ask(mxy,j,1)-ask(a[i].y,j,1)+ask(mxy,j-1,0)-ask(a[i].y,j-1,0);
            ((f[i][j][0]%=p)+=p)%=p;
            ((f[i][j][1]%=p)+=p)%=p;
            add(a[i].y,j,0,f[i][j][0]); add(a[i].y,j,1,f[i][j][1]);
        }
    }
    int ans=0;
    for (int i=1;i<=n;i++)
    {
        ans=(ans+f[i][m][0])%p;
        ans=(ans+f[i][m][1])%p;
    }
    printf("%d\n",ans);
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y),mxy=max(mxy,a[i].y);
    sort(a+1,a+1+n);
    solve();
    return 0;
}