NOIP模拟 gcd 数学

题面不给……题意：动态增减集合中元素个数，动态求$gcd(i,j)==1$的数的个数。

坦白地讲题解并没完全看懂……于是听$lc$开了一会车（快开完时$xyz$强行砸场导致烂尾），结合着自己理解又推一遍想明白了……

在本题之中，我们设$f(x)$为$gcd(i,j)==x$的数的个数，$g(x)$为$gcd(i,j)==kx$的数的个数。很显然$f(x)=g(x)-f(2x)-f(3x)...$，这样推下去，我们就会发现这个系数其实是符合莫比乌斯函数的，最后大力化简就是：

[ sum_{d}^{} {μ(d)*g(d)} =f(d)]。

现在定义$s(i)$为当前选中的集合中为$i$倍数数的个数，可以看出，$g(i)=s(i)*(s(i)-1)/2$。

然后我们就可以预处理莫比乌斯函数，在每一次修改元素个数时候动态维护$s(i)$。

但是我们暴力求$ans$显然会$T$……于是考虑每一次加入或删除元素对于最后结果贡献。

当我们插入一个元素，对于每一个可以整除它的$i$，$Δans=μ(i)*((s(i)+1)*s(i)-s(i)*(s(i)-1))/2=μ(i)s(i)$。同理，删除时，对于每一个这样的$i$，$Δans=μ(i)*((s(i)-2)*(s(i)-1)-s(i)*(s(i)-1))/2=-μ(i)(s(i)-1)$。然后我们根据这个规律转移即可。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=200005,maxx=500005;
bool notprime[maxx];int prime[maxx],cnt,mu[maxx];
void init()
{
    mu[1]=1;
    for(int i=2;i<=500000;i++)
    {
        if(!notprime[i])prime[++cnt]=i,mu[i]=-1;
        for(int j=1;j<=cnt&&i*prime[j]<=500000;j++)
        {
            notprime[i*prime[j]]=1;
            if(i%prime[j])mu[i*prime[j]]=-mu[i];
            else {mu[i*prime[j]]=0;break;}
        }
    }
}
int n,m,a[maxn],s[maxx],maxv,t[maxn];
int haha()
{
    scanf("%d%d",&n,&m);init();
    for(int i=1;i<=n;i++)scanf("%d",&a[i]),maxv=max(maxv,a[i]);
    long long ans=0;
    for(int i=1;i<=m;i++)
    {
        int x;scanf("%d",&x);int sgn=t[x]?-1:1;
        for(int j=1;j*j<=a[x];j++)
            if(!(a[x]%j))
            {
                if(sgn==1)ans+=1ll*mu[j]*s[j];else ans-=1ll*(s[j]-1)*mu[j];s[j]+=sgn*1;
                if(j*j<a[x])
                {
                    if(sgn==1)ans+=1ll*mu[a[x]/j]*s[a[x]/j];else ans-=1ll*(s[a[x]/j]-1)*mu[a[x]/j];
                    s[a[x]/j]+=sgn*1;
                }
            }
        t[x]^=1;
        printf("%lld
",ans);
    }
}
int sb=haha();
int main(){;}