BZOJ 4332 FFT+快速幂
思路:
最裸的方程:f[i][j]=Σf[i-1][j-k]*F[k]
诶呦 这不是卷积嘛
f[i]就可以用f[i-1]卷F 求到
但是这样还是很慢
设p[i] 为Σ f[j](1<=j<=i)
发现p可以倍增推
于是 就 倍增一下 就完了...
http://www.cnblogs.com/Skyminer/p/6561689.html
hz神犇的题解写得非常详细..
//By SiriusRen
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=;
const double pi=acos(-);
int n,M,P,A,O,S,U,L,R[N],F[N],g[N],p[N],t[N],ans;
struct cplxd{double x,y;cplxd(){}cplxd(double X,double Y){x=X,y=Y;}}ca[N],cb[N],cc[N];
cplxd operator+(cplxd a,cplxd b){return cplxd(a.x+b.x,a.y+b.y);}
cplxd operator-(cplxd a,cplxd b){return cplxd(a.x-b.x,a.y-b.y);}
cplxd operator*(cplxd a,cplxd b){return cplxd(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
cplxd operator/(cplxd a,int b){return cplxd(a.x/b,a.y/b);}
void FFT(cplxd *a,int f){
for(int i=;i<n;i++)if(i<R[i])swap(a[i],a[R[i]]);
for(int i=;i<n;i<<=){
cplxd wn=cplxd(cos(pi/i),f*sin(pi/i));
for(int j=;j<n;j+=(i<<)){
cplxd w(,);
for(int k=;k<i;k++,w=w*wn){
cplxd x=a[j+k],y=w*a[j+k+i];
a[j+k]=x+y,a[j+k+i]=x-y;
}
}
}
if(f==-)for(int i=;i<n;i++)a[i]=a[i]/n;
}
void Pw(int *a,int *b,int *c){
for(int i=;i<n;i++)ca[i]=cplxd(a[i],);
for(int i=;i<n;i++)cb[i]=cplxd(b[i],);
FFT(ca,),FFT(cb,);
for(int i=;i<n;i++)cc[i]=ca[i]*cb[i];
FFT(cc,-);
for(int i=;i<=M;i++)c[i]=((int)(0.3+cc[i].x))%P;
}
void pow(int k){
if(k==)return;
pow(k>>);
Pw(p,g,t),Pw(g,g,g);
for(int i=;i<=M;i++)(p[i]+=t[i])%=P;
if(k&){
Pw(g,F,g);
for(int i=;i<=M;i++)(p[i]+=g[i])%=P;
}
}
signed main(){
scanf("%d%d%d%d%d%d",&M,&P,&A,&O,&S,&U);
for(n=;n<=M*;n<<=)L++;
for(int i=;i<=n;i++)R[i]=(R[i>>]>>)|((i&)<<(L-));
for(int i=;i<=M;i++)p[i]=g[i]=F[i]=(i*i*O+S*i+U)%P;
pow(A),printf("%d\n",p[M]);
}
BZOJ 4332 FFT+快速幂的更多相关文章
- BZOJ 3160: 万径人踪灭 FFT+快速幂+manacher
BZOJ 3160: 万径人踪灭 题目传送门 [题目大意] 给定一个长度为n的01串,求有多少个回文子序列? 回文子序列是指从原串中找出任意个,使得构成一个回文串,并且位置也是沿某一对称轴对称. 假如 ...
- 【BZOJ 4332】 4332: JSOI2012 分零食 (FFT+快速幂)
4332: JSOI2012 分零食 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 119 Solved: 66 Description 这里是欢乐 ...
- bzoj 4332 FFT型的快速幂(需要强有力的推导公式能力)
有n个小朋友,m颗糖,你要把所有糖果分给这些小朋友. 规则第 i 个小朋友没有糖果,那么他之后的小朋友都没有糖果..如果一个小朋友分到了 xx 个糖果,那么的他的权值是 f(x) = ox^2 + ...
- 【codeforces 623E】dp+FFT+快速幂
题目大意:用$[1,2^k-1]$之间的证书构造一个长度为$n$的序列$a_i$,令$b_i=a_1\ or\ a_2\ or\ ...\ or a_i$,问使得b序列严格递增的方案数,答案对$10^ ...
- P3321 [SDOI2015]序列统计 FFT+快速幂+原根
\(\color{#0066ff}{ 题目描述 }\) 小C有一个集合S,里面的元素都是小于M的非负整数.他用程序编写了一个数列生成器,可以生成一个长度为N的数列,数列中的每个数都属于集合S.小C用这 ...
- BZOJ 2179: FFT快速傅立叶
2179: FFT快速傅立叶 Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 2923 Solved: 1498[Submit][Status][Di ...
- bzoj(矩阵快速幂)
题意:定义Concatenate(1,N)=1234567……n.比如Concatenate(1,13)=12345678910111213.给定n和m,求Concatenate(1,n)%m. (1 ...
- bzoj 4000 矩阵快速幂优化DP
建立矩阵,跑快速幂 /************************************************************** Problem: 4000 User: idy002 ...
- bzoj 2326 矩阵快速幂
思路:矩阵快速幂搞一搞. #include<bits/stdc++.h> #define LL long long #define fi first #define se second # ...
随机推荐
- Git学习总结二(版本回退)
修改修改仓库中readme.txt文件,然后用git status命令看看结果: $ git status On branch master Changes not staged for commit ...
- python包与模块
Python基础-包与模块 摘要 为重用以及更好的维护代码,Python使用了模块与包:一个Python文件就是一个模块,包是组织模块的特殊目录(包含__init__.py文件). 模块搜索路径,Py ...
- =、==、is、id(内容)
= 赋值 == 比较值是否相等 is 比较.比较的是内存地址 id(内容) 测出内存地址
- 渗透实战(周六):Hydra&Metasploit暴力破解SSH登录口令
一. SSH服务开启前基础配置 1.1 修改配置文件
- Huawei-R&S-网络工程师实验笔记20190530-FTP上传下载、STelnet登录、SFTP登录
>Huawei-R&S-网络工程师实验笔记20190530-FTP上传下载.STelnet登录.SFTP登录 >>实验开始,参考<Huawei-R&S-网络工程 ...
- [HDU3518]Boring counting(后缀数组)
传送门 求出现超过1次的不重叠子串的个数 根据论文中的方法. 枚举子串的长度 k. 用 k 给 height 数组分组,每一组求解,看看当前组的位置最靠后的后缀和位置最靠前的后缀所差个数是否大于长度, ...
- win7安装gmpy2
1.下载地址:https://pypi.python.org/pypi/gmpy2 2.安装python和pip python 安装 下载: https://www.python.org/getit/ ...
- hdu 4971 多校10最大权闭合图
/* 很明显的最大权闭合图题 */ #include<stdio.h> #include<string.h> #include<queue> using names ...
- libcloud代码研究(二)——云服务封装
对Ali OSS和百度云存储的封装 先在libcloud.storage.provider文件中添加Ali OSS provider和Baidu provider,在DRIVERS中添加如下 ...
- python在Linux中安装虚拟环境,区别python2和python3,分别安装模块
安装虚拟环境的时候遇到的问题,解决的过程很闹心,在这里简单直接的分享出来,就是为了解决问题. 安装虚拟环境(须在联网状态下) $ sudo pip install virtualenv $ sudo ...