如题,要求一个式子∑d|n ψ(x)*n/d,这题是个公式推导,然而我连数论是什么都不知道,看别人的推导都是一脸懵,
在网上查了半天定义,性质。。
这个式子是一个狄利克雷卷积出来的
然后根据积性函数的笛利克雷卷积还是积性函数。
积性函数有个性质是f(ab)=f(a)f(b),(a,b互质)
这样把n=(p1^q1)*(p2^q2)*...........(pm^qm)
拆开,因为pi和pj是互质的,所以
ans = (∑d|n1 ψ(x)*n1/d)*(∑d|n2 ψ(x)*n2/d)n2*........
随便拿出一个nk,欧拉函数还有个性质:若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。
就这样,最后结果乘起来就行了
代码,转载自点击打开链接
#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
const int N=10000+10;
#define MOD 998244353
ll power(ll a,ll b)
{
ll res=1;
while(b)
{
if(b&1)
res=(res*a)%MOD;
b>>=1;
a=(a*a)%MOD;
}
return res%MOD;
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
ll p[100],q[100];
int m;
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
scanf("%lld%lld",p+i,q+i);
}
ll ans=1;
for(int i=1;i<=m;i++)
{
ll mu=power(p[i],q[i]-1);
ans=(q[i]*p[i]%MOD+p[i]-q[i]+MOD)*mu%MOD*ans%MOD;
}
printf("%lld\n",ans);
}
return 0;
}