bzoj4011 [HNOI2015]落忆枫音（拓扑序dp+容斥原理+朱刘算法）

最新推荐文章于 2019-08-21 18:03:00 发布

Icefox_zhx

最新推荐文章于 2019-08-21 18:03:00 发布

阅读量378

点赞数

CC 4.0 BY-SA版权

分类专栏： bzoj ---------dp----------- ---------图论--------- 拓扑排序容斥原理

本文链接：https://blog.youkuaiyun.com/Icefox_zhx/article/details/79227930

bzoj 同时被 3 个专栏收录

764 篇文章

订阅专栏

---------dp-----------

304 篇文章

订阅专栏

---------图论---------

260 篇文章

订阅专栏

本文详细介绍了一个基于图论的算法实现，通过使用C++语言进行编码，解决了特定节点间路径问题。该算法首先定义了图的基本结构，并实现了节点间的连接与遍历。通过对特定节点的特殊处理，最终求解出除了指定路径外的所有可能路径的组合数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

大爷题解传送：http://blog.youkuaiyun.com/popoqqq/article/details/45194103

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
#define mod 1000000007
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,X,Y,in[N],du[N],ans=1,f[N],inv[N],h[N],num=0;
struct edge{
    int to,next;
}data[N<<1];
queue<int>q;
inline void add(int x,int y){
    data[++num].to=y;data[num].next=h[x];h[x]=num;in[y]++;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();X=read();Y=read();if(Y!=1) add(X,Y);
    for(int i=1;i<=m;++i){
        int x=read(),y=read();add(x,y);
    }memcpy(du,in,sizeof(in));
    for(int i=2;i<=n;++i) ans=(ll)ans*du[i]%mod;
    if(X==Y||Y==1){printf("%d\n",ans);return 0;}inv[1]=1;
    for(int i=2;i<=n;++i) inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
    f[Y]=(ll)ans*inv[du[Y]]%mod;
    q.push(1);q.push(Y);
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(y==Y) continue;
            f[y]=(f[y]+(ll)f[x]*inv[du[y]])%mod;
            if(--in[y]==0) q.push(y);
        }
    }ans-=f[X];if(ans<0) ans+=mod;else ans%=mod;
    printf("%d\n",ans);
    return 0;
}