bzoj 4929: 第三题
Description
Input
Output
#include<bits/stdc++.h>
typedef std::vector<int> pol;
typedef long long i64;
typedef unsigned long long u64;
const int P=;
int max(int a,int b){return a>b?a:b;}
template<int p>
i64 pw(int a,int n){
if(n<)n+=p-;
i64 v=;
for(;n;n>>=,a=i64(a)*a%p)if(n&)v=v*a%p;
return v;
}
int N,K,r[<<|];
typedef double ld;
struct cplx{
ld a,b;
cplx(ld _x=,ld _y=):a(_x),b(_y){}
cplx operator+(cplx w){return cplx(a+w.a,b+w.b);}
cplx operator-(cplx w){return cplx(a-w.a,b-w.b);}
cplx operator*(cplx w){return cplx(a*w.a-b*w.b,a*w.b+b*w.a);}
}A[<<|],B[<<|],*Es[][],C[<<|],D[<<|];
void pre(int n){
for(N=,K=;N<n;N<<=,++K);
for(int i=;i<N;++i)r[i]=r[i>>]>>|(i&)<<K;
memset(A,,sizeof(cplx)*N);
memset(B,,sizeof(cplx)*N);
memset(C,,sizeof(cplx)*N);
memset(D,,sizeof(cplx)*N);
}
const ld pi=acos(-);
void dft(cplx*a,int t){
for(int i=;i<N;++i)if(i<r[i])std::swap(a[i],a[r[i]]);
for(int i=,z=;i<N;i<<=,++z){
if(!Es[t==][z]){
cplx*E=Es[t==][z]=new cplx[i];
for(int j=;j<i;++j)E[j]=cplx(cos(j*pi/i),t*sin(j*pi/i));
}
cplx*E=Es[t==][z];
for(int j=;j<N;j+=i<<){
cplx*b=a+j,*c=b+i;
for(int k=;k<i;++k){
cplx x=b[k],y=c[k]*E[k];
b[k]=x+y;
c[k]=x-y;
}
}
}
if(t==-)for(int i=;i<N;++i)a[i].a/=N;
}
void mul(cplx*a,cplx*b){
for(int i=;i<N;++i)a[i]=a[i]*b[i];
}
void mov(const pol&x,cplx*a){
for(int i=;i<x.size();++i)a[i].a=x[i];
}
void mov(cplx*a,const pol&_x){
pol&x=const_cast<pol&>(_x);
for(int i=;i<x.size();++i)x[i]=i64(a[i].a+0.49)%P;
}
void mov(const pol&x,int*a){
memcpy(a,x.data(),x.size()*sizeof(int));
}
void mov(int*a,const pol&_x){
pol&x=const_cast<pol&>(_x);
memcpy(x.data(),a,x.size()*sizeof(int));
}
void chk(pol&w){
int p=w.size();
while(p&&!w[p-])--p;
w.resize(p);
}
void rev(pol&w){
std::reverse(w.data(),w.data()+w.size());
}
u64 _c[];
pol operator*(const pol&a,const pol&b){
pol c(a.size()+b.size()-);
if(c.size()<=){
const int*as=a.data(),ap=a.size(),*bs=b.data(),bp=b.size();
int*cs=c.data();
for(int i=;i<c.size();++i)_c[i]=;
for(int i=;i<ap;++i)for(int j=;j<bp;++j)_c[i+j]+=u64(as[i])*bs[j];
for(int i=;i<c.size();++i)cs[i]=_c[i]%P;
}else if(c.size()<=){
pre(c.size());
mov(a,A),mov(b,B);
dft(A,),dft(B,);
mul(A,B);
dft(A,-);
mov(A,c);
}else{
pre(c.size());
for(int i=;i<a.size();++i){
A[i]=a[i]>>;
B[i]=a[i]&;
}
for(int i=;i<b.size();++i){
C[i]=b[i]>>;
D[i]=b[i]&;
}
dft(A,);dft(B,);dft(C,);dft(D,);
for(int i=;i<N;++i){
cplx v1=A[i]*C[i];
cplx v2=A[i]*D[i]+B[i]*C[i];
cplx v3=B[i]*D[i];
A[i]=v1;
B[i]=v2;
C[i]=v3;
}
dft(A,-);dft(B,-);dft(C,-);
for(int i=;i<c.size();++i){
c[i]=((i64(A[i].a+0.49)<<)+(i64(B[i].a+0.49)<<)+i64(C[i].a+0.49))%P;
}
}
chk(c);
return c;
}
pol operator*(const pol&w,int x){
pol a(w);
for(int i=;i<a.size();++i)a[i]=i64(a[i])*x%P;
return a;
}
pol operator-(const pol&a,const pol&b){
pol c(max(a.size(),b.size()));
mov(a,c.data());
for(int i=;i<b.size();++i)c[i]=(c[i]-b[i]+P)%P;
chk(c);
return c;
}
pol inv(pol a){
if(a.size()==){
a.resize();
a[]=pw<P>(a[],P-);
return a;
}
pol b=inv(pol(a.data(),a.data()+(a.size()-)/+));
pol c=b*b*a;
c=b*-c;
c.resize(a.size());
return c;
}
pol operator/(pol a,pol b){
chk(a),chk(b);
rev(a),rev(b);
int sz=a.size()-b.size()+;
b.resize(sz);
pol c=a*inv(b);
c.resize(sz);
rev(c);
chk(c);
return c;
}
pol operator%(const pol&a,const pol&b){
if(a.size()<b.size())return a;
return a-a/b*b;
}
int cal(const pol&w,i64 x){
int y=;
for(int p=w.size()-;p>=;--p)y=(y*x+w[p])%P;
return y;
}
int n,b,c,d,e,a[],as[];
struct Q{int x,id;}qs[];
bool operator<(Q a,Q b){return a.x<b.x;}
pol tr[];
void calc0(int w,int L,int R){
if(L==R){
tr[w]=pol();
tr[w][]=;
tr[w][]=(P-qs[L].x)%P;
return;
}
int M=(L+R)/;
calc0(w<<,L,M);
calc0(w<<^,M+,R);
tr[w]=tr[w<<]*tr[w<<^];
}
void calc1(const pol&a,int w,int L,int R){
if(R-L<=){
for(int i=L;i<=R;++i)as[qs[i].id]=cal(a,qs[i].x);
return;
}
int M=(L+R)/;
calc1(a%tr[w<<],w<<,L,M);
calc1(a%tr[w<<^],w<<^,M+,R);
}
int main(){
scanf("%d%d%d%d%d",&n,&b,&c,&d,&e);
for(int i=;i<n;++i)scanf("%d",a+i);
for(int i=;i<n;++i){
qs[i].x=(b*pw<P>(c,*i)+d*pw<P>(c,*i)+e)%P;
qs[i].id=i;
}
std::sort(qs,qs+n);
tr[]=pol(n);
mov(a,tr[]);
calc0(,,n-);
calc1(tr[],,,n-);
for(int i=;i<n;++i)printf("%d\n",(as[i]+P)%P);
return ;
}
bzoj 4929: 第三题的更多相关文章
- shanquan2的两年三题系列
好像只有2个月就退役啦 不管了,先说一下哪三题:多点求值.lcm.替罪羊树(bzoj3065) [upd0]2016.3.29 多点求值A掉啦,myy卡常数sxbk(不是说好的是shanquan2出的 ...
- Java-集合-第三题 有如下Student 对象, private String name; private int age; private int score; private String classNum; 其中,classNum 表示学生的班号,例如“class05”。 有如下List List list = new ArrayList(); l
第三题 有如下Student 对象, private String name; private int age; private int score; private String classNum; ...
- NOIP 2008提高组第三题题解by rLq
啊啊啊啊啊啊今天已经星期三了吗 那么,来一波题解吧 本题地址http://www.luogu.org/problem/show?pid=1006 传纸条 题目描述 小渊和小轩是好朋友也是同班同学,他们 ...
- NOIP2008提高组(前三题) -SilverN
此处为前三题,第四题将单独发布 火柴棒等式 题目描述 给你n根火柴棍,你可以拼出多少个形如“A+B=C”的等式?等式中的A.B.C是用火柴棍拼出的整数(若该数非零,则最高位不能是0).用火柴棍拼数字0 ...
- 《学习OpenCV》练习题第四章第三题b
#include <highgui.h> #include <cv.h> #include "opencv_libs.h" /* *<学习OpenCV ...
- 《学习OpenCV》练习题第四章第三题a
#include <highgui.h> #include <cv.h> #include "opencv_libs.h" #pragma comment ...
- 阿里聚安全攻防挑战赛第三题Android PwnMe解题思路
阿里聚安全攻防挑战赛第三题Android PwnMe解题思路 大家在聚安全挑战赛正式赛第三题中,遇到android app 远程控制的题目.我们今天带你一探究竟,如何攻破这道题目. 一.题目 购物应用 ...
- 剑指offier第三题
package 剑指office; /* * 第三题二维数组查找 * 在一个二维数组中,每一行都按照从左到右递增的顺序排序, * 每一列都按照从上到下递增的顺序排序.请完成一个函数,输入这样的一个二维 ...
- NOIP2005-普及组复赛-第三题-采药
题目描述 Description 辰辰是个天资聪颖的孩子,他的梦想是成为世界上最伟大的医师.为此,他想拜附近最有威望的医师为师.医师为了判断他的资质,给他出了一个难题.医师把他带到一个到处都是草药的山 ...
随机推荐
- windows下运用批处理实现一键自动开启多个应用
工作时,我每天早上到公司,打开自己的电脑,都会有几个固定的软件(myeclipse,飞信,firefox,foxmail等).文件夹和文件需要打开,每天如此,感到很烦,浪费时间做重复的工作,于是想到一 ...
- iOS UIlabel怎么加载html字符串 富文本的用法
要加载html字符串,用人说,直接用webView啊!但是,有时候我们只需要显示2行文字,如此少的内容却要在复杂的UI排版中加入一个占用资源较多的webview,得不偿失.这里要说的是,我们其实可以用 ...
- DevExpress v17.2新版亮点—WinForms篇(五)
用户界面套包DevExpress v17.2终于正式发布,本站将以连载的形式为大家介绍各版本新增内容.开篇介绍了DevExpress WinForms v17.2 Data Grid Control ...
- slf4j简单使用
一 slf4j+log4j 1.添加依赖 <!-- https://mvnrepository.com/artifact/org.slf4j/slf4j-log4j12 --> <d ...
- 移动端 解决自适应 和 多种dpr (device pixel ratio) 的 [淘宝] 解决方案 lib-flexible
其实H5适配的方案有很多种,网上有关于这方面的教程也非常的多. 不管哪种方法,都有其自己的优势和劣势. 为什么推荐使用Flexible库来做H5页面的终端设备适配呢? 原理 简单易懂 源码疑问 ...
- 结合P2P软件使用Ansible分发大文件
一 应用场景描述 现在我需要向50+数量的服务器分发Logstash新版本的rpm包,大概220MB左右,直接使用Ansible的copy命令进行传输,命令如下: 1 ansible all -m ...
- Redis在.net中的应用学习
在Redis的官网(http://redis.io/clients#c)上可以看到支持Redis C#的客户端. redis的网络连接方式和传统的rdbms相似,一种是长连接,一种是连接池,此处使用长 ...
- Redis学习第三课:Redis Hash类型及操作
Redis hash是一个string类型的field和value的映射表.它的添加.删除操作都是O(1)(平均).hash特别适用于存储对象.相较于对象的每个字段存在单个string类型.将一个对象 ...
- SWIFT显示底部的工具条
有以下页面显示我的讯息,用户可以点击右上角的编辑按钮进入删除状态.点击编辑按钮后,按钮文字改为“取消”,左上角的按钮变为“全选”,同时显示底部工具条带有“删除”按钮 实现起来挺简单的,在正常状态下点击 ...
- Ubuntu install TensorFlow
/******************************************************************************** * Ubuntu install T ...