YbtOJ NOIP2020 模拟赛 B 组 Day10 D. 挺好序列【数论】

---

题目：

传送门

---

题意：

数论题，乖乖去看题面吧

---

分析：

---

代码：

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define LL long long
#define mo 998244353
using namespace std;
inline LL read()
{
	LL s=0,f=1; char c=getchar();
	while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
	while(c>='0'&&c<='9') {s=s*10+c-'0';c=getchar();}
	return s*f;
}
LL n,m,f[20000005],ans,p[80],pri[20000005],cnt,B;
bool v[20000005];
LL power(LL x,LL y)
{
	LL s=1;
	while(y)
	{
		if(y&1) (s*=x)%=mo;
		(x*=x)%=mo; y>>=1;
	}
	return s;
}
int main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	n=read();m=read();B=read();
	f[1]=1;p[1]=1;
	for(LL i=1;i<=60;i++) p[i]=power(i,n);
	for(LL i=2;i<=m;i++)
	{
		if(!v[i])
		{
			LL d=0,P=1,last=1;
			pri[++cnt]=i; 
			while(!(B%i)) B/=i,d++; 
			for(LL j=i,c=1;j<=m;j*=i,c++)
			{
				v[j]=1;
				if(c<=d)
				{
					(f[j]=f[j/i]*i%mo+p[c+1]-p[c]+mo)%=mo;
					last=f[j];P=j;
				}
				else
				{
					f[j]=(last*i%mo+p[c+1]-p[c])%mo-(P*i)%mo*(p[c-d]-p[c-d-1])%mo;
					last=f[j];(f[j]+=P*p[c-d])%=mo;f[j]=(f[j]+mo)%mo; 
				}
			}
		}
		for(LL j=1;j<=cnt;j++)
		{
			if(!(i%pri[j])||i*pri[j]>m) break;
			for(LL k=pri[j];i*k<=m;k*=pri[j])
			  v[i*k]=1,f[i*k]=f[i]*f[k]%mo;
		}
	}
	for(LL i=1;i<=m;i++) (ans+=f[i])%=mo;
	printf("%lld",ans);
	return 0;
}