信息学奥赛一本通 2012：【20NOIP提高组】排水系统

2012：【20NOIP提高组】排水系统

时间限制: 1000 ms 内存限制: 524288 KB
提交数: 287 通过数: 30
【题目描述】
对于一个城市来说，排水系统是极其重要的一个部分。

有一天，小 C 拿到了某座城市排水系统的设计图。排水系统由 nn 个排水结点（它们从 1∼n1∼n 编号）和若干个单向排水管道构成。每一个排水结点有若干个管道用于汇集其他排水结点的污水（简称为该结点的汇集管道），也有若干个管道向其他的排水结点排出污水（简称为该结点的排出管道）。

排水系统的结点中有 mm 个污水接收口，它们的编号分别为 1,2,⋅⋅⋅,m1,2,···,m，污水只能从这些接收口流入排水系统，并且这些结点没有汇集管道。排水系统中还有若干个最终排水口，它们将污水运送到污水处理厂，没有排出管道的结点便可视为一个最终排水口。

现在各个污水接收口分别都接收了 11 吨污水，污水进入每个结点后，会均等地从当前结点的每一个排出管道流向其他排水结点，而最终排水口将把污水排出系统。

现在小 C 想知道，在该城市的排水系统中，每个最终排水口会排出多少污水。该城市的排水系统设计科学，管道不会形成回路，即不会发生污水形成环流的情况。

【输入】
第一个两个用单个空格分隔的整数 nn，mm。分别表示排水结点数与接收口数量。

接下来 nn 行，第 ii 行用于描述结点 i 的所有排出管道。其中每行第一个整数 didi 表示其排出管道的数量，接下来 didi 个用单个空格分隔的整数 a1a1, a2a2, · · · , adiadi 依次表示管道的目标排水结点。

保证不会出现两条起始结点与目标结点均相同的管道。

【输出】
输出若干行，按照编号从小到大的顺序，给出每个最终排水口排出的污水体积。其中体积使用分数形式进行输出，即每行输出两个用单个空格分隔的整数 pp，qq，表示排出的污水体积为 pqpq。要求 pp 与 qq 互素，q=1q=1 时也需要输出 qq。

【输入样例】
5 1
3 2 3 5
2 4 5
2 5 4
0
0
【输出样例】
1 3
2 3
【提示】
【样例 1 解释】

11 号结点是接收口，4、54、5 号结点没有排出管道，因此是最终排水口。

11 吨污水流入 11 号结点后，均等地流向 2、3、52、3、5 号结点，三个结点各流入 1313 吨污水。

22 号结点流入的 1313吨污水将均等地流向 4、54、5 号结点，两结点各流入 1616 吨污水。

33 号结点流入的 1313 吨污水将均等地流向 4、54、5 号结点，两结点各流入 1616 吨污水。

最终，44 号结点排出 16+16=1316+16=13 吨污水，55 号结点排出13+16+16=2313+16+16=23 吨污水。

【数据范围】

测试点编号 n ≤ m ≤
1 ∼ 3 10 1
4 ∼ 6 1000
7 ∼ 8 105105
9 ∼ 10 10
对于所有测试点，保证 1≤n≤1051≤n≤105，1≤m≤101≤m≤10，0≤di≤50≤di≤5。

数据保证，污水在从一个接收口流向一个最终排水口的过程中，不会经过超过 1010个中间排水结点（即接收口和最终排水口不算在内）。

【下载】

下载附加文件

开始蛙蚌住了，但真的不难
废话不多代码奉上
#include
#include
#include
#include
#include
#define ll long long
#define ldb long double
#define N 200005
#define eps 1e-6
using namespace std;
ll f[N5][3],q[N],du[N],len[N],ans[N];
queue h;
ll i,j,k,m,n,o,p,l,s,t;
struct node{
ldb x,y;
}a[N],ad,c;
void read(ll &x)
{
char ch=getchar();x=0;
while (ch<‘0’||ch>‘9’) ch=getchar();
while (ch>=‘0’&&ch<=‘9’) x=(x<<3)+(x<<1)+ch-48,ch=getchar();
}
void insert(ll x,ll y) {f[++t][1]=y,f[t][2]=q[x],q[x]=t;}
ldb mod(ldb x,ldb y) {return fmod(x,y);}
ldb gcd(ldb x,ldb y) {return (mod(x,y)>eps)?gcd(y,mod(x,y)):y;}
node change(ldb x,ldb y) {ll gys=gcd(x,y);return (node){x/gys,y/gys};}
node operator +(const node &a,const node &b) {ll gys=gcd(a.y,b.y);c.y=a.y/gysb.y;c.x=c.y/a.ya.x+c.y/b.yb.x;return c;}
int main()
{
read(n),read(m);
for (i=1;i<=n;i++)
{
read(len[i]);
for (j=1;j<=len[i];j++) read§,insert(i,p),du[p]++;
if (!len[i]) ans[++ans[0]]=i;
}
for (i=1;i<=n;i++)
if (!du[i]) h.push(i),a[i]=(node){1,1};
else a[i]=(node){0,1};
while (!h.empty())
{
ll st=h.front();h.pop();
if (!len[st]) continue;
ad=change(a[st].x,a[st].y*len[st]);
for (ll k=q[st];k;k=f[k][2])
{
a[f[k][1]]=a[f[k][1]]+ad;
du[f[k][1]]–;
if (!du[f[k][1]]) h.push(f[k][1]);
}
}
for (i=1;i<=ans[0];i++) a[ans[i]]=change(a[ans[i]].x,a[ans[i]].y),printf("%.0LF %.0LF\n",a[ans[i]].x,a[ans[i]].y);
return 0;
}