CF915G Coprime Arrays

一、题目

点此看题

二、解法

设$f(i)$为值域在$[1,k]$的$n$个数，最后的$\gcd$为$i$的方案数。

设$F(i)$为值域在$[1,k]$的$n$个数，最后$\gcd$为$i$的倍数的方案数。
$$F(i)=\sum _{i|d}f(d)=(\frac{k}{i}{)}^n$$$$f(i)=\sum _{i|d}\mu (\frac{d}{i})F(d)$$$$f(1)=\sum _{i=1}^k\mu (i)F(i)=\sum _{i=1}^k\mu (i)(\frac{k}{i}{)}^n$$上面的式子显然可以分块，但是不够快，由于异或的存在，我们必须处理出每一个$b_k$，这时我们就要考虑它的共性，可以考虑循环中每一个$i$的贡献，把$k$划分成若干段，每一段内的贡献是一样的（每一段$k/i$一样），我们就可以打差分标记，注意要删除以前的标记，那么时间复杂度就是$O(n\log n)$

#include <cstdio>
#include <iostream>
using namespace std;
#define int long long
const int M = 2000005;
const int jzm = 1e9+7;
int read()
{
    int x=0,flag=1;
    char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int n,k,ans,cnt,pw[M],p[M],b[M],vis[M],mu[M];
int qkpow(int a,int b)
{
    int r=1;
    while(b>0)
    {
        if(b&1) r=r*a%jzm;
        a=a*a%jzm;
        b>>=1;
    }
    return r;
}
void init(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            p[++cnt]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=cnt && i*p[j]<=n;j++)
        {
            vis[i*p[j]]=1;
            if(i%p[j]==0) break;
            mu[i*p[j]]=-mu[i];
        }
    }
}
signed main()
{
    n=read();k=read();
    init(k);
    for(int i=1;i<=k;i++)
        pw[i]=qkpow(i,n);
    for(int i=1;i<=k;i++)
        for(int j=1;j*i<=k;j++)
        {
            b[i*j]=(b[i*j]+mu[i]*pw[j])%jzm;
            b[i*j]=(b[i*j]-mu[i]*pw[j-1])%jzm;
        }
    for(int i=1;i<=k;i++)
    {
        b[i]=((b[i]+b[i-1])%jzm+jzm)%jzm;
        ans=(ans+(b[i]^i)%jzm)%jzm;
    }
    printf("%lld\n",ans);
}