noip2020T1

题意:初始有m个点拥有废水为1,求入度为0的点的废水情况

解法:拓扑排序,注意分母分子呈指数级增长,注意使用__int128或高精,对于任意一个点t子节点i而言,有(ai*bt+at*bi)/bi*bt,注意边处理答案贡献,边简化分数,防止爆128。当入度为0是,将点i放进队列中,注意每一次对于一个点i而言,都要更新一遍贡献

#include<bits/stdc++.h>
#define int __int128
using namespace std;

const int N = 1e5 + 10;
int n,m,cnt[N],gcnt[N];
bool g[N];
struct num{
	int fz,fm;
}dis[N];
struct node{
	int d,fz,fm;
};
queue<node> q;
vector<int> v[N];
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;
}
void get(int &fz,int &fm){
	int p=__gcd(fz,fm);
	fz/=p;fm/=p;
}
void up(int a1,int b1,int &a2,int &b2){
    a2=a1*b2+a2*b1;b2*=b1;
    get(a2,b2);
}
int a_get(int a,int b){
	if(!a)  return b;
	if(!b)  return a;
	else return a*b/__gcd(a,b);
}
void print(int x){
	if(x>9)  print(x/10);
	putchar(x%10+'0');
}
void print(){
	for(int i=1;i<=n;i++)  
	  if(g[i])
	    get(dis[i].fz,dis[i].fm);
	for(int i=1;i<=n;i++)
	  if(g[i]){
	  	print(dis[i].fz);
	  	cout<<" ";
		print(dis[i].fm);
		cout<<endl;
	  }
}
signed main(){
	n=read();m=read();
	for(int i=1,k;i<=n;i++){
		k=read();
		if(!k)  g[i]=true;
		for(int j=1,op;j<=k;j++){
			op=read();v[i].push_back(op);
			++cnt[op];
		}
	}
	
	for(int i=1;i<=m;i++){
		q.push({i,1,1});dis[i]={1,1};
	}
	while(q.size()){
		node t=q.front();
		q.pop();
	
		for(int i:v[t.d]){
			cnt[i]--;
			if(!dis[i].fz){
				dis[i].fz=t.fz;dis[i].fm=dis[t.d].fm*v[t.d].size(); 
			}
			else{
				up(t.fz,t.fm*v[t.d].size(),dis[i].fz,dis[i].fm);
			}
			get(dis[i].fz,dis[i].fm);
			if(!cnt[i])  q.push({i,dis[i].fz,dis[i].fm});
		}
	}
	//cout<<dis[10].fz;
	print();
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值