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
这部分描述就是个烟雾弹。。。
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;
}