Description
当今社会,在社交网络上看朋友的消息已经成为许多人生活的一部分。通常,一个用户在社交网络上发布一条消息
(例如微博、状态、Tweet等)后,他的好友们也可以看见这条消息,并可能转发。转发的消息还可以继续被人转
发,进而扩散到整个社交网络中。在一个实验性的小规模社交网络中我们发现,有时一条热门消息最终会被所有人
转发。为了研究这一现象发生的过程,我们希望计算一条消息所有可能的转发途径有多少种。为了编程方便,我们
将初始消息发送者编号为1,其他用户编号依次递增。该社交网络上的所有好友关系是已知的,也就是说对于A、B
两个用户,我们知道A用户可以看到B用户发送的消息。注意可能存在单向的好友关系,即A能看到B的消息,但B不
能看到A的消息。
还有一个假设是,如果某用户看到他的多个好友转发了同一条消息,他只会选择从其中一个转发,最多转发一次消
息。从不同好友的转发,被视为不同的情况。
如果用箭头表示好友关系,下图展示了某个社交网络中消息转发的所有可能情况。
![](https://www.lydsy.com/JudgeOnline/upload/201804/v1%281%29.png)
初始消息是用户1发送的,加粗箭头表示一次消息转发
Input
输入文件第一行,为一个正整数n,表示社交网络中的用户数:
第二行为一个正整数m.表示社交网络中的好友关系数目。
接下来m行,每行为两个空格分隔的整数ai和bi,表示一组好友关系,即用户ai可以看到用户bi发送的消息。
1≤n≤250,1≤ai,bi≤n,1≤m≤n(n-1)
Output
输出文件共一行,为一条消息所有可能的转发途径的数量,除以1 0007所得的余数。
Sample Input
4
7
2 1
3 1
1 3
2 3
3 2
4 3
4 2
7
2 1
3 1
1 3
2 3
3 2
4 3
4 2
Sample Output
6
Solution
一眼矩阵树定理模板题,省选出这个真的。
用出度矩阵和临接矩阵相减构造基尔霍夫矩阵,求矩阵行列式。
Code
#include <iostream>
#include <cstdio>
#define R register
#define ll long long
using namespace std;
namespace Dntcry
{
inline int read()
{
R int a = 0, b = 1; R char c = getchar();
for(; c < '0'|| c > '9'; c = getchar()) (c == '-') ? b = -1 : 0;
for(; c >= '0' && c <= '9'; c = getchar()) a = (a << 1) + (a << 3) + c - '0';
return a * b;
}
const int Mod = 10007;
int n, m, Map[310][310], inv[10010], Ans;
int Guass()
{
R int res = 1;
for(R int i = 2, j; i <= n; i++)
{
for(j = i; !Map[j][i]; j++) ;
if(j > n) continue ;
if(j != i)
{
res = Mod - res;
for(R int k = i; k <= n; k++) swap(Map[i][k], Map[j][k]);
}
for(j = i + 1; j <= n; j++)
if(Map[j][i])
{
R int tmp = Map[j][i] * inv[Map[i][i]] % Mod;
for(R int k = i; k <= n; k++)
Map[j][k] = (Map[j][k] - tmp * Map[i][k] % Mod + Mod) % Mod;
}
}
for(R int i = 2; i <= n; i++) res = res * max(Map[i][i], 1) % Mod;
return res;
}
int Main()
{
n = read(), m = read();
inv[0] = inv[1] = 1;
for(R int i = 2; i < 10010; i++) inv[i] = (Mod - Mod / i) * inv[Mod % i] % Mod;
for(R int i = 1, u, v; i <= m; i++)
{
v = read(), u = read();
Map[v][v] = (Map[v][v] + 1) % Mod;
Map[u][v] = (Map[u][v] - 1) % Mod;
}
Ans = Guass();
printf("%d\n", Ans);
return 0;
}
}
int main() { return Dntcry :: Main(); }