[hdu5628]Clarke and math（dirichlet卷积）

传送门
狄利克雷卷积定义：\[(f*g)(n)=\sum_{d|n}f(d)*g({\frac n d})\]狄利克雷卷积满足交换律：\[f*g=g*f\]结合律：\[(f*g)*h=f*(g*h)\]还有这么几个性质：\[f*\varepsilon=f\] \[f*1=\sum_{d|n}f(d)\]其中\[1(n)=1,\varepsilon(n)=[n=1]\]我们做这个题就是用的上面那条，题目中的式子可以化成这样：\[1^k*f\]然后快速幂就好了。
代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){int x=0;char ch=' ';int f=1;while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();if(ch=='-')f=-1,ch=getchar();while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();return x*f;
}
const int p=1e9+7;
int n,k;
ll a[100001];
void dirichlet(ll *ans,ll *f){memset(a,0,sizeof(a));for(int i=1;i*i<=n;i++){a[i*i]=(a[i*i]+ans[i]*f[i])%p;for(int j=i+1;j*i<=n;j++){a[i*j]=(a[i*j]+ans[i]*f[j])%p;a[i*j]=(a[i*j]+ans[j]*f[i])%p;}}for(int i=1;i<=n;i++)ans[i]=a[i];
}
ll f[100001],ans[100001],x[100001];
void ksm(){int b=k;while(b){if(b&1){dirichlet(ans,x);}dirichlet(x,x);b>>=1;}
}
void solve(){for(int i=1;i<=n;i++){f[i]=read();x[i]=1;ans[i]=0;}ans[1]=1;ksm();dirichlet(f,ans);for(int i=1;i<n;i++){printf("%lld ",f[i]);}printf("%lld\n",f[n]);
}int T;
int main(){T=read();while(T--){n=read();k=read();solve();}return 0;
}

转载于:https://www.cnblogs.com/stone41123/p/7605431.html

[hdu5628]Clarke and math（dirichlet卷积）相关推荐

积性函数与Dirichlet卷积学习小记
前言首先感谢 XHM 大佬的悉心指导,我懂得了不少~. 链一下他关于这方面的见解.博客--XHM 的Dirichlet卷积学习小记一些定义回归正题,这次我学习了一下狄利克雷卷积方面的知识. 先 ...
NOI数学：狄利克雷（Dirichlet)卷积
从零开始的数学12:Dirichlet卷积(1) 从零开始的数学12:Dirichlet卷积(1)_哔哩哔哩_bilibili 从零开始的数学13:Dirichlet卷积(2) 从零开始的数学13:D ...
Dirichlet卷积学习小记
定义定义数论函数fff和g" role="presentation" style="position: relative;">ggg的狄利克雷 ...
Dirichlet卷积
莫比乌斯反演的前置知识定义设f,gf,gf,g是数论函数,考虑数论函数hhh满足 h(n)=∑d∣nf(d)g(nd)h(n)=\sum_{d|n}f(d)g(\frac{n}{d})h(n)=d ...
基础省选＋NOI 第5部分数论进阶
1.欧拉函数欧拉函数欧拉函数入门_哔哩哔哩_bilibili 听陈老师讲数学家的故事,放羊娃欧拉_哔哩哔哩_bilibili RSA 算法如何计算公钥和私钥_哔哩哔哩_bilibili SWPU- ...
python神经网络案例——CNN卷积神经网络实现mnist手写体识别
分享一个朋友的人工智能教程.零基础!通俗易懂!风趣幽默!还带黄段子!大家可以看看是否对自己有帮助:点击打开全栈工程师开发手册 (作者:栾鹏) python教程全解 CNN卷积神经网络的理论教程参考 ...
【学习小记】狄利克雷卷积+杜教筛
Preface 这东西分明就是玄学暴力用来求简单的数论函数的前缀和,像φ,μφ,\mu这类的东西当然,约数和,约数个数之类的也是可以的 Text 数论函数是指定义域是整数,陪域是复数的函数 Dir ...
《算法竞赛中的初等数论》（三）正文 0x30 积性函数（ACM / OI / MO）（十五万字符数论书）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划写在最前面:本文部分内容来自网上各大博客或是各类图书,由我个人整理,增加些许见解,仅做学习交流使用,无 ...
莫比乌斯，欧拉函数题目练习（完结）
Starttime:2020/11/16Start\ time:2020/11/16Start time:2020/11/16 Lastupdatetime:2020/11/28Last\ updat ...
数论函数 - 莫比乌斯函数与莫比乌斯反演 - 基础杜教筛
原文链接http://www.cnblogs.com/zhouzhendong/p/8627380.html 省选后发现我数学好差.于是先从数论开始学习. 如果发现本文有任何错误,欢迎留言指正. 本文 ...

[hdu5628]Clarke and math（dirichlet卷积）

[hdu5628]Clarke and math（dirichlet卷积）相关推荐

最新文章

热门文章