牛客小白月赛 14

I
积性函数 线性筛 ？？
积性函数的一个博客
积性函数

#include<bits/stdc++.h>
#define mod 1000000007
using namespace std;
int n,k,ans;
int f[10000010],p[10000010];
int prime[2000005],cnt=0;
int ksm(int x,int y)
{
    int res=1;
    while(y)
    {
        if(y&1)
            res=1LL*res*x%mod;
        x=1LL*x*x%mod;
       y>>=1;
    }
    return res;
}

int main()
{
    scanf("%d %d",&n,&k);
    f[1]=1;
    for(int i=2; i<=n; i++)
    {
        if(p[i]==0)
        {
            prime[++cnt]=i;
            f[i]=ksm(i,k);
        }
        for(int j=1; j<=cnt&&i*prime[j]<=n; j++)
        {
            p[i*prime[j]]=1;
            f[i*prime[j]]=1LL*f[i]*f[prime[j]]%mod;
            if(i%prime[j]==0)
                break;
        }
    }
    for(int i=1; i<=n; i++)
    {
        ans+=1LL*f[i]*(n/i)%mod;
        ans%=mod;
    }
    printf("%d\n",ans);

    return 0;
}

G题

#include<bits/stdc++.h>
using namespace std;
const int N=2e6+6;
int a[N];
long long b[N];

int main()
{
    int n,m;
    cin>>n>>a[1]>>m;
    long long ans=0;
    for (int i=2;i<=n;i++)
    {
        a[i]=(a[i-1]+7*i)%m;
        
    }
    for (int i=1;i<=n;i++)
    {
        for (int j=i;j<=n;j+=i)
            b[j]+=a[i];
        
    }
    ans=b[1];
    for (int i=2;i<=n;i++)
        ans=ans^b[i];
    
    cout<<ans<<endl;
    
    return 0;
    
}