10.29 T1 max 單調棧+貢獻法+二階差分
阿新 • • 發佈:2018-11-02
solution
- 對於每個i用單調棧求出小於i的第一個比他大的數l[i],以及大於i的第一個比他大的數r[i]。
- i一定是l[i]-r[i]中最大的數,1-a+1分別加(1-a+1)*a[i],a+2 - b+1加 a+1 a[i] a+2-a+b+1 遞減加
- 用二階差分o(1)完成修改
- 用兩遍字首和求出答案
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; typedef long long ll; const int N=1e7+100; const int inf=0x3f3f3f3f; int a[N]; int n,p; void gen() { int x,y,z; scanf("%d%d%d%d%d",&n,&p,&x,&y,&z); int tmp=1; for(int i=1;i<=n;++i) { tmp=((long long)tmp*x+y)%z; a[i]=tmp; } } int s[N],top; int c[N]; int l[N],r[N]; int main() { freopen("max.in","r",stdin); freopen("max.out","w",stdout); gen(); int i,j; a[n+1]=inf; a[0]=inf; top=1; for(i=1;i<=n+1;i++) { while(top&&a[i]>=a[s[top]]) { r[s[top--]]=i; } s[++top]=i; } top=1; for(i=n;i>=0;i--) { while(top&&a[i]>a[s[top]]) { l[s[top--]]=i; } s[++top]=i; } int x,y; for(i=1;i<=n;i++) { x=i-l[i]; y=r[i]-i; x--;y--; if(x>y) swap(x,y); c[1]=(c[1]+a[i])%p; c[x+2]=(c[x+2]-a[i]+p)%p; c[y+2]=(c[y+2]-a[i]+p)%p; c[x+y+3]=(c[x+y+3]+a[i])%p; } for(int i=1;i<=n;++i)c[i]=(c[i]+p)%p; for(i=1;i<=n;i++) c[i]=(c[i-1]+c[i])%p; for(i=1;i<=n;i++) c[i]=(c[i-1]+c[i])%p; int ans=0; for(i=1;i<=n;i++) ans=ans^c[i]; cout<<ans<<endl; return 0; }