正题

题目链接:https://www.luogu.com.cn/problem/P3700

---

题目大意

一个$n\ast n$个数的数字表格，开始位置$(a,b)$上的是$a\ast b$。数字表格需满足以下条件

1. 对于任意$(a,b)$有$f(a,b)=f(b,a)$
2. 对于任意$(a,b)$有$b\ast f(a,a+b)=(a+b)\ast f(a,b)$

$m$次修改数字表格中的一个数，你需要调整其他数字使其满足条件，然后求前$k$行前$k$列的数字和。

---

解题思路

一道神仙题，首先我们拿出这个不伦不类的式子$b\ast f(a,a+b)=(a+b)\ast f(a,b)$化成
$$\frac{f(a,a+b)}{a(a+b)}=\frac{f(a,b)}{ab}$$
发现这个式子和更相减损法很像，反正最后递归下去可以化成$$\frac{f(a,b)}{ab}=\frac{f(gcd(a,b),gcd(a,b))}{gcd(a,b{)}^2}$$
定义$d=gcd(a,b)$，也就是$f(a,b)=\frac{f(d,d)ab}{gcd(a,b{)}^2}$
那么$k\ast k$的数字和就是$$\sum _{i=1}^k{f}_{i,i}\sum _{x=1}^k\sum _{y=1}^k\frac{xy}{i^2}[gcd(x,y)==i]$$
然后化成$$\sum _{i=1}^k{f}_{i,i}\sum _{x=1}^{\lfloor \frac{k}{i}\rfloor }\sum _{y=1}^{\lfloor \frac{k}{i}\rfloor }xy[gcd(x,y)==1]$$
之后我们有结论是后面那个东西${\sum }_{x=1}^n{\sum }_{y=1}^{n}xy[gcd(x,y)==1]={\sum }_{i=1}^n\phi (i)i^2$
具体证明的话就是

令$H(x)={\sum }_{x=1}^n{\sum }_{y=1}^{n}xy[gcd(x,y)==1]$
那么$H(x)$比$H(x-1)$多了一个$2x{\sum }_{y=1}^n[gcd(x,y)==1]y$
也就是多了一个$x\ast 2\ast \frac{\phi (x)x}{2}$

所以可以线性预处理这个东西

我们就可以对后面那个东西做整除分块，我们就需要处理$f_{i,i}$的前缀和

我们一次操作需要进行$\sqrt{n}$次询问前缀和，但是修改只需要一次，所以我们可以平衡一下时间。也就是降低询问复杂度提升修改复杂度，用分块存储块中每个位置的块内前缀和和块的前缀和即可。

时间复杂度$O(m\sqrt{n})$

---

$code$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
const ll N=4e6+10,P=1e9+7;
ll n,m,T,w[N],v[N],sum[N],L[N],R[N];
ll cnt,pri[N],phi[N],g[N],pos[N];
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
void Prime(){phi[1]=1;for(ll i=2;i<=n;i++){if(!v[i])pri[++cnt]=i,phi[i]=i-1;for(ll j=1;j<=cnt&&i*pri[j]<=n;j++){v[i*pri[j]]=1;if(i%pri[j]==0){phi[i*pri[j]]=phi[i]*pri[j];break;}phi[i*pri[j]]=phi[i]*phi[pri[j]];}}for(ll i=1;i<=n;i++)g[i]=(g[i-1]+phi[i]*i%P*i%P)%P;return;
}
ll GetS(ll x){if(!x)return 0ll;return (sum[pos[x]-1]+w[x])%P;
}
int main()
{scanf("%lld%lld",&m,&n);Prime();T=sqrt(n);for(ll i=1;i<=n;i++)v[i]=i*i%P;for(ll i=1;i<=T;i++)L[i]=R[i-1]+1,R[i]=i*T;if(T*T!=n)T++,L[T]=R[T-1]+1,R[T]=n;for(ll i=1;i<=T;i++){for(ll j=L[i];j<=R[i];j++)pos[j]=i,w[j]=(v[j]+w[j-1]*(j>L[i]))%P;sum[i]=(sum[i-1]+w[R[i]])%P;}while(m--){ll a,b,x,k;scanf("%lld%lld%lld%lld",&a,&b,&x,&k);x%=P;ll d=__gcd(a,b),y=pos[d];v[d]=x*d*d%P*power(a*b%P,P-2)%P;for(ll i=d;i<=R[y];i++)w[i]=((i>L[y])*w[i-1]+v[i])%P;for(ll i=y;i<=T;i++)sum[i]=(sum[i-1]+w[R[i]])%P;ll ans=0;for(ll l=1,r;l<=k;l=r+1){r=k/(k/l);ans=(ans+g[k/l]*(GetS(r)-GetS(l-1)+P)%P)%P;}printf("%lld\n",ans);}return 0;
}