TOJ 4146 seq

题目链接：http://acm.tju.edu.cn/toj/showp4146.html

一道怎么看都很简单的一道题，比赛的时候就是想不到怎么去优化TT。

题意是给我们一串数，然后给我们两种操作，一种是把所有i中对b[v]取模得到的结果为x的a[i]加上一个add，一种是所有i中对b[v]整除得到x的a[i]加上一个add。

小编直接说优化方法了，对于前一个取模，可以先把结果存在一个c[v][x]的二维数组中，表示的是对b[v]取模得到x的数加的值，最后再统一加，对于第二个整除的，显然结果是一段数值，我们可以用前缀和的方法来解决，用s[i]存储下要加的数，每一次s[L] = add,s[R] = -add，然后s[i] += s[i-1]这样就能在O(n)的复杂度下完成这个加法。

#include<cstdio>
#include<algorithm>
using namespace std;
const int N=100010,M=30;
int T,n,m,q,i,j,k,l,r,a
,b[M],op,x,y,z,s
,c[M][1010];
int main()
{
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
for(i=0;i<n;i++) scanf("%d",&a[i]);
for(i=0;i<m;i++) scanf("%d",&b[i]);
for(i=0;i<=n;i++) s[i]=0;
for(i=0;i<m;i++)
for(j=0;j<b[i];j++)
c[i][j]=0;
scanf("%d",&q);
while(q--)
{
scanf("%d%d%d%d",&op,&x,&y,&z);
if(op==0)
c[x][y]+=z;
else
{
x=b[x];
l=max(x*y,0);
r=min(x*(y+1)-1,n-1);
if(l<=r)s[l]+=z,s[r+1]-=z;
}
}
for(i=0;i<m;i++)
for(j=0;j<b[i];j++)
for(k=j;k<n;k+=b[i])
a[k]+=c[i][j];
for(i=0;i<n;i++)
{
if(i)s[i]+=s[i-1];
printf("%d%c",a[i]+s[i],i==n-1?'\n':' ');
}
}
return 0;
}