一个人的数论

2021-08-28 07:31:07 阅读：198 来源： 互联网

莫比乌斯反演神题

换言说，就算你推出柿子还是不会那个神奇的伯努力数。。。。

那么这样，我们先推一波柿子：

$f_d(n)
\\=\sum_{i=1}^{n}[gcd(i,n)=1]i^d
\\ =\sum_{i=1}^{n}i^d\sum_{t|gcd(i,n)}\mu (t)
\\ =\sum_{i=1}^{n}i^d\sum_{t|i \textit{ } and \textit{ } t|n}\mu(t)
\\ =\sum_{t|n}\mu(t)\sum_{i=1}^{\frac{n}{t}} (it)^d
\\ =\sum_{t|n}\mu(t)t^d\sum_{i=1}^{\frac{n}{t}}i^d$

发现柿子里后面的$\sum_{i=1}^{\frac{n}{t}}i^d$极为类似伯努力数的递推公式：

$S^m(n)=\frac{1}{m+1}\sum_{i=0}^{m}C_{m+1}^{i}B_in^{m+1-i}$

其中$S^m(n)=\sum_{i=0}^{n-1}i^d$

那么我们可以知道题目中的关于这一部分的答案。

这也就是伯努力数的应用，将一个看似不可循环的东西（比如说这题里面的$n$）变成了可以循环的东西（这道题里面用伯努力数把循环复杂度降到$d$）

然后题目就会变得可做，详细证明和代码实现可参考$oi-wiki$

我们先把后面的柿子写成多项式的形式代回到原式（里面的$a_i$就是多项式系数：$\frac{1}{m+1}\sum_{i=0}^{m}C_{m+1}^{i}B_i$）：

$\\=\sum_{t|n}\mu(t)t^d\sum_{i=0}^{d+1}a^i(\frac{n}{t})^i
\\=\sum_{i=0}^{d+1}a_i\sum_{t|n}\mu(t)t^d(\frac{n}{t})^i
\\=\sum_{i=0}^{d+1}a_i\sum_{t|n}\mu(t)N^it^{d-i}
\\=\sum_{i=0}^{d+1}a_iN^i\sum_{t|n}\mu(t)t^{d-i}$

这时我们发现后面的一部分$\sum_{t|n}\mu(t)t^{d-i}$是一个积性函数

设$g_n=\sum_{t|n}\mu(t)t^{d-i}=\sum_{t|n}\mu(t)id_{d-i}(t)$，而$\mu$和$id$都是积性函数，那么他们相乘也是积性函数

我看$oi-wiki$上面只说了积性函数卷积的积性，但没说这个东西所以以为自己错了，其实很显然，用乘法结合律和分配律就可以简单证明

就这样我们发现$g_n$是积性函数，那么$g_n=\prod_{j=1}^{\omega} g({p_j}^{c_j})$

题目里面给的是质因数分解，这个不用我们搞了，那么考虑贡献

拿出一个$\prod$里面的项来说事：

$g(p^c)=\sum_{t|p^c}\mu(t)id_{d-i}(t)$

发现莫比乌斯函数很棒，只有在$c$为$0/1$的时候整个柿子的答案不是$0$，那么

$\\g(p^c)=\sum_{t|p^c}\mu(t)id_{d-i}(t)
\\=\mu(1)id_{d-i}(1)+\mu(p)id_{d-i}(p)
\\=1-p^{d-i}$

这样，整个柿子的化简工作就大功告成了！

$f_d(n)=\sum_{i=0}^{d+1}a_iN^i\prod_{j=1}^{\omega}(1-p_j^{d-i})$

$a$可以直接筛出，后面的也可以直接搞，都是超级好实现的那种，

数据范围还贼小，比较舒服的$O(d^2+d\omega)$

 1 #include<bits/stdc++.h>
 2 #define int long long
 3 using namespace std;
 4 namespace AE86{
 5     inline int read(){
 6         int x=0,f=1;char ch=getchar();
 7         while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
 8         while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return x*f;
 9     }inline void write(int x,char opt='\n'){
10         char ch[20];int len=0;if(x<0)x=~x+1,putchar('-');
11         do{ch[len++]=x%10+(1<<5)+(1<<4);x/=10;}while(x);
12         for(int i=len-1;i>=0;--i)putchar(ch[i]);putchar(opt);}
13 }using namespace AE86;
14 const int mod=1e9+7;
15 int d,w,p[1001],t[1001],N=1,ans;
16 int C[150][150],inv[150],B[150],xi[150];
17 inline int qmo(int a,int b){int ans=1,c=mod;a%=c;while(b){if(b&1)ans=ans*a%c;b>>=1;a=a*a%c;}return ans;}
18 inline void getB(){
19     for(int i=0;i<150;i++){C[i][0]=C[i][i]=1;for(int j=1;j<i;j++)C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;}
20     inv[1]=1;for(int i=2;i<150;i++)inv[i]=(mod-mod/i)*inv[mod%i]%mod;
21     B[0]=1;
22     for(int i=1;i<150;i++){
23         int ans=0;if(i==150-1) break;
24         for(int j=0;j<i;j++) (ans+=C[i+1][j]*B[j]%mod)%=mod;
25         ans=(ans*(-inv[i+1])%mod+mod)%mod; B[i]=ans;
26     } for(int i=0;i<=d;i++) xi[d+1-i]=qmo(d+1,mod-2)*C[d+1][i]%mod*B[i]%mod;
27 }
28 namespace WSN{
29     inline short main(){
30         d=read(); w=read(); getB();
31         for(int i=1;i<=w;i++) p[i]=read(),t[i]=read(),N=N*qmo(p[i],t[i])%mod;
32         for(int i=1;i<=d+1;i++){
33             int now=xi[i]*qmo(N,i)%mod,tmp=1;
34             for(int j=1;j<=w;j++) tmp=tmp*(1-qmo(p[j],d-i+mod-1)+mod)%mod;
35             (ans+=now*tmp%mod)%=mod;
36         } write(ans);
37         return 0;
38     }
39 }
40 signed main(){return WSN::main();}

一个人的数论

标签：ch,frac,数论,sum,mu,一个,int,积性
来源： https://www.cnblogs.com/hzoi-wsn/p/15195837.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

一个人的数论