题面

传送门

思路

首先,我们把这个输入的点的生成函数搞出来:

$C=\sum_{i=0}{lim}s_ixi$

其中$lim$为集合里面出现过的最大的数,$s_i$表示大小为$i$的数是否出现过

我们再设另外一个函数$F$,定义$F_k$表示总权值为$k$的二叉树个数

那么,一个二叉树显然可以通过两个子树(可以权值为0,也就是空子树)和一个节点构成

那么有如下求$F$的式子

$F_0=1$

$F_k=\sum_{i=0}^k s_i \sum_{j=0}^{k-i} F_j F_{k-i-j}$

显然这个式子是一个卷积的形式

我们把这个三个东西的卷积放到整个函数上来考虑(也就是把$F_kC_k$看做多项式的系数做卷积),可以得到:

$F=C\ast F\ast F +1$,这里的$\ast$是多项式乘法

由于多项式乘法和普通乘法一样满足各种定律之类的,所以我们可以解个一元二次方程

$C\ast F^2 - F + 1=0$

$F=\frac{1\pm \sqrt{1-4C}}{2C}$

先考虑上面是+的情况

此时${\lim_{x \to 0}}F(x)=+\infty$,舍去

再考虑-的情况

此时${\lim_{x \to 0}}F(x)=1$,OK

其实到这里就可以算了,一个多项式开方+多项式求逆就解决了

但是我们更进一步,再划一划这个式子,可以变成如下形式:

$F=\frac{2}{1+\sqrt{1-4C}}$

多项式开方+多项式求逆解决问题

Code:

这题超!级!卡!常!我这种人傻常数大的,只能参照爷稳稳的代码写一写了qwq

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cassert>
#include<cmath>
#include<ctime>
#define ll long long
#define MOD 998244353
using namespace std;
inline int read(){
int re=0,flag=1;char ch=getchar();
while(ch>'9'||ch<'0'){
if(ch=='-') flag=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9') re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
return re*flag;
}
int n,m,K;
ll A[600010],B[600010],C[600010],x[600010],y[600010];
int g=3,ginv,inv2,r[600010],lim,cnt;
ll qpow(ll a,ll b){
ll re=1;
while(b){
if(b&1) re=re*a%MOD;
a=a*a%MOD;b>>=1;
}
return re;
}
ll add(ll a,ll b){
a+=b;
return (a>=MOD)?a-MOD:a;
}
ll dec(ll a,ll b){
a-=b;
return (a<0)?a+MOD:a;
}
void ntt(ll *a,ll type){
int i,j,k;ll x,y,w,wn,inv,mid;
for(i=0;i<lim;i++) if(i<r[i]) swap(a[i],a[r[i]]);
for(mid=1;mid<lim;mid<<=1){
wn=qpow((type==1)?g:ginv,(MOD-1)/(mid<<1));
for(j=0;j<lim;j+=(mid<<1)){
w=1;
for(k=0;k<mid;k++,w=w*wn%MOD){
x=a[j+k];y=a[j+k+mid]*w%MOD;
a[j+k]=add(x,y);
a[j+k+mid]=dec(x,y);
}
}
}
if(~type) return;
inv=qpow(lim,MOD-2);
for(i=0;i<lim;i++) a[i]=a[i]*inv%MOD;
}
void getinv(ll *a,ll *b,int len,int lena){
if(len==1){b[0]=qpow(a[0],MOD-2);return;}
assert(len==((len>>1)<<1));
int i,mid=len>>1;
getinv(a,b,mid,lena);
lim=1;cnt=0;r[0]=0;
while(lim<=(len<<1)) lim<<=1,cnt++;
for(i=0;i<lim;i++) r[i]=((r[i>>1]>>1)|((i&1)<<(cnt-1))); for(i=0;i<len;i++) x[i]=a[i];
for(i=len;i<lim;i++) x[i]=0;
for(i=0;i<mid;i++) y[i]=b[i];
for(i=len;i<lim;i++) y[i]=0; ntt(x,1);ntt(y,1);
for(i=0;i<lim;i++) x[i]=y[i]*(2ll-y[i]*x[i]%MOD+MOD)%MOD;
ntt(x,-1);
for(i=0;i<len;i++) b[i]=x[i];
}
void getrt(ll *a,ll *b,int len){
if(len==1){b[0]=1;return;}
int i,mid=len>>1; getrt(a,b,mid);
memset(C,0,sizeof(C));getinv(b,C,len,len); lim=1;cnt=0;r[0]=0;
while(lim<=(len<<1)) lim<<=1,cnt++;
for(i=0;i<lim;i++) r[i]=((r[i>>1]>>1)|((i&1)<<(cnt-1))),x[i]=a[i];
for(i=len;i<lim;i++) x[i]=C[i]=b[i]=0; ntt(x,1);ntt(C,1);
for(i=0;i<lim;i++) x[i]=x[i]*C[i]%MOD;
ntt(x,-1); for(i=0;i<len;i++) b[i]=add(b[i],x[i])*inv2%MOD;
}
void init(){
ginv=qpow(g,MOD-2);
inv2=qpow(2,MOD-2);
}
int main(){
init();
n=read();m=read();int i,t1;A[0]=K=1;
for(i=1;i<=n;i++) t1=read(),A[t1]=dec(A[t1],4);
while(K<=m) K<<=1; getrt(A,B,K);
B[0]=add(B[0],1);
memset(C,0,sizeof(C));
getinv(B,C,K,K);
for(i=1;i<=m;i++) printf("%lld\n",(2*C[i]+MOD)%MOD);
}

[Codeforces438E][bzoj3625] 小朋友和二叉树 [多项式求逆+多项式开根]的更多相关文章

  1. FFT模板 生成函数 原根 多项式求逆 多项式开根

    FFT #include<iostream> #include<cstring> #include<cstdlib> #include<cstdio> ...

  2. 【BZOJ3625】【codeforces438E】小朋友和二叉树 生成函数+多项式求逆+多项式开根

    首先,我们构造一个函数$G(x)$,若存在$k∈C$,则$[x^k]G(x)=1$. 不妨设$F(x)$为最终答案的生成函数,则$[x^n]F(x)$即为权值为$n$的神犇二叉树个数. 不难推导出,$ ...

  3. 2019.01.01 bzoj3625:小朋友和二叉树(生成函数+多项式求逆+多项式开方)

    传送门 codeforces传送门codeforces传送门codeforces传送门 生成函数好题. 卡场差评至今未过 题意简述:nnn个点的二叉树,每个点的权值KaTeX parse error: ...

  4. bzoj 3625小朋友和二叉树 多项式求逆+多项式开根 好题

    题目大意 给定n种权值 给定m \(F_i表示权值和为i的二叉树个数\) 求\(F_1,F_2...F_m\) 分析 安利博客 \(F_d=F_L*F_R*C_{mid},L+mid+R=d\) \( ...

  5. NTT+多项式求逆+多项式开方(BZOJ3625)

    定义多项式$h(x)$的每一项系数$h_i$,为i在c[1]~c[n]中的出现次数. 定义多项式$f(x)$的每一项系数$f_i$,为权值为i的方案数. 通过简单的分析我们可以发现:$f(x)=\fr ...

  6. 【learning】多项式相关(求逆、开根、除法、取模)

    (首先要%miskcoo,这位dalao写的博客(这里)实在是太强啦qwq大部分多项式相关的知识都是从这位dalao博客里面学的,下面这篇东西是自己对其博客学习后的一些总结和想法,大部分是按照其博客里 ...

  7. 【BZOJ3456】轩辕朗的城市规划 无向连通图计数 CDQ分治 FFT 多项式求逆 多项式ln

    题解 分治FFT 设\(f_i\)为\(i\)个点组成的无向图个数,\(g_i\)为\(i\)个点组成的无向连通图个数 经过简单的推导(枚举\(1\)所在的连通块大小),有: \[ f_i=2^{\f ...

  8. bzoj 3456 城市规划——分治FFT / 多项式求逆 / 多项式求ln

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3456 分治FFT: 设 dp[ i ] 表示 i 个点时连通的方案数. 考虑算补集:连通的方 ...

  9. P6295-有标号 DAG 计数【多项式求逆,多项式ln】

    正题 题目链接:https://www.luogu.com.cn/problem/P6295 题目大意 求所有\(n\)个点的弱联通\(DAG\)数量. \(1\leq n\leq 10^5\) 解题 ...

随机推荐

  1. a链接打开另外的新页面

    在a标签添加target = "_blank" 属性即可

  2. IntelliJ IDEA 方法注释教程

    首先Ctrl +Alt +S ,打开Settings ,找到Live Template,然后点击右侧的绿色的“+”,选择Template Group 然后给新建的Group随便命个名 选中自己刚才创建 ...

  3. 设置vim tab为4个空格

    Vim 编辑器默认tab为8个空格,但对于pythoner来说,必须要调整到4个空格. 方法如下: 在~/.vimrc文件中加入下面设置: set ts=4 #设置tabstop为4个空格 set e ...

  4. PHP实现的一个时间帮助类

    背景 工作的过程中经常会遇到各种时间类的操作,因此封装了一个帮助工具类,提高代码的复用率 主要功能 根据相差的天数获取连续的时间段 /** * 根据相差的天数获取所有连续的时间段 * @param $ ...

  5. 两种js方法发起微信支付:WeixinJSBridge,wx.chooseWXPay区别

    原文链接:https://www.2cto.com/weixin/201507/412752.html 1.为什么会有两种JS方法可以发起微信支付? 当你登陆微信公众号之后,左边有两个菜单栏,一个是微 ...

  6. PHP审计(一)

    一.php中常见的危险函数和审计要点 危险函数(功能过于强大)    参数是否外部可控,有没有正确的过滤. PHP获取外界传入参数是通过下面几个全局函数的形式,所以审计参数传入经常要和下面几个变量打交 ...

  7. 用PHP关于Jquery表单插件ajaxForm里success不返回问题

    简单说一下吧,在用ajaxForm的时候,sucess突然之间不返回了,直接转到error里面去, 网页代码 ................. $('#add-type').ajaxForm({ d ...

  8. C语言实例解析精粹学习笔记——26

    实例26:阿拉伯数字转换为罗马数字,将一个整数n(1~9999)转换为罗马数字,其中数字和罗马数字的对应关系如下: 原书中的开发环境很老,我也没有花心思去研究.自己在codeblocks中进行开发的, ...

  9. 2015 acm taipei L-Reward the Troop(uvalive 7465)(找规律)

    原题链接 就大概说的是一个将军要给部下发勋章,他的部下以和别人不一样的勋章为荣,但是他没这么多钱,所以问你最少要多少钱 要求是每个人的上司是他的上两级,他的下两级是他的部下,每个人的勋章不能和他的上司 ...

  10. python-8错误调试测试

    1-错误处理 import logging try: print('try.......') r = 10/0 except ValueError as e: print('result:', e) ...