题意:bc round 72 中文题面

分析(官方题解):

如果学过Dirichlet卷积的话知道这玩意就是g(n)=(f*1^k)(n),

由于有结合律,所以我们快速幂一下1^k就行了。

当然,强行正面刚和式也是能搞的(反正我不会)。

一次Dirichlet卷积复杂度是O(nlogn)的,所以总时间复杂度为O(nlognlogk).

注:普及Dirichlet卷积概念,设f,g为两个数论函数,

那么规定 (f∗g)=∑d|nf(d)g(n/d)为f,g的Dirichlet卷积

Dirichlet卷积有以下性质,有交换律和结合律,对加法有分配律
f∗(g∗h)=(f∗g)∗h
f∗(g+h)=f∗g+f∗h
f∗g=g∗f

然后对于这个题解上所说的一次Dirichlet卷积是O(nlogn)这个我不是很懂

但是看网上大神写的估计也差不多,感觉和埃氏筛的复杂度差不多甚至还要小应该是在O(nlognlogn)左右

所以下面的代码的复杂度大概应该是O(nlog^2nlogk)

然后顺带提一下,对于Dirichlet卷积来说,规定有一个幺元e的话

即:f*e=f,

这个幺元函数有一个性质e[1]=1,e[k]=0,k=1,2,3.....

这个在快速幂的时候会用到

然后上代码

#include <iostream>
#include <cstdio>
using namespace std;
typedef long long LL;
const int N=1e5+;
const LL mod=1e9+;
int n,k;
LL f[N],s[N],ans[N],tmp[N];
void solve()
{
while(k)
{
if(k&)
{
for(int i=; i<=n; ++i)tmp[i]=;
for(int i=; i*i<=n; ++i)
{
tmp[i*i]=(tmp[i*i]+ans[i]*s[i])%mod;
for(int j=i+; j*i<=n; ++j)
tmp[i*j]+=(ans[i]*s[j]%mod+ans[j]*s[i]%mod)%mod,tmp[i*j]%=mod;
}
for(int i=; i<=n; ++i)ans[i]=tmp[i];
}
for(int i=; i<=n; ++i)tmp[i]=;
for(int i=; i*i<=n; ++i)
{
tmp[i*i]=(tmp[i*i]+s[i]*s[i])%mod;
for(int j=i+; j*i<=n; ++j)
tmp[i*j]+=(s[i]*s[j]%mod+s[j]*s[i]%mod)%mod,tmp[i*j]%=mod;
}
for(int i=;i<=n;++i)s[i]=tmp[i];
k>>=;
}
for(int i=; i<=n; ++i)tmp[i]=;
for(int i=; i*i<=n; ++i)
{
tmp[i*i]=(tmp[i*i]+f[i]*ans[i])%mod;
for(int j=i+; j*i<=n; ++j)
tmp[i*j]+=(f[i]*ans[j]%mod+f[j]*ans[i]%mod)%mod,tmp[i*j]%=mod;
}
for(int i=;i<=n;++i)ans[i]=tmp[i];
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&k);
for(int i=; i<=n; ++i)
scanf("%I64d",&f[i]),s[i] = ,ans[i] = ;
ans[] = ;
solve();
for(int i=; i<n; ++i)
printf("%I64d ",ans[i]);
printf("%I64d\n",ans[n]);
}
return ;
}

HDU 5628 Clarke and math Dirichlet卷积+快速幂的更多相关文章

  1. HDU.5628.Clarke and math(狄利克雷卷积 快速幂)

    \(Description\) \[g(i)=\sum_{i_1|i}\sum_{i_2|i_1}\sum_{i_3|i_2}\cdots\sum_{i_k|i_{k-1}}f(i_k)\ mod\ ...

  2. HDU 5628 Clarke and math——卷积,dp,组合

    HDU 5628 Clarke and math 本文属于一个总结了一堆做法的玩意...... 题目 简单的一个式子:给定$n,k,f(i)$,求 然后数据范围不重要,重要的是如何优化这个做法. 这个 ...

  3. HDU 5628 Clarke and math dp+数学

    Clarke and math 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5628 Description Clarke is a patient ...

  4. HDU 1757 A Simple Math Problem (矩阵快速幂)

    题目 A Simple Math Problem 解析 矩阵快速幂模板题 构造矩阵 \[\begin{bmatrix}a_0&a_1&a_2&a_3&a_4&a ...

  5. hdu 1757 A Simple Math Problem_矩阵快速幂

    题意:略 简单的矩阵快速幂就行了 #include <iostream> #include <cstdio> #include <cstring> using na ...

  6. hdu 5187 zhx's contest [ 找规律 + 快速幂 + 快速乘法 || Java ]

    传送门 zhx's contest Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others ...

  7. HDU 2855 斐波那契+矩阵快速幂

    http://acm.hdu.edu.cn/showproblem.php?pid=2855 化简这个公式,多写出几组就会发现规律 d[n]=F[2*n] 后面的任务就是矩阵快速幂拍一个斐波那契模板出 ...

  8. HDU 5950:Recursive sequence(矩阵快速幂)

    http://acm.hdu.edu.cn/showproblem.php?pid=5950 题意:给出 a,b,n,递推出 f(n) = f(n-1) + f(n-2) * 2 + n ^ 4. f ...

  9. HDU1757 A Simple Math Problem 矩阵快速幂

    A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Ot ...

随机推荐

  1. Xcode文档下载与安装路径

    https://developer.apple.com/library/downloads/docset-index.dvtdownloadableindex ~/Library/Developer/ ...

  2. multimap和multiset 认知和使用

    之前只是在C++ Primer里面看过关联容器,可能因为没有实际用过,只是看看,所以导致用的时候并不熟悉: 在这之前,map和set的特性应该要了解,map是关联数组,也就是由键值对组成的,而set只 ...

  3. js解决checkbox全选和反选的问题

    function SelectAll() { var checkboxs=document.getElementsByName("chk_list"); for (var i=0; ...

  4. linux fork函数与vfork函数

    一.fork1. 调用方法#include <sys/types.h>#include <unistd.h> pid_t fork(void);正确返回:在父进程中返回子进程的 ...

  5. 【leetcode】4Sum(middle)

    Given an array S of n integers, are there elements a, b, c, and d in S such that a + b + c + d = tar ...

  6. 在线学习SQL语句?没问题~~

    以前弄得少,没注意.. http://sqlfiddle.com/ CREATE TABLE Presidents ( Id INT UNSIGNED NOT NULL AUTO_INCREMENT, ...

  7. HDU4628+状态压缩DP

    /* 状态压缩DP dp[ i ]:达到i状态的最小step. 题意:每次可以去掉一个回文串,求最少几步能取完. */ #include<stdio.h> #include<stri ...

  8. Redis hash数据类型操作

    Redis hash是一个string类型的field和value的映射表.一个key可对应多个field,一个field对应一个value.将一个对象存储 为hash类型,较于每个字段都存储成str ...

  9. fork产生子进程利用pipe管道通信

    http://siqun.blog.163.com/blog/static/213496001201341231121720/ 转载链接:http://hi.baidu.com/hj11yc/item ...

  10. [java线段树]2015上海邀请赛 D Doom

    题意:n个数 m个询问 每个询问[l, r]的和, 再把[l, r]之间所有的数变为平方(模为9223372034707292160LL) 很明显的线段树 看到这个模(LLONG_MAX为922337 ...