这道题目题面真长，废话一堆。
另外：这大概是我第一道独立做出来的HNOI2011年以后的题目了吧。像我水平这么差的都能做出来，dalao您不妨试一下自己想想？
题目大意：给一个DAG，其中1号点没有入度，现在新加入一条不重复的边，使得它可能有环。求它的生成子图个数，使得子图正好包含N-1条边且1号点与其它的所有点连通。
题目分析：
我们首先要发现这是一个树的结构！有向的树。
分析树的特点，树的父亲只有一个，我们不妨从这里入手。
在这一个生成子图中，谁是谁的父亲？
我们知道1号点一定是root，这是为什么呢？
原因在于一号点没有入度，我们就认为了它是root。
那么假设有一条x->y的路径，那么x就是y的父亲，这样的定义是合理的。
考虑不加边之前的情况，图是一个DAG，我们求它的生成树(姑且这么叫吧)。
根据我们上面的分析，我们猜想：除一号点外所有点的入度乘积为该图的生成树数量。
证明如下：
首先我们考虑这个是怎么来的。每个点选择一个父亲，则根据乘法原理得到上面的猜想。
那么我们的命题是：每个点选一个父亲，则这个方案必定合法。
我们一共对N-1个点选了到fa的边，每次选择都合并了两个点，共合并了N-1次，因此共N个点被合并，又由于不存在环的情况，所以这样的方案是合法的。
考虑多了一条边的情况。那么我们由于已经求得了所有不包含这条边的情况，我们只需要分析多了这条边的情况。
假设这是一条从x到y的边。当这条边必选的时候，y点必定不被考虑。因为它的父亲已经确定。
那么我们对剩下的点重新做一遍乘法，得到了ans2。
考虑答案多了什么。
当选出的边与当前边成环的时候，这个并不是一个生成树，我们没办法解决了吗？当然不是。
我们不妨单独拿出一个环考虑。当形成了这个环的时候，对答案的影响是多少。实际上这不难统计，把剩下的点的入度全部乘起来就是答案。嘿嘿，这不就是总和除以现在的点数吗。
那么我们考虑求y到x的路径的点的入度的逆，DP一发轻松解决。然后乘法分配率证明正确性。
后记：这题我先写了个错的然后发现我想错了，答案小了，调了调变大了，然后又调，结果搜索出了正解233。
代码：
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 102000;
const int mod = 1000000007;
int n,m,x,y;
ll ans,ans2;
int in[maxn],arr[maxn];
ll dist[maxn],inv[maxn];
vector <int> g[maxn];
vector <int> ng[maxn];
ll fast_pow(int now,int p){
if(p == 1) return now;
if(p == 0) return 1;
ll z = fast_pow(now,p/2); z*=z; z %= mod;
if(p & 1){z*=now;z%=mod;}
return z;
}
void read(){
scanf("%d%d%d%d",&n,&m,&x,&y);
for(int i=1;i<=m;i++){
int a,b; scanf("%d%d",&a,&b);
g[a].push_back(b);in[b]++;
ng[b].push_back(a);
}
for(int i=2;i<=n;i++) inv[i] = fast_pow(in[i],mod-2);
}
void dfs(int now){
if(arr[now]) return;
arr[now] = 1;
for(int i=0;i<ng[now].size();i++){
dfs(ng[now][i]);
dist[now] += (dist[ng[now][i]]*inv[now])%mod;
dist[now] %= mod;
}
}
void work(){
ans = ans2 = 1;
for(int i=2;i<=n;i++) ans = (ans*in[i]) % mod;
if(x == y){printf("%lld\n",ans);return;}
if(y == 1){printf("%lld\n",ans);return;}
for(int i=2;i<=n;i++) if(i != y) ans2 = (ans2*in[i]) % mod;
in[y]++;inv[y] = fast_pow(in[y],mod-2);dist[y] = inv[y];arr[y] = 1;dfs(x);
ans += ans2; ans %= mod;
ans -= (ans*dist[x])%mod; ans += mod; ans %= mod;
printf("%lld\n",ans);
}
int main(){
read();
work();
return 0;
}
原文地址：https://www.cnblogs.com/1-1-1-1/p/8596753.html