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[N],b[M],op,x,y,z,s[N],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;
}