HDU 5628 Clarke and math Dirichlet卷积+快速幂
题意: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卷积+快速幂的更多相关文章
- 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\ ...
- HDU 5628 Clarke and math——卷积,dp,组合
HDU 5628 Clarke and math 本文属于一个总结了一堆做法的玩意...... 题目 简单的一个式子:给定$n,k,f(i)$,求 然后数据范围不重要,重要的是如何优化这个做法. 这个 ...
- HDU 5628 Clarke and math dp+数学
Clarke and math 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5628 Description Clarke is a patient ...
- HDU 1757 A Simple Math Problem (矩阵快速幂)
题目 A Simple Math Problem 解析 矩阵快速幂模板题 构造矩阵 \[\begin{bmatrix}a_0&a_1&a_2&a_3&a_4&a ...
- hdu 1757 A Simple Math Problem_矩阵快速幂
题意:略 简单的矩阵快速幂就行了 #include <iostream> #include <cstdio> #include <cstring> using na ...
- hdu 5187 zhx's contest [ 找规律 + 快速幂 + 快速乘法 || Java ]
传送门 zhx's contest Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others ...
- HDU 2855 斐波那契+矩阵快速幂
http://acm.hdu.edu.cn/showproblem.php?pid=2855 化简这个公式,多写出几组就会发现规律 d[n]=F[2*n] 后面的任务就是矩阵快速幂拍一个斐波那契模板出 ...
- 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 ...
- HDU1757 A Simple Math Problem 矩阵快速幂
A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Ot ...
随机推荐
- 第一个js面向对象的小实验
$.extend({ cal: function (num1,num2,oper,aftercal) { this.n1 = num1; ...
- 类似nike+、香蕉打卡的转场动画效果-b
首先,支持并感谢@wazrx 的 http://www.jianshu.com/p/45434f73019e和@onevcat 的https://onevcat.com/2013/10/vc-tran ...
- HTML5 中的块级链接
英文叫做 “Block-level” links,我以为只有我厂那些鸟毛不知道,没想到不知道的还挺多, 需要普及一下. 最近看了 kejun 的 PPT 前端开发理论热点面对面:从怎么看,到怎么做?, ...
- java对象数组
问题描述: java 对象数组的使用 问题解决: 数组元素可以是任何类型(只要所有元素具有相同的类型) 数组元素可以是基本数据类型 数组元素也可以是类对象,称这样的数组为对象数组.在这种情况下 ...
- [Oracle]any, all解析
因为很少用到, 所以几乎忘记了这几个函数, 不过它们还是很有用的使用它们可以大大简化一些SQL文的语法, 至于效率问题, 如CCW所说它们和EXISTS, IN 之类没有什么差别, 而且要具体问题具体 ...
- html5移动web开发实战必读书记
原文 http://itindex.net/detail/50689-html5-移动-web 主题 HTML5 一.配置移动开发环境 1.各种仿真器.模拟器的下载安装 http://www.mob ...
- C++11新特性:右值引用和转移构造函数
问题背景 #include <iostream> using namespace std; vector<int> doubleValues (const vector< ...
- 删除appcompat_v7会出很多错误
创建工程的时候会出现appcompat_v7这个文件夹 手贱删除后,发现出错了 说明test项目是依赖于appcompat_v7包的,所以这个appcompat_v7包是不能被删除的. appcomp ...
- ArcGIS学习记录—KMZ KML与SHP文件互相转换
1.在google earth中绘制边界 工具栏中选择"Add Polygon".随意绘制一个多边形. 右击添加的图层名(左侧)保存位置为,选择保存为kmz或kml文件. ...
- 理解extern char s[100]与extern char *s
在x.c中定义了一个字符数组 char s[100],在l.c中进行引用extern char s[200], 有些c程序新手经常把它写成extern char *s. 这两种写法的含义一样吗? 首先 ...