题目描述

Alice 和 Bob 又在玩游戏。

对于一次游戏,首先 Alice 获得一个长度为 ​ 的序列 ​,Bob 获得一个长度为 ​ 的序列 bb。之后他们各从自己的序列里随机取出一个数,分别设为 ​,定义这次游戏的 ​次价值为​。

由于他们发现这个游戏实在是太无聊了,所以想让你帮忙计算对于 ​一次游戏​ 次价值的期望是多少。

由于答案可能很大,只需要求出模 ​下的结果即可。

输入输出格式

输入格式:

第一行两个整数 ​,分别表示 Alice 和 Bob 序列的长度。

接下来一行 ​ 个数,第 ​ 个数为 ​,表示 Alice 的序列。

接下来一行 ​ 个数,第 ​ 个数为 ​,表示 Bob 的序列。

接下来一行一个整数 ​,意义如上所述。

输出格式:

共 ​ 行,第 ​ 行表示一次游戏 ​ 次价值的期望。

输入输出样例

输入样例#1:

复制

1 1

1

2

3

输出样例#1:

复制

3

9

27

输入样例#2:

复制

2 8

764074134 743107904

663532060 183287581 749169979 7678045 393887277 27071620 13482818 125504606

6

输出样例#2:

复制

774481679

588343913

758339354

233707576

36464684

461784746

神仙题啊!!

前置知识:

生成函数。

Taylor展开。

NTT。

多项式求​。

​次价值得期望就是​。

用二项式定理将分子展开:

考虑将​展开:

我们设:

答案就是​

问题的关键就在于预处理出​。

预处理出这个多项式的方法是一个经典套路(根本不会)。























通过​的系数可以预处理出A和B。

解释一下,最后一步变换:

根据Taylor展开:

​求导如下:

复合函数求导遵循链式法则:

于是:

我们再​处展开,也就是令​

得到:

于是我们求出​然后就可以用第​项得系数来算出​和​。

代码(算​的逆元的时候乘爆了,调了一晚上):

include<bits/stdc++.h>

define ll long long

define N 300005

define mod 998244353

using namespace std;

inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch'-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return xf;}



ll n,m;

int t;

ll a[N],b[N];

ll fac[N<<2],inv[N<<2];



ll ksm(ll t,ll x) {

ll ans=1;

for(;x;x>>=1,t=t
t%mod)

if(x&1) ans=ans*t%mod;

return ans;

}

ll Mod(ll a) {

if(a<0) return a+mod;

if(a>=mod) return a-mod;

return a;

}

int rev[N<<2];

void NTT(ll *a,int d,int flag) {

static const ll G=3;

int n=1<<d;

for(int i=0;i<n;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<d-1);

for(int i=0;i<n;i++) if(i<rev[i]) swap(a[i],a[rev[i]]);

for(int s=1;s<=d;s++) {

int len=1<<s,mid=len>>1;

ll w=flag1?ksm(G,(mod-1)/len):ksm(G,mod-1-(mod-1)/len);

for(int i=0;i<n;i+=len) {

ll t=1;

for(int j=0;j<mid;j++,t=tw%mod) {

ll u=a[i+j],v=a[i+j+mid]
t%mod;

a[i+j]=Mod(u+v);

a[i+j+mid]=Mod(u-v+mod);

}

}

}

if(flag-1) {

ll inv=ksm(n,mod-2);

for(int i=0;i<n;i++) a[i]=a[i]*inv%mod;

}

}



struct poly {

ll g[N<<2],f[N<<2];

ll a[N<<2],ans[N<<2];

void solve(int l,int r) {

if(lr) return ;

int len=r-l+1,mid=l+r>>1;

solve(l,mid),solve(mid+1,r);

int d=ceil(log2(len+1));

for(int i=0;i<(1<<d);i++) g[i]=f[i]=0;

g[0]=f[0]=1;

for(int i=l;i<=mid;i++) f[i-l+1]=a[i];

for(int i=mid+1;i<=r;i++) g[i-mid]=a[i];

NTT(f,d,1),NTT(g,d,1);

for(int i=0;i<(1<<d);i++) f[i]=f[i]g[i]%mod;

NTT(f,d,-1);

for(int i=l;i<=r;i++) a[i]=f[i-l+1];

}

ll inv[N<<2],A[N<<2];

void Inverse(ll a,int d,ll f) {

f[0]=ksm(a[0],mod-2);

for(int s=1;s<=d;s++) {

int len=1<<s,Len=len<<1;

for(int i=0;i<len;i++) A[i]=a[i],A[i+len]=0;

NTT(A,s+1,1),NTT(f,s+1,1);

for(int i=0;i<Len;i++) f[i]=(2
f[i]-f[i]
f[i]%mod
A[i]%mod+mod)%mod;

NTT(f,s+1,-1);

for(int i=len;i<Len;i++) f[i]=0;

}

}

void DER(ll a,int n) {

for(int i=0;i<n;i++) a[i]=a[i+1]
(i+1)%mod;

a[n]=a[n-1]=0;

}

void INT(ll a,int n) {

for(int i=n-1;i>=0;i--) a[i+1]=a[i]
ksm(i+1,mod-2)%mod;

a[0]=0;

}

void Ln(ll *a,int d,ll f) {

Inverse(a,d,inv);

DER(a,1<<d);

NTT(inv,d+1,1),NTT(a,d+1,1);

for(int i=0;i<(1<<d+1);i++) f[i]=inv[i]
a[i]%mod;

NTT(f,d+1,-1);

INT(f,1<<d);

}

void mian(int n,int lim,ll *orig,ll f) {

for(int i=1;i<=n;i++) a[i]=orig[i];

solve(1,n);

int d=ceil(log2(lim));

a[0]=1;

Ln(a,d,ans);

for(int i=1;i<=(1<<d);i++) {

f[i]=ans[i];

f[i]=f[i]
::inv[i-1]%mod;

if(!(i&1)) f[i]=(mod-f[i])%mod;

}

}

}pre[2];

ll A[N<<2],B[N<<2];

int main() {

n=Get(),m=Get();

for(int i=1;i<=n;i++) a[i]=Get();

for(int i=1;i<=m;i++) b[i]=Get();

t=Get();

fac[0]=1;
for(int i=1;i<=t;i++) fac[i]=fac[i-1]*i%mod; inv[t]=ksm(fac[t],mod-2);
for(int i=t-1;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod; int lim=max(n,m)+t;
pre[0].mian(n,lim,a,A);
pre[1].mian(m,lim,b,B); int d=ceil(log2(t*2));
A[0]=n,B[0]=m;
for(int i=t+1;i<(1<<d);i++) A[i]=B[i]=0; NTT(A,d,1),NTT(B,d,1);



for(int i=0;i<(1<<d);i++) A[i]=A[i]*B[i]%mod;

NTT(A,d,-1);

ll invnm=ksm(1ll*n*m%mod,mod-2);

for(int i=1;i<=t;i++) {
cout<<fac[i]*A[i]%mod*invnm%mod<<"\n";
} return 0;

}

Luogu P4705 玩游戏的更多相关文章

  1. 洛谷 P4705 玩游戏 解题报告

    P4705 玩游戏 题意:给长为\(n\)的\(\{a_i\}\)和长为\(m\)的\(\{b_i\}\),设 \[ f(x)=\sum_{k\ge 0}\sum_{i=1}^n\sum_{j=1}^ ...

  2. Luogu 4705 玩游戏

    看见这个题依稀想起了$5$月月赛时候的事情,到现在仍然它感觉非常神仙. 游戏$k$次价值的期望答案 $$ans_k = \frac{1}{nm}\sum_{i = 1}^{n}\sum_{j = 1} ...

  3. P4705 玩游戏

    思路 超级麻烦... 写了一堆最后常数太大T飞了... 真的难受 发现solve函数可以不用把下一层复制上来,直接传指针就可以,下次再说写不写叭 思路 \[ ans_k=\sum_{i=1}^n\su ...

  4. 洛谷P4705 玩游戏 [生成函数,NTT]

    传送门 这是两个月之前写的题,但没写博客.现在回过头来看一下发现又不会了-- 还是要写博客加深记忆. 思路 显然期望可以算出总数再乘上\((nm)^{-1}\). 那么有 \[ \begin{alig ...

  5. [洛谷P4705]玩游戏

    题目大意:对于每个$k\in[1,t]$,求:$$\dfrac{\sum\limits_{i=1}^n\sum\limits_{j=1}^m(a_i+b_j)^k}{nm}$$$n,m,t\leqsl ...

  6. 洛谷 P4705 玩游戏

    题目分析 题目要求的是: \[ \sum_{i=1}^n\sum_{j=1}^m(a_i+b_j)^x(x\in [1,T]) \] 利用二项式定理化式子, \[ \begin{aligned} &a ...

  7. 洛谷P4705 玩游戏(生成函数+多项式运算)

    题面 传送门 题解 妈呀这辣鸡题目调了我整整三天--最后发现竟然是因为分治\(NTT\)之后的多项式长度不是\(2\)的幂导致把多项式的值存下来的时候发生了一些玄学错误--玄学到了我\(WA\)的点全 ...

  8. 【洛谷5月月赛】玩游戏(NTT,生成函数)

    [洛谷5月月赛]玩游戏(NTT,生成函数) 题面 Luogu 题解 看一下要求的是什么东西 \((a_x+b_y)^i\)的期望.期望显然是所有答案和的平均数. 所以求出所有的答案就在乘一个逆元就好了 ...

  9. [luogu]P1070 道路游戏[DP]

    [luogu]P1070 道路游戏 题目描述小新正在玩一个简单的电脑游戏.游戏中有一条环形马路,马路上有 n 个机器人工厂,两个相邻机器人工厂之间由一小段马路连接.小新以某个机器人工厂为起点,按顺时针 ...

随机推荐

  1. 第一个Quartz程序 (二)

    1 我们使用maven项目 2 创建一个job类,在execute()方法里写上业务逻辑代码. 3 在另外一个类中创建触发器,调度器,并且绑定job. 首先在项目的pom.xml引入需要的jar包. ...

  2. apache伪静态配置(URL重写)

    1.打开apache配置文件 httpd.conf .2.开启rewrite模块,去掉注释# #LoadModule rewrite_module modules/mod_rewrite 3.让apa ...

  3. 【手记】.net正则行尾匹配符$的问题

    本来想用正则Split一下sql语句中简单场景的的GO,于是用^GO$(配合忽略大小写和多行模式),可居然连这种情况都搞不掂: go 如果删掉$就能匹配了,但这显然不是办法,遂又在VS的C#交互窗口. ...

  4. C#编译错误 CS0009:未能打开元数据文件

    C#编译错误 CS0009:未能打开元数据文件 编译错误 说明: 在编译向该请求提供服务所需资源的过程中出现错误.请检查下列特定错误详细信息并适当地修改源代码. 编译器错误消息: CS0009: 未能 ...

  5. python基础小结

    1.常用的格式符号 2.input函数 3.在程序中,看到了%这样的操作符,这就是python中的格式化输出 age = 18 name = "xiaohua" print(&qu ...

  6. 洛谷P3600 随机数生成器(期望dp 组合数)

    题意 题目链接 Sol 一条重要的性质:如果某个区间覆盖了另一个区间,那么该区间是没有用的(不会对最大值做出贡献) 首先不难想到枚举最终的答案\(x\).这时我们需要计算的是最大值恰好为\(x\)的概 ...

  7. 【读书笔记】iOS-解析JSON

    JSON相比XML最显著的优点是不需要使用重量级的解析库,因为其本身就是面向数据的,而且非常容易转换成哈希字典.除此之外,JSON文档相比同样的XML文档更小.在网络宽带有限的情况下,你很容易在Iph ...

  8. 【读书笔记】iOS-nonatomic

    原子性atomicity(nonatomic),是关系线程安全的,但是会影响性能.如果确定不考虑线程安全问题可以使用nonatomic. 参考资料:<iPhone与iPad开发实战-iOS经典应 ...

  9. HTML中令人惊喜的全局属性

    1.accesskey 属性 : 规定激活元素的快捷键. 浏览器支持:几乎所有浏览器均 accesskey 属性,除了 Opera. 定义和用法 accesskey 属性规定激活(使元素获得焦点)元素 ...

  10. loadrunner 脚本开发-url解码

    url解码 by:授客 QQ:1033553122 脚本结构如下: Action.c中的代码如下: int htoi(char *s) { int value = 0; int c = 0; c = ...