【JZOJ5740】【20190706】幻想世界
题目
小 $\omega $ 想要进行烟火表演,她一开始有\(n\)颗彗星和\(n\)颗陨石
如果小 \(\omega\) 有\(i\)颗彗星而没有陨石,那么她会消耗\(i\)颗彗星并得到\(a_i\)的火焰
如果小 \(\omega\) 有\(j\)颗陨石而没有彗星,那么她会消耗\(j\)颗陨石并得到\(a_j\)的火焰
否则小 \(\omega\) 有 \(p\) 的概率消耗一颗陨石,得到\(\alpha\)的火焰,有 \(q\) 的概率消耗一颗彗星,得到\(\beta\)的火焰
有\((1-p-q)\)的概率丢弃掉所有材料
设\(f_{i,j}\)表示小 \(\omega\) 有 \(i\) 颗彗星, \(j\) 颗陨石期望得到的火焰数量
求
\[\sum_{i=0}^{n}\sum_{j=0}^{n} h^{i(n+1)+j} \ f_{i,j}
\]
$ n \le 200000 $
题解
设\(\gamma = p \alpha + q \beta\),得到\(f_{i,j}\)的转移式
\[\begin{align}
f_{i,j} \ = \ pf_{i-1,j} + qf_{i,j-1} + \gamma
\end{align}
\]相当于向下贡献 $\times p $ ,向右贡献 \(\times q\) ,分别考虑\(a_i,b_i,\gamma\)的贡献和
\[\begin{align}
ans &= \\
&\gamma \sum_{i=1}^{n}\sum_{j=1}^{n} \sum_{x=i}^{n}\sum_{y=j}^{n} (^{x-i+y-j}_{x-i})p^{x-i}q^{y-j} h^{x(n+1)+y}\\
+&\sum_{i=1}^{n}a_i\sum_{x=i}^{n}\sum_{y=1}^{n}(^{x-i+y-1}_{x-i})p^{x-i}q^{y}h^{x(n+1)+y}\\
+&\sum_{j=1}^{n}b_j\sum_{x=1}^{n}\sum_{y=j}^{n}(^{x-1+y-j}_{y-j})p^{x}q^{y-j}h^{x(n+1)+y}\\
+&\sum_{i=1}^{n}a_ih^{i(n+1)} + b_{i}h^i \\
\end{align}
\]考虑前三项,交换求和顺序分别是:
\[\begin{align}
&\sum_{x=0}^{n}\sum_{y=0}^{n}(^{x+y}_x)p^xq^y(\sum_{i=x+1}^{n}h^{i(n+1)})(\sum_{j=y+1}^{n}h^{j}) \\
&\sum_{x=0}^{n}\sum_{y=1}^{n}(^{x+y-1}_{x}) p^xq^y (\sum_{i=x+1}^{n}a_{i-x}h^{i(n+1)})h^y \\
&\sum_{x=1}^{n}\sum_{y=0}^{n}(^{x-1+y}_{y}) p^xq^y h^{x(n+1)} (\sum_{j=y+1}^{n}b_{j-y}h^j) \\
\end{align}
\]拆开组合数可以直接看成关于\(x,y\)的卷积
第一个式子的系数直接预处理,第二个式子的系数可以预先做一遍卷积
话说我有些细节推错了好多遍........
Code
//JZOJ锅掉了,暂时只知道这份代码能跑过大样例
#include<bits/stdc++.h>
#define ll long long
#define mod 998244353
#define G 3
using namespace std;
const int N=1<<19;
int n,h,p,q,ta,tb,r,a[N],b[N],tx,ty,dx,dy,A[6][N],B[3][N],fac[N],inv[N],ny[N];
int hx[N],hy[N],sx[N],sy[N],sa[N],sb[N],rev[N],len,L,ans;
char gc(){
static char*p1,*p2,s[1000000];
if(p1==p2)p2=(p1=s)+fread(s,1,1000000,stdin);
return(p1==p2)?EOF:*p1++;
}
int rd(){
int x=0;char c=gc();
while(c<'0'||c>'9')c=gc();
while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+c-'0',c=gc();
return x;
}
int pw(int x,int y){
int re=1;if(y<0)y+=mod-1;
for(;y;y>>=1,x=1ll*x*x%mod)
if(y&1)re=1ll*re*x%mod;
return re;
}
void inc(int&x,int y){x+=y;if(x>=mod)x-=mod;}
int pls(int x,int y){int re=x+y;return re>=mod?re-mod:re;}
int mis(int x,int y){int re=x-y;return re<0?re+mod:re;}
void ntt(int*F,int sg){
for(int i=0;i<len;++i)if(i<rev[i])swap(F[i],F[rev[i]]);
for(int i=1;i<len;i<<=1){
int wn=pw(G,sg*(mod-1)/i/2);
for(int j=0;j<len;j+=i<<1){
int w=1;
for(int k=0;k<i;++k,w=(ll)w*wn%mod){
int x=F[j+k],y=(ll)w*F[j+k+i]%mod;
F[j+k]=pls(x,y);F[j+k+i]=mis(x,y);
}
}
}
if(!~sg){
int iv=pw(len,mod-2);
for(int i=0;i<len;++i)F[i]=(ll)F[i]*iv%mod;
}
}
int main(){
freopen("dream.in","r",stdin);
freopen("dream2.out","w",stdout);
n=rd();h=rd();ta=rd();tb=rd();
p=rd();p=(ll)p*pw(rd(),mod-2)%mod;
q=rd();q=(ll)q*pw(rd(),mod-2)%mod;
r=((ll)p*ta+(ll)q*tb)%mod;
for(int i=1;i<=n;++i)a[i]=rd();
for(int i=1;i<=n;++i)b[i]=rd();
ty=pw(h,n+1),tx=pw(ty,n+1);
dy=pw(h,mod-2),dx=pw(ty,mod-2);
for(int i=n;i;--i){
sx[i]=hx[i]=tx=(ll)tx*dx%mod;
sy[i]=hy[i]=ty=(ll)ty*dy%mod;
inc(sx[i],sx[i+1]);
inc(sy[i],sy[i+1]);
inc(ans,((ll)hx[i]*a[i]+(ll)hy[i]*b[i])%mod);
}
ny[1]=1;for(int i=2;i<=n<<1;++i)ny[i]=(ll)(mod-mod/i)*ny[mod%i]%mod;
for(int i=fac[0]=inv[0]=1;i<=n<<1;++i){
inv[i]=1ll*inv[i-1]*ny[i]%mod;
fac[i]=1ll*fac[i-1]*i%mod;
}
for(L=0,len=1;len<=n<<1;len<<=1,L++);
for(int i=0;i<len;++i)rev[i]=(rev[i>>1]>>1)|((i&1)<<(L-1));
for(int i=1,tp=1,tq=1;i<=n;++i){
A[2][i-1]=(ll)tp*sx[i]%mod*inv[i-1]%mod;
A[3][i-1]=(ll)tq*sy[i]%mod*inv[i-1]%mod;
tp=(ll)tp*p%mod;tq=(ll)tq*q%mod;
A[0][i-1]=(ll)tp*hx[i]%mod*inv[i-1]%mod;
A[1][i-1]=(ll)tq*hy[i]%mod*inv[i-1]%mod;
}
reverse(a,a+n+1);reverse(b,b+n+1);
ntt(a,1);ntt(b,1);ntt(hx,1);ntt(hy,1);
for(int i=0;i<len;++i){
a[i]=(ll)a[i]*hx[i]%mod;
b[i]=(ll)b[i]*hy[i]%mod;
}
ntt(a,-1);ntt(b,-1);//ntt(hx,-1);ntt(hy,-1);
for(int i=0,tp=1,tq=1;i<=n;++i){
A[4][i]=(ll)tp*a[n+i]%mod*inv[i]%mod;
A[5][i]=(ll)tq*b[n+i]%mod*inv[i]%mod;
tp=(ll)tp*p%mod;tq=(ll)tq*q%mod;
}
for(int i=0;i<6;++i)ntt(A[i],1);
for(int i=0;i<len;++i){
// B[0][i]=(ll)A[2][i]*A[3][i]%mod;
// B[1][i]=(ll)A[4][i]*A[1][i]%mod;
// B[2][i]=(ll)A[0][i]*A[5][i]%mod;
B[0][i]=((ll)A[2][i]*A[3][i]%mod*r+(ll)A[4][i]*A[1][i]+(ll)A[0][i]*A[5][i])%mod;
}
//for(int i=0;i<3;++i)ntt(B[i],-1);
ntt(B[0],-1);
for(int i=0;i<=n<<1;++i){
//inc(ans , ((ll)r*B[0][i] + B[1][i] + B[2][i])%mod * fac[i] %mod);
inc(ans , (ll)B[0][i]*fac[i]%mod);
}
cout<<ans<<endl;
return 0;
}
【JZOJ5740】【20190706】幻想世界的更多相关文章
- 【JZOJ5739】【20190706】毒奶
题目 有\(n\)个现实城市,另有\(n\)个幻想城市 原图中在现实城市存在\(m\)条边,在幻想城市存在\(m-1-n\)条边 一个排列是合法的当且进当显示城市 \(i\) 向幻想城市 \(p_i\ ...
- 【JZOJ5738】【20190706】锁屏杀
题目 $n \le 2000 $ 题解 \(B\)的数字可以对1440取模,对三个图分别进行\(dp\)即可 时间复杂度\(O(n\times 1440 \times 10)\) Code #incl ...
- 2019-07-06 sql 连续出现次数
由手机通讯记录界面想到的问题 SELECT CASE WHEN AA.num=1 THEN AA.Tel ELSE AA.Tel+'('+CASt(AA.num AS VARCHAR(4))+')' ...
- 2019-07-06 sql备忘 连续取最大
连续最大: SELECT M.* FROM #temp MINNER JOIN (SELECT ISNULL(A.score,0)-b.score AS score,B.id FROM #temp A ...
- QQ 腾讯QQ(简称“QQ”)是腾讯公司开发的一款基于Internet的即时通信(IM)软件
QQ 编辑 腾讯QQ(简称“QQ”)是腾讯公司开发的一款基于Internet的即时通信(IM)软件.腾讯QQ支持在线聊天.视频通话.点对点断点续传文件.共享文件.网络硬盘.自定义面板.QQ邮箱等多种功 ...
- js+CSS实现模拟华丽的select控件下拉菜单效果
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- HTML5 CSS3专题 诱人的实例 CSS3打造百度贴吧的3D翻牌效果
首先感谢w3cfuns的老师~ 今天给大家带来一个CSS3制作的翻牌效果,就是鼠标移到元素上,感觉可以看到元素背后的信息.大家如果制作考验记忆力的连连看.扑克类的游戏神马的,甚至给女朋友写一些话语,放 ...
- DOIS 2019 DevOps国际峰会北京站来袭~
DevOps 国际峰会是国内唯一的国际性 DevOps 技术峰会,由 OSCAR 联盟指导.DevOps 时代社区与高效运维社区联合主办,共邀全球80余名顶级专家畅谈 DevOps 体系与方法.过程与 ...
- 模拟select控件,css模拟下拉
<!DOCTYPE html > <head> <meta http-equiv="Content-Type" content="t ...
随机推荐
- thinkphp的普通查询与表达式查询
一.普通查询方式 a.字符串:$arr=$m->where("sex=0 and username='gege'")->find();//字符串需要加引号 b.数组 $ ...
- 四种软件开发模式:tdd、bdd、atdd和ddd的概念
看一些文章会看到TDD开发模式,搜索后发现有主流四种软件开发模式,这里对它们的概念做下笔记. TDD:测试驱动开发(Test-Driven Development) 测试驱动开发是敏捷开发中的一项核心 ...
- 一次kuberneets evicted的历险
一.概述 kubernetes 的eviction检测diskpresure,检测的是kubelet的root-dir.kubelet的默认root-dir是/var/lib/kubelet,可以使用 ...
- MySQL如何定位慢sql
MySQL如何定位慢sql MySQL"慢SQL"定位 数据库调优我个人觉得必须要明白两件事 1.定位问题(你得知道问题出在哪里,要不然从哪里调优呢) 2.解决问题(这个没有基本的 ...
- iviewUI 前端静态页面实现增删改查分页
完整代码部分 (仅供参考哈): <template> <div> <label prop="name"> 姓名: </label> ...
- Python: 把txt文件转换成csv
最近在项目上需要批量把txt文件转成成csv文件格式,以前是手动打开excel文件,然后导入txt来生产csv文件,由于这已经变成每周需要做的事情,决定用python自动化脚本来实现,思路: 读取文件 ...
- Service Mesh服务网格新生代--Istio
原文: 数人云|万字解读:Service Mesh服务网格新生代--Istio 参考: istio 简介 Istio是啥?一文带你彻底了解! 使用Istio治理微服务入门 Istio 流量管理 ist ...
- 基于vue+springboot+docker网站搭建【三】安装docker
安装docker 我们整个工程都会基于docker来跑,因此首先安装docker.在安装所有软件之前,我们需要先下载一个下载工具. 1.下载工具安装 yum-utils yum install -y ...
- Ivanti的垃圾软件landesk
landesk是Ivanti公司推出的终端管理工具,这个工具垃圾就垃圾在无法卸载,进程杀不死.文件删不掉,奉劝大家千万不要安装这个软件.前些天公司的IT部门一直在催促员工安装这个软件,我一时糊涂安装了 ...
- xss学习
1.了解xss的定义 2.理解xss的原理:反射型和存储型 3.理解xss的攻击方式 4.掌握xss的防御措施