HDU 3030

dp+树状数组

递推公式很简单 

dp[i]表示以第 i 个数为结尾的上升子序列的个数。最后的答案即是dp[1]+...+dp[n]

对速度进行离散化处理后，通过树状数组可以方便的求出

for i = 0 to n-1
print A[i mod m]
A[i mod m] = (X * A[i mod m] + Y * (i + 1)) mod Z

这部分描述就是个烟雾弹。。。

#include<stdio.h>
#include<iostream>
#include<memory.h>
#include<algorithm>
using namespace std;

#define Maxn 500010
#define MOD 1000000007
int n,m,MaxSp;
long long X,Y,Z,a[Maxn],sl[Maxn],tree[Maxn];
struct node
{
    long long s;
    int num;
}tmpsl[Maxn];

bool cmp(node a,node b)
{
    return a.s<b.s;
}

//求2^k
int lowbit(int t)
{
    return t & ( t ^ ( t - 1 ) );
}
//求前n项和
long long Sum(long long T[],int end)
{
   long long sum = 0;
   while(end > 0)
  {
     sum = (sum+T[end])%MOD;
     end -= lowbit(end);
  }
  return sum;
}
//增加某个元素的大小
void Add(long long T[],int pos,long long num)
{
   while(pos <= MaxSp)
  {
     T[pos] = (T[pos]+num)%MOD;
     pos += lowbit(pos);
  }
}

int main()
{
    int T;
    scanf("%d",&T);
    for(int Case=1;Case<=T;Case++)
    {
        memset(tree,0,sizeof(tree));
        scanf("%d %d %I64d %I64d %I64d",&n,&m,&X,&Y,&Z);
        for(int i=0;i<m;i++)
            scanf("%I64d",&a[i]);
        for(int i=0;i<n;i++)
        {
            sl[i]=a[i%m];
            tmpsl[i+1].s=sl[i],tmpsl[i+1].num=i;
            a[i%m]=(X*a[i%m]+Y*(i+1))%Z;
        }
        sort(tmpsl+1,tmpsl+1+n,cmp);
        sl[tmpsl[1].num]=MaxSp=1;
        for(int i=2;i<=n;i++)
        {
            if(tmpsl[i].s==tmpsl[i-1].s)
                sl[tmpsl[i].num]=MaxSp;
            else
                sl[tmpsl[i].num]=++MaxSp;
        }

        long long ans=0,add;
        for(int i=0;i<n;i++)
        {
            add=Sum(tree,sl[i]-1)+1;
            ans=(ans+add)%MOD;
            Add(tree,sl[i],add);
        }
        printf("Case #%d: %I64d\n",Case,ans);
    }
    return 0;
}