hdu5628 Clarke and math（狄利克雷卷积快速幂）

最新推荐文章于 2019-08-24 13:24:48 发布

SC.ldxcaicai

最新推荐文章于 2019-08-24 13:24:48 发布

阅读量382

点赞数 2

分类专栏： # 数论 # 数学

本文链接：https://blog.csdn.net/dreaming__ldx/article/details/98186147

版权

博客介绍了如何利用狄利克雷卷积快速幂解决HDU5628问题，指出常规方法的时间复杂度为O(nlogn^2)，并提出一种更优的解决方案，将单次卷积复杂度降低到O(nn*lnn)，最终总时间复杂度达到O(nlogn)。

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

传送门
一眼题。。。
注意不能有行末空格233
分析一下显然有 $g=I^k*f$
直接上狄利克雷卷积+快速幂即可。。。
复杂度 $O(n\log n^2)$
代码：

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int rlen=1<<18|1;
inline char gc(){
   
    static char buf[rlen],*ib,*ob;
    (ib==ob)&&(ob=(ib=buf)+fread(buf,1,rlen,stdin));
    return ib==ob?-1:*ib++;
}
inline int read(){
   
    int ans=0;
    char ch=gc();
    while(!isdigit(ch))ch=gc();
    while(isdigit(ch))ans=((ans<<2)+ans<<1)+(ch^48),ch=gc();
    return ans;
}
typedef long long ll;
const int mod=1e9+7;
inline int add(const int&a,const int&b){
   return a+b>=mod?a+b-mod:a+b;}
inline int dec(const int&a,const int&b){
   return a>=b?a-b:a-b+mod;}
inline int mul(const int&a,const int&b){
   return (ll)a*b%mod;}
inline void Add(int&a,const int&b){
   a=a+b>=mod?a+b-mod:a+b;}
inline void Dec(int&a,const int&b){
   a=a>=b?a-b:a-b+mod;}
inline void Mul(int&a,const int&b){
   a=(ll)a*b%mod;}
inline int ksm(int a,int p){
   int ret=1;for(;p;p>>=1,a=mul(a,a))if(p&1)Mul(ret,a);return ret;}\
const int N=1e5+5;
struct poly{
   
    vector<int>a;
    poly