[容斥 反演] 51nod1678 lyk与gcd

只会无脑反演…

$$\sum_{j}[gcd(i,j)==1]*a[j]=\sum_{j} \sum_{k|i,k|j} \mu(k) *a[j]=\sum_{k|i} \mu(k) \sum_{j} a[k*j]$$
设 $g[k]=\sum_{k|j} a[j]$
一个 $log$ 预处理约数。修改时维护 $g(x)$ ，询问就暴力好了。
$O(n \sqrt n)$

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=100005,maxv=100005,N=100000;
typedef long long LL;
int p[maxv],mu[maxv];
vector<int> lst[maxv];
bool vis[maxv];
void Pre(){
    mu[1]=1;
    for(int i=2;i<=N;i++){
        if(!vis[i]) p[++p[0]]=i, mu[i]=-1;
        for(int j=1;j<=p[0]&&(LL)i*p[j]<=N;j++){
            vis[i*p[j]]=true;
            if(i%p[j]==0){ mu[i*p[j]]=0; break; }
            mu[i*p[j]]=-mu[i];
        }
    }
    for(int i=1;i<=N;i++)
     for(int j=1;(LL)i*j<=N;j++) lst[i*j].push_back(i);
}
int n,Q,a[maxn],g[maxv];
void Updata(int x,int val){
    for(int i=0;i<lst[x].size();i++) g[lst[x][i]]+=val;
}
int main(){
    freopen("51nod1678.in","r",stdin);
    freopen("51nod1678.out","w",stdout);
    Pre();
    scanf("%d%d",&n,&Q);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]), Updata(i,a[i]);
    while(Q--){
        int _type,x,y; scanf("%d%d",&_type,&x);
        if(_type==1) scanf("%d",&y), Updata(x,y-a[x]), a[x]=y; else{
            LL res=0;
            for(int i=0;i<lst[x].size();i++) res+=g[lst[x][i]]*mu[lst[x][i]];
            printf("%lld\n",res);
        }
    }
    return 0;
}