华华给月月出题---线性筛优化

最新推荐文章于 2025-12-22 09:47:08 发布

原创最新推荐文章于 2025-12-22 09:47:08 发布 · 289 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#算法

数学专栏收录该内容

15 篇文章

订阅专栏

这段代码实现了一个计算质数及其乘积的程序。它首先定义了一个快速幂运算函数，然后通过遍历并标记合数来找出质数。接着，计算每个质数的模逆并乘积，最后得到一个异或结果。程序适用于大整数处理，具有较高的计算效率。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
#define int long long
const int N = 15000000,mod = 1e9+7;
int qmi(int a,int b){int res=1;while(b){if(b&1)res=res*a%mod;a=a*a%mod;b>>=1;}return res%mod;}
int s[N];
bool st[N]={0};
int primes[1000000],idx=0;
signed main()
{
    int n;cin>>n;
    int res=1;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])primes[idx++]=i,s[i]=qmi(i,n);
        for(int j=0;primes[j]*i<=n;j++)
        {
            st[primes[j]*i]=true;
            s[primes[j]*i]=s[i]*s[primes[j]]%mod;
            if(i%primes[j]==0)break;
        }
    }
    for(int i=1;i<=n;i++)res^=s[i];
	cout<<res<<'\n';
}