P5493 【模板】质数前缀统计

最新推荐文章于 2023-04-05 12:04:16 发布

_Shmily

最新推荐文章于 2023-04-05 12:04:16 发布

阅读量191

点赞数

分类专栏：数论基本算法

本文链接：https://blog.csdn.net/zhaoxinxin1234/article/details/104297724

版权

该博客介绍了如何使用拉格朗日插值法计算质数前缀统计，并强调了在处理过程中对质数贡献进行min25筛法优化以及预处理的重要性，以降低算法的时间复杂度。

摘要由CSDN通过智能技术生成

这题需要用拉格朗日插值法求前n项的 i^k 的和。
然后用min25筛一下质数的贡献。
有点卡常，能预处理的都预处理。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#define ll long long
#define llu unsigned ll
using namespace std;
const int maxn=1001000;
ll mod;
ll pre[maxn],suf[maxn],fac[maxn],inv[maxn],a[maxn];
ll y[maxn];
ll inv1;

ll p[maxn],id1[maxn],id2[maxn],pc[maxn];
ll g[maxn],sum[maxn],wi[maxn];
ll cnt,sq,n,k;
bool ha[maxn];

ll mypow(ll a,ll b)
{
   
    ll ans=1;
    while(b)
    {
   
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}

void init(ll k)
{
   
    int m=k+2;
    inv1=mypow(mod-1,mod-2