[多项式] COGS 有标号的二分图计数系列

最新推荐文章于 2024-09-29 17:24:23 发布
最新推荐文章于 2024-09-29 17:24:23 发布
阅读量691 收藏
点赞数
分类专栏： 多项式
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/70664725
版权

我是抄的 orzz

QAQ_bipartite_one

∑k=0nCkn∗2(n−k)∗k

注意这里是 2(n−k)∗k 不是 2(n−k)+k 怎么卷？
2(n−k)∗k=(2√)n2(2√)k2∗(2√)(n−k)2 

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef long long ll;

const int P=998244353;
const int G=3;
const int SQRT=116195171;
const int INV=(P+1)>>1;
const int N=300005;

inline ll Pow(ll a,int b){
  ll ret=1; for (;b;b>>=1,a=a*a%P) if (b&1) ret=ret*a%P; return ret;
}
inline ll Inv(ll a){
  return Pow(a,P-2);
}
ll fac[N];
inline void Init(int n){
  fac[0]=1; for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%P;
}


int num;
int w[2][N];
inline void Pre(int m){
  num=m;
  int g=Pow(G,(P-1)/m);
  w[0][0]=w[1][0]=1;
  for (int i=1;i<num;i++) w[1][i]=(ll)w[1][i-1]*g%P;
  for (int i=1;i<num;i++) w[0][i]=w[1][num-i];
}

int R[N];

inline void FFT(int n,int *a,int r){
  for (int i=0;i<n;i++) if (i<R[i]) swap(a[i],a[R[i]]);
  for (int i=1;i<n;i<<=1)
    for (int j=0;j<n;j+=(i<<1))
      for (int k=0;k<i;k++){
    int x=a[j+k],y=(ll)a[j+k+i]*w[r][num/(i<<1)*k]%P;
    a[j+k]=(x+y)%P; a[j+i+k]=(x+P-y)%P;
      }
  if (!r) for (int i=0,inv=Pow(n,P-2);i<n;i++) a[i]=(ll)a[i]*inv%P;
}

int n,m;
int a[N];

int main(){
  freopen("QAQ_bipartite_one.in","r",stdin);
  freopen("QAQ_bipartite_one.out","w",stdout);
  scanf("%d",&n); Init(n);
  int L=0; m=1;
  while (m<=(n<<1)) m<<=1,L++;
  for (int i=1;i<m;i++) R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
  Pre(m);
  for (int i=0;i<=n;i++)
    a[i]=Inv(fac[i]*Pow(SQRT,(ll)i*i%(P-1))%P);
  FFT(m,a,1);
  for (int i=0;i<m;i++) a[i]=(ll)a[i]*a[i]%P;
  FFT(m,a,0);
  printf("%d\n",(ll)a[n]*fac[n]%P*Pow(SQRT,(ll)n*n%(P-1))%P);
  return 0;
}
QAQ_bipartite_two

设上一题的也就是有色的指数生成函数是 F(x) ，这题的无色的指数生成函数为 G(x) ，无色的强制必须连通的指数生成函数是 H(x) ，那么

G=eH

F=∑i=0∞2iHii!=e2H

因为 k 个联通块染色会有2k

所以 F=G2 直接多项式开根

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;

const int P=998244353;
const int G=3;
const int SQRT=116195171;
const int INV=(P+1)>>1;
const int N=600005;

inline ll Pow(ll a,int b){
  ll ret=1; for (;b;b>>=1,a=a*a%P) if (b&1) ret=ret*a%P; return ret;
}
inline int Root(int a){
  if (a==1) return 1;
}
inline ll Inv(ll a){
  return Pow(a,P-2);
}
ll fac[N];
inline void Init(int n){
  fac[0]=1; for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%P;
}


int num;
int w[2][N];
inline void Pre(int m){
  num=m;
  int g=Pow(G,(P-1)/m);
  w[0][0]=w[1][0]=1;
  for (int i=1;i<num;i++) w[1][i]=(ll)w[1][i-1]*g%P;
  for (int i=1;i<num;i++) w[0][i]=w[1][num-i];
}

int R[N];

inline void FFT(int *a,int n,int r){
  for (int i=0;i<n;i++) if (i<R[i]) swap(a[i],a[R[i]]);
  for (int i=1;i<n;i<<=1)
    for (int j=0;j<n;j+=(i<<1))
      for (int k=0;k<i;k++){
    int x=a[j+k],y=(ll)a[j+k+i]*w[r][num/(i<<1)*k]%P;
    a[j+k]=(x+y)%P; a[j+i+k]=(x+P-y)%P;
      }
  if (!r) for (int i=0,inv=Pow(n,P-2);i<n;i++) a[i]=(ll)a[i]*inv%P;
}


inline void GetInv(int *a,int *b,int n){
  static int tmp[N];
  if (n==1) return void(b[0]=Inv(a[0]));
  GetInv(a,b,n>>1);
  for (int i=0;i<n;i++) tmp[i]=a[i],tmp[n+i]=0;
  int L=0; while (!(n>>L&1)) L++;
  for (int i=1;i<(n<<1);i++) R[i]=(R[i>>1]>>1)|((i&1)<<L);
  FFT(tmp,n<<1,1); FFT(b,n<<1,1);
  for (int i=0;i<(n<<1);i++)
    tmp[i]=(ll)b[i]*(2+P-(ll)tmp[i]*b[i]%P)%P;
  FFT(tmp,n<<1,0);
  for (int i=0;i<n;i++) b[i]=tmp[i],b[n+i]=0;
}
const int INV2=(P+1)/2;
inline void GetRoot(int *a,int *b,int n){
  static int tmp[N],invb[N];
  if (n==1) return void(b[0]=Root(a[0]));
  GetRoot(a,b,n>>1);
  memset(invb,0,sizeof(int)*n); GetInv(b,invb,n);
  int L=0; while (!(n>>L&1)) L++;
  for (int i=1;i<(n<<1);i++) R[i]=(R[i>>1]>>1)|((i&1)<<L);
  for (int i=0;i<n;i++) tmp[i]=a[i],tmp[n+i]=0;
  FFT(tmp,n<<1,1); FFT(b,n<<1,1); FFT(invb,n<<1,1);
  for (int i=0;i<(n<<1);i++)
    tmp[i]=(ll)((ll)b[i]*b[i]+tmp[i])%P*invb[i]%P*INV2%P;
  FFT(tmp,n<<1,0);
  for (int i=0;i<n;i++) b[i]=tmp[i],b[n+i]=0;
}

int n,m;
int a[N],b[N];
ll pw[N];

int main(){
  freopen("QAQ_bipartite_two.in","r",stdin);
  freopen("QAQ_bipartite_two.out","w",stdout);
  scanf("%d",&n); Init(n);
  int L=0; m=1;
  while (m<=(n<<1)) m<<=1,L++;
  for (int i=1;i<m;i++) R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
  Pre(m<<1);
  for (int i=0;i<=n;i++)
    a[i]=Inv(fac[i]*(pw[i]=Pow(SQRT,(ll)i*i%(P-1)))%P);
  FFT(a,m,1);
  for (int i=0;i<m;i++) a[i]=(ll)a[i]*a[i]%P;
  FFT(a,m,0);
  for (int i=0;i<=n;i++) a[i]=a[i]*pw[i]%P;
  for (int i=n+1;i<m;i++) a[i]=0;
  m=1; while (m<=n) m<<=1;
  GetRoot(a,b,m);
  printf("%d",(ll)b[n]*fac[n]%P);
  return 0;
}
QAQ_bipartite_thr

F=e2H 多项式ln

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;

const int P=998244353;
const int G=3;
const int SQRT=116195171;
const int INV2=(P+1)>>1;
const int N=600005;

inline ll Pow(ll a,int b){
  ll ret=1; for (;b;b>>=1,a=a*a%P) if (b&1) ret=ret*a%P; return ret;
}
inline int Root(int a){
  if (a==1) return 1;
}
inline ll Inv(ll a){
  return Pow(a,P-2);
}
ll fac[N],inv[N];
inline void Init(int n){
  fac[0]=1; for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%P;
  inv[1]=1; for (int i=2;i<=n;i++) inv[i]=(ll)(P-P/i)*inv[P%i]%P;
  inv[0]=1; for (int i=1;i<=n;i++) inv[i]=inv[i-1]*inv[i]%P;
}


int num;
int w[2][N];
inline void Pre(int m){
  num=m;
  int g=Pow(G,(P-1)/m);
  w[0][0]=w[1][0]=1;
  for (int i=1;i<num;i++) w[1][i]=(ll)w[1][i-1]*g%P;
  for (int i=1;i<num;i++) w[0][i]=w[1][num-i];
}

int R[N];

inline void FFT(int *a,int n,int r){
  for (int i=0;i<n;i++) if (i<R[i]) swap(a[i],a[R[i]]);
  for (int i=1;i<n;i<<=1)
    for (int j=0;j<n;j+=(i<<1))
      for (int k=0;k<i;k++){
    int x=a[j+k],y=(ll)a[j+k+i]*w[r][num/(i<<1)*k]%P;
    a[j+k]=(x+y)%P; a[j+i+k]=(x+P-y)%P;
      }
  if (!r) for (int i=0,inv=Pow(n,P-2);i<n;i++) a[i]=(ll)a[i]*inv%P;
}


inline void GetInv(int *a,int *b,int n){
  static int tmp[N];
  if (n==1) return void(b[0]=Inv(a[0]));
  GetInv(a,b,n>>1);
  for (int i=0;i<n;i++) tmp[i]=a[i],tmp[n+i]=0;
  int L=0; while (!(n>>L&1)) L++;
  for (int i=1;i<(n<<1);i++) R[i]=(R[i>>1]>>1)|((i&1)<<L);
  FFT(tmp,n<<1,1); FFT(b,n<<1,1);
  for (int i=0;i<(n<<1);i++)
    tmp[i]=(ll)b[i]*(2+P-(ll)tmp[i]*b[i]%P)%P;
  FFT(tmp,n<<1,0);
  for (int i=0;i<n;i++) b[i]=tmp[i],b[n+i]=0;
}

int n,m;
int a[N],b[N],c[N],f[N];
ll pw[N];

int main(){
  freopen("QAQ_bipartite_thr.in","r",stdin);
  freopen("QAQ_bipartite_thr.out","w",stdout);
  scanf("%d",&n); Init(n);
  int L=0; m=1;
  while (m<=(n<<1)) m<<=1,L++;
  for (int i=1;i<m;i++) R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
  Pre(m<<1);
  for (int i=0;i<=n;i++)
    a[i]=Inv(fac[i]*(pw[i]=Pow(SQRT,(ll)i*i%(P-1)))%P);
  FFT(a,m,1);
  for (int i=0;i<m;i++) a[i]=(ll)a[i]*a[i]%P;
  FFT(a,m,0);
  for (int i=0;i<=n+1;i++) f[i]=(ll)a[i]*pw[i]%P*fac[i]%P,a[i]=(ll)a[i]*pw[i]%P;
  for (int i=0;i<=n;i++) c[i]=f[i+1]*inv[i]%P;
  for (int i=n+1;i<=m;i++) a[i]=0;
  GetInv(a,b,m>>1);
  FFT(b,m,1); FFT(c,m,1);
  for (int i=0;i<m;i++) b[i]=(ll)c[i]*b[i]%P;
  FFT(b,m,0);
  printf("%d\n",(ll)b[n-1]*fac[n-1]%P*INV2%P);
  return 0;
}
里阿奴摩西
关注
专栏目录
[离散化][状态压缩][集合][数论][图论][容斥原理]二分图计数
~NULL~
07-07 475
题目描述 分析 很明显是要求二分图的完备匹配数 但是n个点都与m-1个点相连，很难运用这个条件求出，但是可以想相反条件，即容斥原理 那么容易得出： F(S)=Σ(-1)^|S|*h(S1,S) S1⊆S h(S1,S)表示S1与和自己无法匹配的点匹配，S任意匹配。 都很容易求，S1只有一种匹配方案，S只用乘上另一边有多少个点即可 然后重复计算问题思考一下，(-1)^|...
生成函数，多项式算法与图的计数1
08-03
生成函数、多项式算法与图的计数是数学和计算机科学中的重要概念，尤其是在组合优化、图论和算法设计等领域。这些工具可以帮助我们系统地分析和解决复杂的问题，特别是涉及计数问题时。 生成函数是一种表示序列的...
有标号的二分图计数 [生成函数 多项式]
weixin_34337381的博客
05-03 172
有标号的二分图计数 题目也在COGS上 [HZOI 2015]有标号的二分图计数 I [HZOI 2015]有标号的二分图计数 II [HZOI 2015]有标号的二分图计数 III I 求n个点的二分图(可以不连通)的个数。\(n \le 10^5\) 其中二分图进行了黑白染色，两个二分图不同：边不同 或 点的颜色不同 水题啊，只有黑白之间连边。 \[ \sum_{k=0}^n \bino...
自定义的pow()函数（C++）
最新发布
qq_54844497的博客
09-29 535
C++中有自带的pow()函数，具有求指定底数的指定幂值。通常使用该函数求解幂，但是如果不使用原C++自带的pow()函数，有能如何呢？
有标号二分图计数
weixin_30492047的博客
04-29 275
T1 不要求联通，进行黑白染色的方案数 水题 T2 要求联通 不进行黑白染色的方案数 设S(x)为上一道题的指数生成函数，H(x)为这一道题的指数生成函数 T2 不要求联通 不进行黑白染色的方案数 把T2 exp过去即可 转载于:https://www.cnblogs.com/Creed-qwq/p/10788100.html...
COGS 2392 2393 2395 有标号的二分图计数
weixin_34015860的博客
04-22 142
有黑白关系： 枚举左部点（黑色点），然后$2^{i*(n-i)}$处理方法同：COGS 2353 2355 2356 2358 有标号的DAG计数 无关系： 发现，假设$f(i)$是一个连通块，对于一个连通块，变成无颜色的，除以二即可 由结论COGS 2353 2355 2356 2358 有标号的DAG计数：G,F为EGF，$G=ln F$ 所以方案就是：$e^{\frac{lnF}{...
COGS 有标号的二分图计数系列
afghbv0316的博客
03-22 198
其实这三道题都是不错的……(虽然感觉第三题略套路了……) 分别写一下做法好了…… COGS2392 有标号的二分图计数 I 这个就很简单了，Noip难度。 显然可以直接认为黑点和白点分别位于二分图两侧，枚举二分图左侧的点数，如果左侧的点数为$k$，那么就有$C_n^k$种选择方案，并且有$k(n-k)$条边可选，因为每条边都可以选或不选，因此答案就是 \begin{align}...
光子计数位置灵敏探测器畸变多项式校正
02-09
采用多项式校正方法对光子计数位置灵敏探测器成像畸变进行校正。介绍光子计数位置灵敏探测器工作原理并分析其畸变产生原因；介绍多项式校正原理，并给出光子计数位置灵敏探测器畸变多项式校正流程；采用该方法对两种...
jisuanqi.rar_一元 稀疏多项式 简单 计数_一元 稀疏多项式 简单计算器_简单计算器
09-22
在数学中，多项式通常表示为一系列项的和，其中每个项由一个系数和一个变量的幂次组成。当多项式中有许多系数为零时，我们称其为稀疏多项式。对于这样的多项式，存储所有零系数会浪费大量的空间，因此在编程实现时，...
有标号的二分图计数问题123
weixin_30417487的博客
03-11 205
因为OJ挂了, 所以三个代码实际上都是假代码, 应该不能过题. 有标号的二分图计数1 /* 求n个点的二分图(可以不连通)的个数。n ≤10^5 其中二分图进行了黑白染色，两个二分图不同：边不同 或 点的颜色不同 然后答案就是枚举有几个点是黑的， 然后黑白点之间连边即可 */ #include<cstdio> #include<algorithm> #include&...
[前缀和][容斥原理] Jzoj P4196 二分图计数
weixin_30399821的博客
07-06 132
Description Input Output Sample Input 样例11 20样例23 3012 Sample Output 样例11样例270 Data Constraint 题解 求二分图完备匹配的个数，左边每个点只和右边一个点不相连 ...
JZOJ4196. 二分图计数(容斥)
puppywolf的博客
07-06 339
题意： 二分图，左边n个点，右边m个点 左边每个点都与右边某一点不相连 令S表示选择左边点的状态集合，g（S）表示选择里面的点，且里面的点与右边某一点连接，且右边的点最多只能跟一个左边点相连，总共的方案数 求∑2n−1s=1S∗G(S)∑s=12n−1S∗G(S)\sum _{s=1}^{2^n-1} S*G(S) 想法： 比赛看懂题意，想到容斥...
2018.07.06【2018提高组】模拟B组 二分图计数
RainbowCrown_CGH's blog
07-06 234
From jzoj4196 Description Input Output Sample Input 样例1 1 2 0 样例2 3 3 0 1 2 Sample Output 样例1 1 样例2 70 Data Constraint 首先，这题题目大意真的很难看懂，因为这图片里面很多错字。（非良心出题人） 比如：“其中...
[HDU 5729] Rigid Frameworks （二分图联通方案计数）
MisDeer's blog
07-24 838
HDU - 5729 给定一个 N*M的矩形格子，问在一个格子里只添加一条斜线， 将整个图形变成一个刚性的图形的方案数有多少种 根据 Project Euler 434，问题可以转化为 一个左边有 NN个点，右边有 MM个点的二分图的联通方案计数 首先对于一个 N∗MN*M的格子图，每个格子要么不加， 要么向左加一条斜线，要么向右加一条斜线，所以总方案数 2N∗M2^{N*M} 接
cogs [HZOI 2015]有标号的二分图计数
weixin_30700977的博客
12-03 114
题目分析 n个点的二分染色图计数 很显然的一个式子 \[ \sum_{i=0}^n\binom{n}{i}2^{i(n-i)} \] 很容易把\(2^{i(n-i)}\)拆成卷积形式，前面讲过，不再赘述。 n个点的二分图计数 设\(f_n\)表示n个点的二分染色图个数。 设\(g_n\)表示n个点的二分连通图个数。 设\(h_n\)表示n个点的二分图个数。 分别构造f,g,h的EGF\(F,G,H...
jzoj4196. 二分图计数
某司机的黑车
07-08 241
题目描述 70% 看到这种限制条件特别少的题目，首先考虑用容斥来解 （然而我考试时就想到了但是每写出来） 我以为40min写容斥很充足 先枚举选了的点集i，之后在枚举这些点中那些选了不能选的 不能选的点被固定，之后从剩余的点数往下乘到(m-点集大小+1) 要预先处理出那些情况是不可能的，即选了的点不重复 乘积和状态中1的个数都可以预处理，所以复杂度是O(22n)O(22n)...
【JZOJ4196】二分图计数（容斥+dfs）
qq_36551189的博客
07-06 474
Problem Hint Solution 　　看到这道题时我也想了一会容斥，只不过当时比较傻逼，想到的是2n2n2^n暴枚一个S，再2n2n2^n暴枚一个S1⊆S，然后就误以为会gg。。。竟然没有想到直接O(3n)O(3n)O(3^n)枚举。。。 　　下面是TJ： 　　这个容斥应该很好理解。当S1为空集时，S-S1=S，则S中所有点都可以任意匹配；这样可能会使某个点...
二分图染色与匈牙利算法：学习笔记解析
1. 定义：匈牙利算法是由Edmonds提出的一种在多项式时间内解决二分图最大匹配问题的算法。 2. 基本思路：通过寻找增广路径的方式逐步增加匹配边数，直到找到最大匹配。 3. 实现步骤： a. 从任意一边开始，尝试为每...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值