4815: [Cqoi2017]小Q的表格
Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 832 Solved: 342
[Submit][Status][Discuss]Description
小Q是个程序员。作为一个年轻的程序员，小Q总是被老C欺负，老C经常把一些麻烦的任务交给小Q来处理。每当小Q不知道如何解决时，就只好向你求助。为了完成任务，小Q需要列一个表格，表格有无穷多行，无穷多列，行和列都从1开始标号。为了完成任务，表格里面每个格子都填了一个整数，为了方便描述，小Q把第a行第b列的整数记为f(a,b),为了完成任务，这个表格要满足一些条件：(1)对任意的正整数a,b,都要满足f(a,b)=f(b,a)；(2)对任意的正整数a,b,都要满足b×f(a,a+b)=(a+b)*f(a,b)。为了完成任务，一开始表格里面的数很有规律，第a行第b列的数恰好等于a*b，显然一开始是满足上述两个条件的。为了完成任务，小Q需要不断的修改表格里面的数，每当修改了一个格子的数之后，为了让表格继续满足上述两个条件，小Q还需要把这次修改能够波及到的全部格子里都改为恰当的数。由于某种神奇的力量驱使，已经确保了每一轮修改之后所有格子里的数仍然都是整数。为了完成任务，小Q还需要随时获取前k行前k列这个有限区域内所有数的和是多少，答案可能比较大，只需要算出答案mod1,000,000,007之后的结果。Input
第1行包含两个整数m,n,表示共有m次操作，所有操作和查询涉及到的行编号和列编号都不超过n。接下来m行，每行4个整数a,b,x,k表示把第a行b列的数改成x，然后把它能够波及到的所有格子全部修改，保证修改之后所有格子的数仍然都是整数，修改完成后计算前k行前k列里所有数的和。1<=m<=10000,1<=a,b,k<=N<=4*10^6,0<=x<=10^18Output
输出共m行，每次操作之后输出1行，表示答案mod 1,000,000,007之后的结果。
Sample Input
3 3
1 1 1 2
2 2 4 3
1 2 4 2Sample Output
9
36
14
一开始,表格的前3行前3列如图中左边所示,前2次操
作后表格没有变化,第3次操作之后的表格如右边所示。
1 2 3 2 4 6
2 4 6 4 4 12
3 6 9 6 12 9
HINT
Source
[Submit][Status][Discuss]
https://blog.sengxian.com/solutions/bzoj-4815
题目里的式子有点像辗转相除，所以可以猜测与gcd有关。
最后的式子利用的是i和j的对称性。
最后要求的是$sumlimits_{d=1}^{k}f(d)g(frac{k}{d})$，其中$g(n)=sumlimits_{i=1}^ni^2phi(i)$。
因为f(d)是在线修改的，所以需要维护前缀和，容易想到树状数组，但是会T。
注意到对于单次操作，修改个数是$O(1)$的，询问个数是$O(sqrt n)$的，所以我们需要的是修改$O(sqrt n)$，询问$O(1)$的算法。
分块可以满足这个需求，注意需要前缀和与整块的差分。
1 #include<cmath> 2 #include<cstdio> 3 #include<algorithm> 4 #define rep(i,l,r) for (int i=l; i<=r; i++) 5 typedef long long ll; 6 using namespace std; 7 8 const int N=4000100,mod=1000000007; 9 bool pr[N]; 10 ll x; 11 int n,m,Q,bl,a,b,k,tot,c[N],pos[N],phi[N],p[N],v[N],f[N],g[N]; 12 struct B{ int l,r; }B[2010]; 13 int gcd(int a,int b){ return (!b) ? a : gcd(b,a%b); } 14 15 int ksm(int a,int b){ 16 int res; 17 for (res=1; b; a=1ll*a*a%mod,b>>=1) 18 if (b & 1) res=1ll*a*res%mod; 19 return res; 20 } 21 22 void pre(){ 23 phi[1]=1; 24 rep(i,2,n){ 25 if (!pr[i]) p[++tot]=i,phi[i]=i-1; 26 for (int j=1; j<=tot && i*p[j]<=n; j++){ 27 int t=i*p[j]; pr[t]=1; 28 if (i%p[j]) phi[t]=phi[i]*(p[j]-1); 29 else { phi[t]=phi[i]*p[j]; break; } 30 } 31 } 32 rep(i,1,n) g[i]=(g[i-1]+1ll*i*i%mod*phi[i])%mod; 33 } 34 35 void init(){ 36 bl=sqrt(n); m=(n-1)/bl+1; 37 rep(i,1,n) pos[i]=(i-1)/bl+1,v[i]=1ll*i*i%mod,f[i]=(f[i-1]+v[i])%mod; 38 rep(i,1,m) B[i].l=(i-1)*bl+1,B[i].r=i*bl; 39 } 40 41 void mdf(int x,int k){ 42 int r=B[pos[x]].r,d=(k-v[x]+mod)%mod; v[x]=k; 43 if (!k) return; 44 rep(i,x,r) f[i]=(f[i]+d)%mod; 45 rep(i,pos[x]+1,m) c[i]=(c[i]+d)%mod; 46 } 47 48 int que(int x){ return (c[pos[x]]+f[x])%mod; } 49 50 int solve(int n){ 51 int ans=0,lst=0,r=0,now=0; 52 for (int i=1; i<=n; i=r+1,lst=now){ 53 r=n/(n/i); now=que(r); 54 ans=(ans+1ll*(now-lst+mod)*g[n/i]%mod)%mod; 55 } 56 return ans; 57 } 58 59 int main(){ 60 freopen("bzoj4815.in","r",stdin); 61 freopen("bzoj4815.out","w",stdout); 62 scanf("%d%d",&Q,&n); pre(); init(); 63 while (Q--){ 64 scanf("%d%d%lld%d",&a,&b,&x,&k); x%=mod; 65 int d=gcd(a,b); mdf(d,1ll*x*d%mod*d%mod*ksm((1ll*a*b)%mod,mod-2)%mod); 66 printf("%d ",solve(k)); 67 } 68 return 0; 69 }