HDU 5552 Bus Routes(NTT+分治)
题意
给定 \(n\) 个点,任意连边,每条边有 \(m\) 种颜色可选,求带环连通图的方案数。
\(1\leq n\leq 10000\)
\(1\leq m < 2^{31}\)
思路
直接求带环连通图显然比较难求,正难则反,考虑容斥。用连通图的个数减去无环连通图(树)的个数。
\(n\) 个节点的无根树,每个节点有区别,可以直接套用公式 \(n^{n-2}\) 。而再考虑边的颜色,就是 \(m^{n-1}n^{n-2}\) 。
我们设 \(n\) 个点,考虑边的颜色,构成不同连通图的方案数为 \(f(n)\) 。
直接求连通图还是不方便,那么我们再容斥:用图的个数减不连通图的个数,\(n\) 个点,考虑边的颜色,可以有 \((m+1)^{n(n+1)\over2}\) 种情况,设之为 \(g(n)\)。
有一个小 \(\text{trick}\) ,我们固定一个点,选一些点和它构成一个连通块,剩下的点任意构图,显然这样是可以不重不漏的,转移式如下
\]
化简得
\]
这样就是一个 \(n^2\) 的 \(dp\) 式,并且形式上满足多项式乘法的形式,只是 \(f\) 在右边出现了。
那我们只能考虑左边对右边的转移,不难想到\(\text{CDQ}\)分治。
void CDQ(int l,int r)
{
if(l==r){/*转移常量给dp[l]*/return;}
int mid=(l+r)>>1;
CDQ(l,mid);
/*处理[l,mid]的多项式和转移给[mid+1,r]的多项式*/
_Polynomial::multiply(/**/);
/*转移结果给dp[mid+1,r]*/
CDQ(mid+1,r);
return;
}
代码流程如上,在分治过程中考虑左边转移给右边,需保证在转移前,左边的值以计算完毕。
\(\text{dp}\)式一般写成 \(dp_i=A_i\cdot \sum dp_jf_{i-j}+B_i\) 看的会比较清晰。
代码
#include<bits/stdc++.h>
#define FOR(i,x,y) for(int i=(x),i##END=(y);i<=i##END;++i)
#define DOR(i,x,y) for(int i=(x),i##END=(y);i>=i##END;--i)
using namespace std;
template<typename T,typename _T>inline bool chk_min(T &x,const _T y){return y<x?x=y,1:0;}
template<typename T,typename _T>inline bool chk_max(T &x,const _T y){return x<y?x=y,1:0;}
typedef long long ll;
const int P=152076289;
const int N=1<<14|5;
namespace _Maths
{
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
void exgcd(ll a,ll b,ll &x,ll &y)
{
if(!b){x=1,y=0;return;}
exgcd(b,a%b,y,x),y-=a/b*x;
}
ll Pow(ll a,ll p,ll P)
{
ll res=1;
for(;p>0;p>>=1,(a*=a)%=P)if(p&1)(res*=a)%=P;
return res;
}
ll inv(ll a,ll P){ll x,y;exgcd(a,P,x,y);return (x%P+P)%P;}
};
using namespace _Maths;
namespace _Polynomial
{
const int g=106;
int A[N<<1],B[N<<1];
int w[N<<1],r[N<<1];
void DFT(int *a,int op,int n)
{
FOR(i,0,n-1)if(i<r[i])swap(a[i],a[r[i]]);
for(int i=2;i<=n;i<<=1)
for(int j=0;j<n;j+=i)
for(int k=0;k<i/2;k++)
{
int u=a[j+k],t=(ll)w[op==1?n/i*k:(n-n/i*k)&(n-1)]*a[j+k+i/2]%P;
a[j+k]=(u+t)%P,a[j+k+i/2]=(u-t)%P;
}
if(op==-1)
{
int I=inv(n,P);
FOR(i,0,n-1)a[i]=(ll)a[i]*I%P;
}
}
void multiply(const int *a,const int *b,int *c,int n1,int n2)
{
int n=1;
while(n<n1+n2-1)n<<=1;
FOR(i,0,n1-1)A[i]=a[i];
FOR(i,0,n2-1)B[i]=b[i];
FOR(i,n1,n-1)A[i]=0;
FOR(i,n2,n-1)B[i]=0;
FOR(i,0,n-1)r[i]=(r[i>>1]>>1)|((i&1)*(n>>1));
w[0]=1,w[1]=Pow(g,(P-1)/n,P);
FOR(i,2,n-1)w[i]=(ll)w[i-1]*w[1]%P;
DFT(A,1,n),DFT(B,1,n);
FOR(i,0,n-1)A[i]=(ll)A[i]*B[i]%P;
DFT(A,-1,n);
FOR(i,0,n1+n2-2)c[i]=(A[i]+P)%P;
}
};
int A[N],B[N],C[N<<1];
int fac[N],ifac[N],f[N],g[N];
int n;ll m;
void CDQ(int l,int r)
{
if(l==r){f[l]=(g[l]-(ll)fac[l-1]*f[l]%P)%P;return;}
int mid=(l+r)>>1;
CDQ(l,mid);
FOR(i,l,mid)A[(i)-l]=(ll)f[i]*ifac[i-1]%P;
FOR(i,1,r-l)B[(i)-1]=(ll)g[i]*ifac[i]%P;
_Polynomial::multiply(A,B,C,mid-l+1,r-l);
FOR(i,mid+1,r)f[i]=((ll)f[i]+C[(i)-l-1])%P;
CDQ(mid+1,r);
}
int main()
{
fac[0]=fac[1]=1;FOR(i,2,N-1)fac[i]=(ll)fac[i-1]*i%P;
ifac[0]=ifac[1]=1;FOR(i,2,N-1)ifac[i]=(ll)(P-P/i)*ifac[P%i]%P;
FOR(i,2,N-1)ifac[i]=(ll)ifac[i-1]*ifac[i]%P;
int T;
scanf("%d",&T);
FOR(Ti,1,T)
{
scanf("%d%lld",&n,&m);
FOR(i,1,n)f[i]=0;
FOR(i,1,n)g[i]=Pow(m+1,(ll)i*(i-1)/2,P);
CDQ(1,n);
printf("Case #%d: %lld\n",Ti,(((ll)f[n]-Pow(n,n-2,P)*Pow(m,n-1,P))%P+P)%P);
}
return 0;
}
HDU 5552 Bus Routes(NTT+分治)的更多相关文章
- hdu 5552 Bus Routes
hdu 5552 Bus Routes 考虑有环的图不方便,可以考虑无环连通图的数量,然后用连通图的数量减去就好了. 无环连通图的个数就是树的个数,又 prufer 序我们知道是 $ n^{n-2} ...
- HDU 5552 Bus Routes(2015合肥现场赛A,计数,分治NTT)
题意 给定n个点,任意两点之间可以不连边也可以连边.如果连边的话可以染上m种颜色. 求最后形成的图,是一个带环连通图的方案数. 首先答案是n个点的图减去n个点能形成的树. n个点能形成的树的方案数比 ...
- HDU 5322 Hope ——NTT 分治 递推
发现可以推出递推式.(并不会) 然后化简一下,稍有常识的人都能看出这是一个NTT+分治的情况. 然而还有更巧妙的方法,直接化简一下递推就可以了. 太过巧妙,此处不表,建议大家找到那篇博客. 自行抄写 ...
- URAL 1137 Bus Routes(欧拉回路路径)
1137. Bus Routes Time limit: 1.0 secondMemory limit: 64 MB Several bus routes were in the city of Fi ...
- hdu 3842 Machine Works(cdq分治维护凸壳)
题目链接:hdu 3842 Machine Works 详细题解: HDU 3842 Machine Works cdq分治 斜率优化 细节比较多,好好体会一下. 在维护斜率的时候要考虑x1与x2是否 ...
- [LeetCode] Bus Routes 公交线路
We have a list of bus routes. Each routes[i] is a bus route that the i-th bus repeats forever. For e ...
- [Swift]LeetCode815. 公交路线 | Bus Routes
We have a list of bus routes. Each routes[i]is a bus route that the i-th bus repeats forever. For ex ...
- LeetCode解题报告—— Bus Routes
We have a list of bus routes. Each routes[i] is a bus route that the i-th bus repeats forever. For e ...
- [LeetCode] 815. Bus Routes 公交路线
We have a list of bus routes. Each routes[i] is a bus route that the i-th bus repeats forever. For e ...
随机推荐
- 11.vue 数据交互
vue new Vue({ el,选择器 string/obj 不能选择html/body data, methods, template string/obj //生命周期 -- 虚拟DOM 1.初 ...
- 8、路由 router
路由:router 用户功能 /user ----> index.html /user/login ----> login.html /user/reg ----> reg.html ...
- curl HTTP Header
对于"User-Agent", "Cookie", "Host"这类标准的HTTP头部字段,通常会有另外一种设置方法.curl命令提供了 ...
- 简述TCP三次握手和四次挥手过程
TCP握手协议 在TCP/IP协议中,TCP协议提供可靠的连接服务,采用三次握手建立一个连接.第一次握手:建立连接时,客户端发送syn包(syn=j)到服务器,并进入SYN_SEND状态,等待服务器确 ...
- linux基本介绍
Linux介绍 操作系统: 主要作用是管理好硬件设备,并为用户和应用程序提供简单的接口,以便于使用.作为中间人链接软件和硬件. 不同领域的操作系统: 1.桌面操作系统 Windows(用户群大).ma ...
- mysql提权常用方法。 hack某某
一般是root权限,知道mysql root权限,root账号密码 启动项提权:原理:利用高权限的root写入一个vbs脚本到启动项,再通过一些方法如ddos,社工管理员之类的方法来让服务器重启,运行 ...
- Laravel开发采坑系列问题
2017年12月22日17:40:03 不定时更新 版本5.4.X 一下是可能会遇到的坑 1,必须的写路由转发才能访问控制器,当然你可以自动路由访问,但是需要些匹配规则,其实还是转发了 好多人讨论过自 ...
- java Int 和 String 之间的转换
String 转换成 int Integer.parseInt(formParams.get("id")) int 转换成 string Integer.toString(id);
- 转:ArcGIS API For JavaScript官方文档(二十)之图形和要素图层——①Graphics概述
原文地址:ArcGIS API For JavaScript官方文档(二十)之图形和要素图层——①Graphics概述 ArcGIS JavaScript API允许在地图上绘制graphic(图形) ...
- 目标检测(六)YOLOv2__YOLO9000: Better, Faster, Stronger
项目链接 Abstract 在该论文中,作者首先介绍了对YOLOv1检测系统的各种改进措施.改进后得到的模型被称为YOLOv2,它使用了一种新颖的多尺度训练方法,使得模型可以在不同尺寸的输入上运行,并 ...