题链:
http://codeforces.com/problemset/problem/549/F
题解:
分治,链表。
考虑对于一个区间[L,R],其最大值在p位置,
那么答案的贡献就可以分为3部分:
1.[L,p-1]中合法的区间
2.[p+1,R]中合法的区间
3.[L,R]中经过p的合法区间。
前两个因为是相同子问题,可以递归解决。
考虑如何求出3.的答案。
令S[i]表示$\sum_{k=1}^{i}A[i]$(即前缀和)
显然对于一组(l,r)表示的区间[l,r],
当且仅当满足L<=l<=p且p<r<=R(或L<=l<p且p<=r<=R),
并且 S[r]-S[l-1]-A[p]≡0 (mod K)时,才是合法区间。
但是如果直接递归实现,复杂度高达O(N^2)
所以预处理出每个元素在哪个包含它的最大区间里是最大值。
可以用单调栈O(N)维护出。
然后对于每个元素以及刚刚得出来的区间去计算贡献。
同时在计算贡献时,采用的方法是:
枚举一边的元素,然后用形如Query(ql,qr,qv)的方式去询问另一边S[*]==qv的个数
至于枚举哪一边,就是哪边短就枚举哪一边,这样可以保证询问次数为nlogn次的。
然后考虑如何询问。
如果采用主席树,则需要把复杂度再乘上一个log的询问代价,总的复杂度为O(N logN logN)
或者把询问离线成差分形式,插入链表,最后扫描链表O(N logN)得出答案(链表元素有N logN个),总的复杂度为O(N logN)
比如说对于询问Query(l,r,v)拆成两个(v,sign)的形式,我们在链表l-1中加入(v,-1),在链表r中加入(v,1)
最后从i=0到N扫描一遍,首先cnt[S[i]]++,(cnt[x]表示前缀i中x这个值出现了多少次)
然后再扫描i的询问链表,把ans+=cnt[v]*sign即可。
最后的ans既是答案。
代码:
#include<bits/stdc++.h> #define MAXN 300050 #define INF 0x3f3f3f3f using namespace std; int N,K,more; int A[MAXN],S[MAXN],L[MAXN],R[MAXN]; struct LINK{ int lnt; int nxt[MAXN*20],val[MAXN*20],sign[MAXN*20],head[MAXN]; LINK(){lnt=2;} void Add(int u,int v,int s){ if(u<0) return; val[lnt]=v; sign[lnt]=s; nxt[lnt]=head[u]; head[u]=lnt++; } void Query(int l,int r,int v){ Add(l-1,v,-1); Add(r,v,1); } long long Getans(){ long long ret=0; static int cnt[1000050]; for(int i=0;i<=N;i++){ cnt[S[i]]++; for(int j=head[i];j;j=nxt[j]) ret+=sign[j]*cnt[val[j]]; } return ret; } }Q; void prework(){ static int stk[MAXN],top; A[0]=A[N+1]=INF; stk[top=1]=0; for(int i=1;i<=N;i++){ while(A[stk[top]]<=A[i]) top--; L[i]=stk[top]+1; stk[++top]=i; } stk[top=1]=N+1; for(int i=N;i>=1;i--){ while(A[stk[top]]<A[i]) top--; R[i]=stk[top]-1; stk[++top]=i; } for(int i=1;i<=N;i++) S[i]=(S[i-1]+A[i])%K; } void solve(){ int al,ar,qv,ql,qr,sign; for(int i=1;i<=N;i++){ if(L[i]==R[i]) continue; more++; if(i-L[i]+1<=R[i]-i+1) al=L[i]-1,ar=i-1,ql=i,qr=R[i],sign=1; else al=i,ar=R[i],ql=L[i]-1,qr=i-1,sign=-1; for(int j=al;j<=ar;j++){ qv=(1ll*S[j]+sign*A[i]%K+K)%K; Q.Query(ql,qr,qv); } } } int main(){ scanf("%d%d",&N,&K); for(int i=1;i<=N;i++) scanf("%d",&A[i]); prework(); solve(); long long ans=Q.Getans(); printf("%lld\n",ans-more); return 0; }