【51nod1376】最长递增子序列的数量

【51nod1376】最长递增子序列的数量
好题。
首先nlgn的方法求以每个位置的数为结尾的最长递增子序列，记为dp[i]
显然dp值为x的数的答案数是由dp值为x-1的数贡献的。
例如现在要求dp值为x的数i的答案数，从dp值为x-1的数中找一个数j，当满足j的坐标小于i，且j的值小于i的值时j的答案数就可以加到i的答案里，对i产生贡献。由此转化为偏序问题，用cdq分治求解（树状数组也行）。

#include<cstdio>
#include<vector>
#include<algorithm>
std::vector<int>s[50005];
struct node{int x,y,cnt,id;}p[50005],q[50005];
bool comp(node a,node b){return a.x<b.x;}
int n,f[50005],mx,ans[50005],a[50005];const int mo=1000000007;
int max(int a,int b){return a>b?a:b;}
void cdq(int l,int r)
{
    if(l==r)return;
    int mid=l+r>>1,i,j,cnt=l-1,sum=0;
    cdq(l,mid);cdq(mid+1,r);
    for(i=l,j=mid+1;i<=mid&&j<=r;)
     if(q[i].y<q[j].y)
     {
        if(!q[i].id)(sum+=q[i].cnt)%=mo;
        p[++cnt]=q[i++];
     }else
     {
        if(q[j].id)(ans[q[j].id]+=sum)%=mo;
        p[++cnt]=q[j++];
     }
    while(j<=r){if(q[j].id)(ans[q[j].id]+=sum)%=mo;p[++cnt]=q[j++];}
    while(i<=mid)p[++cnt]=q[i++];
    for(int k=l;k<=r;k++)q[k]=p[k];
}
int find(int l,int r,int x)
{
    int ans=0;
    while(l<=r)
    {
        int mid=l+r>>1;
        if(f[mid]<x)ans=mid,l=mid+1;else r=mid-1;
    }
    return ans; 
}
int main()
{
    scanf("%d%d",&n,&f[1]);a[1]=f[1];s[1].push_back(1);int t=1;
    for(int i=2,x;i<=n;i++)
    {
        scanf("%d",&x);int j=find(1,t,x);//找最大的小于x的数
        if(f[j+1]>x)f[j+1]=x;if(j==t)f[++t]=x;
        s[j+1].push_back(i);a[i]=x;
    }
    for(int i=0;i<s[1].size();i++)ans[s[1][i]]=1;
    for(int i=2;i<=t;i++)
    {
        int cnt=0;
        for(int j=0;j<s[i-1].size();j++)
         q[++cnt]=(node){s[i-1][j],a[s[i-1][j]],ans[s[i-1][j]],0};//修改 
        for(int j=0;j<s[i].size();j++)
         q[++cnt]=(node){s[i][j],a[s[i][j]],0,s[i][j]};//查询 
        std::sort(q+1,q+1+cnt,comp);//按坐标排 
        cdq(1,cnt); 
    }
    int sum=0;
    for(int i=0;i<s[t].size();i++)(sum+=ans[s[t][i]])%=mo;
    printf("%d\n",sum); 
}