Clarke and math——狄利克雷卷积

行走天涯的豆沙包

于 2020-03-31 18:37:24 发布

阅读量159

点赞数

分类专栏：数学

本文链接：https://blog.csdn.net/weixin_42979819/article/details/105229014

版权

数学专栏收录该内容

97 篇文章 0 订阅

订阅专栏

博客探讨了狄利克雷卷积的概念，通过举例解释了当k分别为1和2时，函数g与f的关系，指出g=f∗(I)k，其中I表示单位函数。

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

传送门

题解：
当k=1的时候 $g(i)=∑_{i1|i}f(i1)$ 写成狄利克雷卷积形式为 $g = f * I$
k=2 $g = (f * I) * I$
以此类推： $g=f*(I)^k$

#include <bits/stdc++.h>
using namespace std;
#define  LL long long
const int N=1e5+7;
int a[N],tmp[N],g[N],I[N],n,k;
const int mod = 1e9 + 7;
void mul(int *a,int *b)
{
    for(int i=1;i<=n;i++) tmp[i]=0;
    for(int i=1;i<=n;i++){
        for(int j=1;j*i<=n;j++){
            tmp[i*j]=(tmp[i*j]+1ll*a[i]*b[j])%mod;
        }
    }
    for(int i=1;i<=n;i++) a[i]=tmp[i];
}
int main()
{
    int t; scanf("%d",&t);
    while(t--){
        scanf("%d%d",&n,&k);
        for(int i=1;i<=n;i++) scanf("%d",&g[i]),I[i]=1;
        while(k){
            if(k&1) mul(g,I);
            mul(I,I);
            k>>=1;
        }
        for(int i=1;i<=n;i++) printf("%d%c",g[i],i==n?'\n':' ');
    }
}