「NOI十联测」黑暗

\(n\) 个点的无向图,每条边都可能存在,一个图的权值是连通块个数的 \(m\) 次方,求所有可能的图的权值和。(n≤30000,m≤15)

令\(ans[n][m]\)为n个点,图的权值是连通块个数的 m 次方,最终的答案.

\(g[i]\)为\(i\)个点无向联通图数目.

一个图,若连通块数T增加1,那么贡献为\(\displaystyle (T+1)^m=\sum_{i=0}^m C_m^i T^i\).

这里选取\(i\)个点组成一个连通块(包含一个指定的点),联通块数目增加1.

\[ans[n][m]=\sum_{i=1}^{n} g_i C_{n-1}^{i-1} (\sum_{j=0}^m C_m^j ans[n-i][j])
\]
\[ans[n][m]=m!(n-1)! \sum_{i=1}^{n} {g_i \over (i-1)!}{{\sum_{j=0}^m {ans[n-i][j]\over j! (m-j)!}} \over (n-i)!}
\]

令 \(\displaystyle T[i][j]=\sum_{k=0}^j {ans[i][k]\over k! (j-k)!}\)

\[ans[n][m]=m!(n-1)! \sum_{i=1}^{n} {g_i \over (i-1)!}{ T[n-i][m] \over (n-i)!}
\]

枚举m.

m=0的时候ans即为图的个数.

之后分治FFT...

复杂度\(O(nmlog(n^2)+nm^2)\).

代码里把\(i=n\)单独提出,其实可以不提出。

#include<bits/stdc++.h>
#define rep(q,a,b) for(int q=a,q##_end_=b;q<=q##_end_;++q)
#define dep(q,a,b) for(int q=a,q##_end_=b;q>=q##_end_;--q)
#define mem(a,b) memset(a,b,sizeof a )
#define debug(a) cerr<<#a<<' '<<a<<"___"<<endl
using namespace std;
void in(int &r) {
static char c;
r=0;
while(c=getchar(),!isdigit(c));
do r=(r<<1)+(r<<3)+(c^48);
while(c=getchar(),isdigit(c));
}
const int mn=30005;
const int mod=998244353;
int mul[mn],inv[mn];
int all[mn],un_able[mn],_able[mn];
int mlv(int x,int v){
int ans=1;
while(v){
if(v&1)ans=1LL*ans*x%mod;
x=1LL*x*x%mod,v>>=1;
}
return ans;
}
namespace NTT{
const int g=3;
const int gg=332748118;
int n,lim,to[66000];
void DFT(int* a,int iv){
rep(q,0,n-1)if(to[q]<q)swap(a[q],a[to[q]]);
for(int len=2;len<=n;len<<=1){
int mv=1,ml=mlv(iv?g:gg,(mod-1)/len),sp=len>>1;
for(int k=0;k<sp;++k){
for(int* p=a;p!=a+n;p+=len){
int mid=1LL*p[k+sp]*mv%mod;
p[k+sp]=(p[k]-mid)%mod;
p[k]=(p[k]+mid)%mod;
}
mv=1LL*mv*ml%mod;
}
}
}
void solve(int* a,int len1,int* b,int len2,int* c){
n=1,lim=0;
while(n<(len1+len2))n<<=1,++lim;
rep(q,len1,n-1)a[q]=0;
rep(q,len2,n-1)b[q]=0;
rep(q,0,n-1)to[q]=to[q>>1]>>1 | (q&1)<<(lim-1);
DFT(a,1),DFT(b,1);
rep(q,0,n-1)c[q]=1LL*a[q]*b[q]%mod;
DFT(c,0);
int inv_n=mlv(n,mod-2);
rep(q,0,n-1)c[q]=1LL*c[q]*inv_n%mod;
rep(q,0,n-1)c[q]=(c[q]+mod)%mod;
}
}
int dp[mn],dp1[mn],A[66000],B[66000],C[66000];
void solve(int l,int r){
if(l>r)return;
if(l==r){
un_able[l]=1LL*un_able[l]*mul[l-1]%mod;
_able[l]=((long long)all[l]-un_able[l]+mod)%mod;
return;
}
int mid=l+r>>1;
solve(l,mid);
rep(q,l,mid)A[q-l]=1LL*_able[q]*inv[q-1]%mod;
rep(q,1,r-l)B[q-1]=1LL*all[q]*inv[q]%mod;
NTT::solve(A,mid-l+1,B,r-l,C);
rep(q,mid+1,r)un_able[q]=(un_able[q]+C[q-l-1])%mod;
solve(mid+1,r);
}
int ans[16][30005],mid[16][30005],nw;
int *md,*md1;
void cdq(int l,int r){
if(l>r)return;
if(l==r){
md[l]=1LL*mul[nw]*mul[l-1]%mod*md[l]%mod;
md[l]=(md[l]+_able[l])%mod;
rep(q,0,nw)md1[l]=(md1[l]+1LL*ans[q][l]*inv[q]%mod*inv[nw-q])%mod;
return;
}
int mid=l+r>>1;
cdq(l,mid);
rep(q,l,mid)A[q-l]=1LL*md1[q]*inv[q]%mod;
rep(q,1,r-l)B[q-1]=1LL*_able[q]*inv[q-1]%mod;
NTT::solve(A,mid-l+1,B,r-l,C);
rep(q,mid+1,r)md[q]=(md[q]+C[q-l-1])%mod;
cdq(mid+1,r);
}
struct qr{
int a,b;
}st[1005];
int main(){
freopen("dark.in","r",stdin);
freopen("dark.out","w",stdout); int T,m,n=15;
in(T);
rep(q,1,T)in(st[q].a),in(st[q].b),n=max(n,st[q].a); mul[0]=1;
rep(q,1,n)mul[q]=1LL*mul[q-1]*q%mod;
inv[0]=inv[1]=1;
rep(q,2,n)inv[q]=1LL*(mod-mod/q)*inv[mod%q]%mod;
rep(q,1,n)inv[q]=1LL*inv[q]*inv[q-1]%mod; all[0]=1;
rep(q,1,n)all[q]=mlv(2,q*(q-1)>>1); solve(1,n); md=ans[0],md1=mid[0];
rep(q,1,n)md[q]=md1[q]=mlv(2,q*(q-1)>>1);
rep(q,1,15){
md=ans[q],md1=mid[q],nw=q;
cdq(1,n);
} rep(q,1,T)printf("%d\n",ans[st[q].b][st[q].a]);
return 0;
}

「NOI十联测」黑暗的更多相关文章

  1. 「NOI十联测」深邃

    「NOI十联测」深邃 要使得最大的连通块最小,显然先二分答案. 先固定1结点为根. 对于一个果实,显然是先处理子树中未分配的点,再向外延伸. 每个结点记录一个\(si[]\),表示子树中未分配的点数, ...

  2. 「NOI十联测」奥义商店

    「NOI十联测」奥义商店 若lzz想花费最少的钱,那么显然要选择数目较少的颜色. 先考虑暴力的写法. 每次向两边统计,每个物品要求被买的概率可以由上一个物品推出. now=1;//now 被买概率 M ...

  3. 「NOI十联测」反函数

    30pts 令(为1,)为-1: 暴力枚举每个点为起始点的路径,一条路径是合法的当且仅当路径权值和为0且路径上没有出现过负数. 将所有答案算出. 100pts 使用点分治. 要求知道经过重心root的 ...

  4. HHHOJ #153. 「NOI模拟 #2」Kotomi

    抽代的成分远远大于OI的成分 首先把一个点定为原点,然后我们发现如果我们不旋转此时答案就是所有位置的\(\gcd\) 如果要选择怎么办,我们考虑把我们选定的网格边连同方向和大小看做单位向量\(\vec ...

  5. HHHOJ #151. 「NOI模拟 #2」Nagisa

    计算几何板子题(我才没有拷板子的说--) 众所周知,三角形的重心坐标是\((\frac{x_1+x_2+x_3}{3},\frac{y_1+y_2+y_3}{3})\) 然后我们发现如果我们有一个点集 ...

  6. Solution -「SHOI2016」「洛谷 P4336」黑暗前的幻想乡

    \(\mathcal{Description}\)   link.   有一个 \(n\) 个结点的无向图,给定 \(n-1\) 组边集,求从每组边集选出恰一条边最终构成树的方案树.对 \(10^9+ ...

  7. Solution -「NOI 模拟赛」彩色挂饰

    \(\mathcal{Description}\)   给定一个含 \(n\) 个点 \(m\) 条边的简单无向图,设图中最大点双的大小为 \(s\),则保证 \(s\le6\).你将要用 \(k\) ...

  8. Solution -「NOI 模拟赛」出题人

    \(\mathcal{Description}\)   给定 \(\{a_n\}\),求一个 \(\{b_{n-1}\}\),使得 \(\forall x\in\{a_n\},\exists i,j\ ...

  9. 「洛谷 P1801」黑匣子

    好像很久没有更过博客了,因为博主这几周很忙.其实是在搞颓. 题意很难懂,所以就不重复了.其实是懒. 一眼看上去这是个 \(Splay\) 裸题,直接插入一个数,查询区间第 \(K\) 大,但是这样太不 ...

随机推荐

  1. 对XSS的插入的新了解,灵感来自天驿安全

    此次针对的是通过Get请求进行插入的XSS语句,或者dom型的xss,也算是了解到的新的插入方式 首先,JavaScript语言中存在拼接性 可以通过代审后闭合前置语句进行self测试是否可以拼接 s ...

  2. 主流的 API 架构

    1. RPC:调用另一个系统的函数 2. SOAP:使数据作为服务可用 3. REST:使数据作为资源可用 4. GraphQL:仅请求所需要的数据 1. RPC:调用另一个系统的函数 远程过程调用是 ...

  3. Docker如何制作镜像-Dockerfile的使用

    1:什么是Dockerfile Dockerfile是一个文本文档,可以通过docker build 命令构建成一个镜像. 我们可以在Dockerfile中定义一系列的命令,构建出我们想要的镜像. 想 ...

  4. Zookeeper基础教程(三):Zookeeper连接使用—zkCli

    上一篇介绍Zookeeper的安装,并介绍了使用ZooInspector连接Zookeeper,这里主要介绍以命令行的形式介绍Zookeeper 假如我们已经安装了Zookeeper集群,集群中的安装 ...

  5. Oracle打怪升级之路二【视图、序列、游标、索引、存储过程、触发器】

    前言 在之前 <Oracle打怪升级之路一>中我们主要介绍了Oracle的基础和Oracle常用查询及函数,这篇文章作为补充,主要介绍Oracle的对象,视图.序列.同义词.索引等,以及P ...

  6. html基础 有语义的布局标签

  7. 双buffer实现无锁切换

    大家好,我是雨乐! 在我们的工作中,多线程编程是一件太稀松平常的事.在多线程环境下操作一个变量或者一块缓存,如果不对其操作加以限制,轻则变量值或者缓存内容不符合预期,重则会产生异常,导致进程崩溃.为了 ...

  8. Unity3D开发入门教程(二)—— Lua入门

    五邑隐侠,本名关健昌,12年游戏生涯. 本教程以 Unity 3D + VS Code + C# + tolua 为例. 如果你还没有编程基础,建议你先学习一些编程基础.本文不是完全菜鸟教程,主要针对 ...

  9. 初识python: 回调函数

    回调函数 简单理解就是:将一个函数通过参数的形式传递给另一个函数 #!/user/bin env python # author:Simple-Sir # time:2019/8/9 10:49 # ...

  10. 初识python: 自定义函数

    什么是函数? 函数是组织好的,可重复使用的,用来实现单一,或相关联功能的代码段.函数能提高应用的模块性,和代码的重复利用率. 函数的定义方法: def test(x): '函数定义方法' x+=1 r ...