[带标号DAG计数 容斥原理 多项式求逆 多项式求ln] COGS 有标号的DAG计数系列

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

题目链接 传送门I 传送门II 传送门III 传送门IV
题解戳这里
大致解法和思想是同这里一样的

DAG1

找出一个入度为 0 的点 那么答案为C1n∗21∗(n−1)∗fn−1 但是一张DAG里会有多个入度为 0 的点会重复计数
考虑容斥 得到

fn=∑k=1n(−1)k−1∗Ckn∗2k∗(n−k)∗fn−k

直接递推求解

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

const int N=5005;
const int P=10007;
const int Root2=2641;

inline int Pow(int a,int b){
  int ret=1; for (;b;b>>=1,a=a*a%P) if (b&1) ret=ret*a%P; return ret;
}

int _pow[P],C[N][N];
inline void Pre(){
  C[0][0]=1;
  for (int i=1;i<=5000;i++){
    C[i][0]=1;
    for (int j=1;j<=i;j++)
      C[i][j]=(C[i-1][j-1]+C[i-1][j])%P;
  }
  _pow[0]=1;
  for (int i=1;i<P;i++) _pow[i]=_pow[i-1]*2%P;
}

int n;
ll f[N];

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&n); Pre();
  f[0]=1;
  for (int i=1;i<=n;i++){
    for (int j=1;j<=i;j++){
      int tmp=(ll)C[i][j]*_pow[j*(i-j)%(P-1)]*f[i-j]%P;
      if (j&1) f[i]+=tmp; else f[i]-=tmp;
    }
    f[i]=(f[i]%P+P)%P;
  }
  printf("%d\n",f[n]);
  return 0;
}

DAG2

考虑如何拆分 2k∗(n−k)

2k∗(n−k)=2n222k22∗2(n−k)22

这里我们需要求 2 的二次剩余
那么把DAG1中的柿子改为分治FFT或多项式求逆就可以了
但是在DAG1中不能这么做 因为模数不对 汗

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

const int N=600005;
const int P=998244353;
const int G=3; const int Root2=116195171;

inline int 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;
}

int num;
int w[2][N];
ll fac[N],inv[N];

inline void Pre(int n){
  num=n;
  ll g=Pow(G,(P-1)/num);
  w[1][0]=1; for (int i=1;i<num;i++) w[1][i]=g*w[1][i-1]%P;
  w[0][0]=1; 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++){
    ll x=a[j+k],y=(ll)w[r][num/(i<<1)*k]*a[j+i+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]=Pow(a[0],P-2));
  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];

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&n);
  fac[0]=1; for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%P;
  for (m=1;m<=n;m<<=1); Pre(m<<1);
  A[0]=1;
  for (int i=1;i<=n;i++){
    ll tmp=fac[i]*Pow(Root2,(ll)i*i%(P-1))%P;
    tmp=Pow(tmp,P-2);
    A[i]=i&1?P-tmp:tmp;
  }
  GetInv(A,B,m);
  int Ans=fac[n]*B[n]%P*Pow(Root2,(ll)n*n%(P-1))%P;
  printf("%d\n",Ans);
  return 0;
}

DAG3

设fn为 n 个点的DAG的个数
设gn为 n 个点的连通DAG的个数
容斥的柿子

gn=fn−∑k=1n−1Ck−1n−1∗gk∗fn−k

大力递推 

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

const int N=5005;
const int P=10007;
const int Root2=2641;

inline int Pow(int a,int b){
  int ret=1; for (;b;b>>=1,a=a*a%P) if (b&1) ret=ret*a%P; return ret;
}

int _pow[P],C[N][N];
inline void Pre(){
  C[0][0]=1;
  for (int i=1;i<=5000;i++){
    C[i][0]=1;
    for (int j=1;j<=i;j++)
      C[i][j]=(C[i-1][j-1]+C[i-1][j])%P;
  }
  _pow[0]=1;
  for (int i=1;i<P;i++) _pow[i]=_pow[i-1]*2%P;
}

int n;
ll f[N],g[N];

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&n); Pre();
  f[0]=1;
  for (int i=1;i<=n;i++){
    for (int j=1;j<=i;j++){
      int tmp=(ll)C[i][j]*_pow[j*(i-j)%(P-1)]*f[i-j]%P;
      if (j&1) f[i]+=tmp; else f[i]-=tmp;
    }
    f[i]=(f[i]%P+P)%P;
  }
  for (int i=0;i<=n;i++){
    g[i]=f[i];
    for (int j=1;j<i;j++)
      g[i]-=(ll)C[i-1][j-1]*g[j]*f[i-j];
    g[i]=(g[i]%P+P)%P;
  }
  printf("%d\n",g[n]);
  return 0;
}

DAG4

不论是由容斥的柿子还是多项式exp ln的组合意义都不难得出多项式的做法

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
typedef long long ll;

const int N=600005;
const int P=998244353;
const int G=3; const int Root2=116195171;

inline int 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;
}

int num;
int w[2][N];
ll fac[N],inv[N];

inline void Pre(int n){
  num=n;
  ll g=Pow(G,(P-1)/num);
  w[1][0]=1; for (int i=1;i<num;i++) w[1][i]=g*w[1][i-1]%P;
  w[0][0]=1; 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++){
    ll x=a[j+k],y=(ll)w[r][num/(i<<1)*k]*a[j+i+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]=Pow(a[0],P-2));
  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,F[N];
int A[N],B[N],invB[N];

inline void Calc(){
  A[0]=1;
  for (int i=1;i<=n;i++){
    ll tmp=fac[i]*Pow(Root2,(ll)i*i%(P-1))%P;
    tmp=Pow(tmp,P-2);
    A[i]=i&1?P-tmp:tmp;
  }
  cl(B); GetInv(A,B,m);
  for (int i=0;i<=n;i++)
    F[i]=fac[i]*B[i]%P*Pow(Root2,(ll)i*i%(P-1))%P;
}

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&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]=inv[P%i]*(P-P/i)%P;
  inv[0]=1; for (int i=1;i<=n;i++) inv[i]=(inv[i]*inv[i-1])%P;
  for (m=1;m<=n;m<<=1); Pre(m<<1);
  Calc();
  cl(A); cl(B);
  for (int i=1;i<=n;i++) A[i]=F[i]*inv[i-1]%P;
  for (int i=0;i<=n;i++) B[i]=F[i]*inv[i]%P;
  GetInv(B,invB,m);
  FFT(A,m<<1,1); FFT(invB,m<<1,1);
  for (int i=0;i<(m<<1);i++) A[i]=(ll)A[i]*invB[i]%P;
  FFT(A,m<<1,0);
  printf("%d\n",(ll)A[n]*fac[n-1]%P);
  return 0;
}
里阿奴摩西
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
有标号DAG计数
Freopen的博客
05-23 492
题目 设fif_ifi​为iii个点的有标号DAG\rm DAGDAG数量。 枚举有jjj个点使其入度为000的DAG\rm DAGDAG数量： gj=(ij)2j(i−j)fi−jg_j = \binom ij 2^{j(i-j)} f_{i-j}gj​=(ji​)2j(i−j)fi−j​ 设恰好有jjj个点入度为000的DAG\rm DAGDAG数量为hjh_jhj​ 有gk=∑j=ki(jk)hjg_k = \sum_{j=k}^i\binom jk h_jgk​=∑j=ki​(kj​)hj​ 二项式
【题解】有标号的DAG计数4
weixin_30267697的博客
07-28 125
[HZOI 2015] 有标号的DAG计数 IV 我们已经知道了\(f_i\)表示不一定需要联通的\(i\)节点的dag方案，考虑合并 参考【题解】P4841 城市规划(指数型母函数+多项式Ln)，然后答案\(h_i\)母函数\(H(x)\)就这样解 由于 \[ H(x)=\sum_{i=0}^{\inf} \dfrac {(F(x))^i} {i!} \] 则 \[ H(x)=e^{F(x)}...
EGF中多项式exp的组合意义
最新发布
sophilex的博客
10-24 329
仔细理解一下，众所周知，一个大小为n的排列一定可以拆成若干个环，每一个环内部的排列数就是一个圆排列的方案数，所以大小为n的排列的方案数就是把。EGF一般用来处理多重集的排列问题，在其上可以定义多项式的exp运算，在处理一类问题的时候有独特的作用。分成若干个非空集合，每一个集合的圆排列方案数之积，这与我们上面讲到的exp的组合意义相符合。的个数，这是一个只跟集合大小有关的值），那么根据生成函数的定义，我们不难得到下式。是因为这k个集合是无序的，而原本的多个多项式卷积显然是有序的。考虑大小为n的排列的个数是。
指数型生成函数（二）：贝尔数、无向连通图计数、弱连通DAG计数
hans774882968的博客
08-09 421
贝尔数：n个元素的集合的划分的个数 注意到S2的定义就是n个元素放入k个无区别的盒子的个数，“集合”是一种无区别的盒子，所以贝尔数等于一行S2的和。 定义dp[n]为所求，则考虑1号元素所在集合的元素个数： dp[n]=∑i=1ndp[n−i]∗Cn−1i−1=∑i=0n−1dp[n−1−i]∗Cn−1i,dp[0]=1 dp[n]=\sum_{i=1}^ndp[n-i]*C_{n-1}^{i-1}=\sum_{i=0}^{n-1}dp[n-1-i]*C_{n-1}^i,dp[0]=1 dp[n]=i=.
多项式对数函数（多项式求ln）小结
a_forever_dream的博客
05-19 1737
多项式对数函数 顾名思义，给出一个多项式 A(x)A(x)A(x)，求出一个多项式 G(x)G(x)G(x)，满足 G(x)=ln⁡A(x)(modxn)G(x)=\ln A(x)\pmod {x^n}G(x)=lnA(x)(modxn)。 设 F(x)=ln⁡(x)F(x)=\ln (x)F(x)=ln(x)，那么有 G(x)=F(A(x))G(x)=F(A(x))G(x)=F(A(x))，对等式两边同时求导，得到 G′(x)=F′(A(x))A′(x) G'(x)=F'(A(x))A'(x) G′(x)
随机生成有向无环图 DAG源代码示例（C语言）
01-02
随机生成有向无环图 DAG 源代码示例（C 语言） 有向无环图（Directed Acyclic Graph，DAG）是一种常见的数据结构，用于表示有向边和节点之间的关系。在计算机科学和信息技术领域，DAG 广泛应用于数据分析、机器学习...
dag:Golang中的又一个有向无环图（DAG）实现
05-28
有向无环图（DAG）的实现。 该实现是快速且线程安全的。 它可以防止添加循环或重复，从而始终保持有效的DAG。 该实现缓存后代和祖先，以加快后续调用的速度。 快速开始 跑步： package main import ( "fmt" ...
dag-solve:基于有向无环图（DAG）的求解器
05-12
达格溶 基于有向无环图（DAG）的求解器。 此处构建的图始终只有一个根节点，并且这个根节点只能有一个度数，而我没有度数。 在根节点以下，除防止循环外，没有其他限制。 可以将图中创建的每个节点配置为： 数学求解...
airflow-maintenance-dags:一系列DAG工作流程可帮助维持Airflow的运行
05-07
一系列DAG /工作流程可帮助维持Airflow的运行 DAG /工作流程 备份配置您可以部署到Airflow中的维护工作流程，以定期获取各种Airflow配置和文件的备份。 清除缺失的数据 您可以部署到Airflow中的维护工作流程，以...
dag:Haskell 中类型良好的有向无环图
06-29
import Data.Graph.DAG.Edge -- | Edges are statically defined: edges = ECons ( Edge :: EdgeValue " foo " " bar " ) $ ECons ( Edge :: EdgeValue " bar " " baz " ) $ ECons ( Edge :: EdgeValue " foo " " ...
有标号的DAG计数系列问题
Cyhlnj
12-06 1165
传送门 II 设 fif_ifi​ 表示 iii 个点的答案 那么枚举至少 jjj 个点的出度为 000 ∑j=0i(−1)j(ij)fi−j2(i−j)j=0\sum_{j=0}^{i}(-1)^j\binom{i}{j}f_{i-j}2^{(i-j)j}=0j=0∑i​(−1)j(ji​)fi−j​2(i−j)j=0 所以 fi=∑j=1i(−1)j+1(ij)fi−j2(i−j)jf_i=\...
COGS2353[HZOI 2015] 有标号的DAG计数 I
ShadyPi的IT笔记
08-22 459
原题链接：http://cogs.pro:8080/cogs/problem/problem.php?pid=2353 有标号的DAG计数 I 【题目描述】 给定一正整数n，对n个点有标号的有向无环图(可以不连通)进行计数，输出答案mod 10007的结果。 【输入格式】 一个正整数n 【输出格式】 一个数，表示答案 【样例输入】 3 【样例输出】 25 ...
有标号的DAG计数I~IV
weixin_30640769的博客
12-11 146
题解： https://www.cnblogs.com/zhoushuyu/p/10077241.html 看到这么一篇，发现挺不错的。。 转载于:https://www.cnblogs.com/yinwuxiao/p/10105089.html
【省选模拟】20/04/07 考试 （组合意义）（生成函数）（多项式）
FSYo 的博客
04-07 270
有 nnn 个独立的变量（实数）xi∈[0,ai]x_i\in[0,a_i]xi​∈[0,ai​]，求 (∑xi)m(\sum x_i)^m(∑xi​)m 的期望 直接考虑组合意义：就是 nnn 个里面任选 mmm 个乘起来，先不考虑顺序，那么我们可以枚举每个 xix_ixi​ 被选的次数，基于此构造 EGFEGFEGF Fi(x)=∑j≥0∫0aitjdtj!xt=∑j≥0aij+1(j...
BZOJ 3684 大朋友与多叉树 多项式求幂/求exp+拉格朗日反演
-Claude
09-14 2847
BZOJ 3684 大朋友与多叉树 多项式求幂/求exp+拉格朗日反演
BZOJ 3812: 主旋律 （容斥状压+DAG计数理论）
Freopen的博客
03-02 325
题意： 有多少有向边的子集删去之后整个图仍然强联通。 n&amp;lt;=15, m&amp;lt;=n * (n-1) 上一次做弱化版：有多少有向边的子集使整个图没有环。 当时的做法是。。。。。。 ...
一类划分关系和指数级生成函数，多项式exp的关系
Cyhlnj
12-08 522
划分关系 姑且这么叫着 设满足性质 AAA 的集合为 SAS_ASA​，每个元素有标号 如果 SBS_BSB​ 是由若干个 SAS_ASA​ 组成的一个大集合 设 aia_iai​ 表示大小为 iii 的 SAS_ASA​ 的个数 设 bib_ibi​ 表示大小为 iii 的 SBS_BSB​ 的个数 构造指数级生成函数 A(x)=∑i=0∞aixii!A(x)=\sum_{i=0}^{\inft...
大数据处理之-DAG计算
姜居士的居室
05-19 7404
DAG是有向无环图（Directed Acyclic Graph）的简称。在大数据处理中，DAG计算常常指的是将计算任务在内部分解成为若干个子任务，将这些子任务之间的逻辑关系或顺序构建成DAG（有向无环图）结构。 DAG在分布式计算中是非常常见的一种结构，在各个细分领域都可以看见它，比如Dryad,Flumejava和Tez，都是明确构建DAG计算模型的典型，再如流式计算的Storm等系统或机器...
多项式求ln
gigo_64 sink in my world
08-06 170
求 两边求导 所以对于g’求导，乘g的逆元，然后得到f'，再积分回去。 关于求导： 关于积分： 关于代码 #include<bits/stdc++.h> using namespace std; #define in read() #define int long long int in{ int cnt=0,f=1;char ch=0; while(!isdigit...
DAG求拓扑排序的算法
06-11
DAG的拓扑排序算法可以用来确定每个节点在拓扑排序中的位置。具体算法如下： 1. 找到DAG中没有前置节点的节点，并将其加入结果列表中。 2. 删除该节点以及其出边。 3. 重复1和2，直到所有的节点都在结果列表中。 如果DAG中存在环路，则该算法无法完成排序。 以下是一个Python实现： ```python def topo_sort(graph): # 统计每个节点的前置节点数 in_degree = {u: 0 for u in graph} for u in graph: for v in graph[u]: in_degree[v] += 1 # 初始化队列 queue = [u for u in graph if in_degree[u] == 0] # 拓扑排序 result = [] while queue: # 取出队首节点 u = queue.pop(0) result.append(u) # 删除该节点及其出边 for v in graph[u]: in_degree[v] -= 1 if in_degree[v] == 0: queue.append(v) # 如果存在环路，则排序失败 if len(result) != len(graph): return None else: return result ``` 其中，参数graph是一个字典，表示DAG中每个节点及其出边。例如，如果DAG有三个节点0、1、2，其中0指向1和2，1指向2，则可以表示为{0: [1, 2], 1: [2], 2: []}。函数返回拓扑排序后的节点列表，如果存在环路则返回None。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值