2 条题解

  • 0
    @ 2026-9-2 0:51:15

    如果说去掉高精度的话,还是一道非常好的 拓扑 排序题目。

    本篇代码无高精,请自行添加(因为太复杂了qwq)。

    拓扑排序:

    拓扑排序主要思路在一个有向无环图中,先统计出每个点的入度个数,然后将入度为0的点入队,接着把队中每个点向它的出边做一个运算(本题中是将水分流到与其相连节点),然后断边(相连的点入度-1),最后就会得出排水节点的水量。

    有不明白的同学可以看 神经网络 车站分级 旅行计划 都是很好的拓扑排序模板题。

    拓扑排序函数:

    void tp(){
    	for(int i=1;i<=n;i++)//所有入度为0的点入队(1-m)
    		if(!in[i]){
    			book[i]=1;
    			q.push(i);
    			xx[i]=1,yy[i]=1;
    		}
    	while(!q.empty()){
    		int p=q.front();
    		q.pop();
    		if(out[p])
    			continue;
    		for(int i=0;i<a[p].size();i++){
    			add(a[p][i],xx[p],yy[p]*(1ll*a[p].size()));
    			if(book[a[p][i]])
    				continue;
    			in[a[p][i]]--;
    			if(in[a[p][i]]==0){
    				book[a[p][i]]=1;
    				q.push(a[p][i]);
    			}
    		}
    	}
    	return;
    }
    

    注意几点:

    1. 如果是 vector 存边,一定不要访问排水节点的 size() 这样可能会炸,要事先存一下。

    2.本题是前 m 个点是入水口,一定要注意审好题(虽然只有前 m 个点入度为0)。


    分数处理:

    我主要运用的是通分思想:

    $\frac{a}{b}+\frac{c}{d}=\frac{ad}{bd}+\frac{cb}{bd}=\frac{ad+cb}{bd}$

    紧接着用 gcd 化简(考场差点写挂):

    ll gcd(ll x,ll y){
    	if(y==0)
    		return x;
    	return gcd(y,x%y);
    }
    

    下面是通分代码:

    void add(int u,ll x,ll y){
    	if(y==0)
    		return;
    	if(yy[u]==0){
    		xx[u]=x;
    		yy[u]=y;
    		return;
    	}
    	ll p1=xx[u]*y+yy[u]*x;
    	ll p2=yy[u]*y;
    	ll p3=gcd(p1,p2);
    	xx[u]=p1/p3;
    	yy[u]=p2/p3;
    	return;
    }
    

    注意事项:

    1.一定要判出要添加的分母是否为0,如果为0 直接赋值即可。

    2.最后输出保险在约分一下。

    全部代码(考试代码):

    #include<iostream>
    #include<cstdio>
    #include<cmath>
    #include<algorithm>
    #include<queue>
    #include<stack>
    #include<vector>
    #include<map>
    #define ll long long
    using namespace std;
    inline ll read(){
    	ll 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,in[100001],out[100001],book[100001];
    ll xx[100001],yy[100001];
    ll gcd(ll x,ll y){
    	if(y==0)
    		return x;
    	return gcd(y,x%y);
    }
    void add(int u,ll x,ll y){
    	if(y==0)
    		return;
    	if(yy[u]==0){
    		xx[u]=x;
    		yy[u]=y;
    		return;
    	}
    	ll p1=xx[u]*y+yy[u]*x;
    	ll p2=yy[u]*y;
    	ll p3=gcd(p1,p2);
    	xx[u]=p1/p3;
    	yy[u]=p2/p3;
    	return;
    }
    vector<int> a[500001];
    queue<int> q;
    void tp(){
    	for(int i=1;i<=n;i++)
    		if(!in[i]){
    			book[i]=1;
    			q.push(i);
    			xx[i]=1,yy[i]=1;
    		}
    	while(!q.empty()){
    		int p=q.front();
    		q.pop();
    		if(out[p])
    			continue;
    		for(int i=0;i<a[p].size();i++){
    			add(a[p][i],xx[p],yy[p]*(1ll*a[p].size()));
    			if(book[a[p][i]])
    				continue;
    			in[a[p][i]]--;
    			if(in[a[p][i]]==0){
    				book[a[p][i]]=1;
    				q.push(a[p][i]);
    			}
    		}
    	}
    	return;
    }
    int main()
    {
    	//freopen("water.in","r",stdin);
    	//freopen("water.out","w",stdout);
    	n=read(),m=read();
    	for(int i=1;i<=n;i++){
    		int d=read();
    		if(d==0){
    			out[i]=1;
    			continue;
    		}
    		while(d--){
    			int v;
    			v=read();
    			a[i].push_back(v);
    			in[v]++;
    		}
    	}
    	tp();
    	for(int i=1;i<=n;i++){
    		if(out[i]){
    			add(i,0,1);
    			printf("%lld %lld\n",xx[i],yy[i]);
    		}
    	}
    	return 0;
    }
    

    谢谢大家!

    • 0
      @ 2026-6-14 8:42:24

      // 拓扑排序 O(n)
      #include<bits/stdc++.h>
      #define i8 __int128
      using namespace std;
      
      int read(){
        int x=0; char c=getchar();
        while(c<'0'||c>'9') c=getchar();
        while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
        return x;
      }
      const int N=100005;
      int h[N],idx,ne[N*5],to[N*5];
      void add(int u,int v){
        to[++idx]=v;ne[idx]=h[u];h[u]=idx;
      }
      int n,m,cd[N],rd[N]; //出度,入度
      i8 fz[N],fm[N]; //分子,分母
      
      i8 gcd(i8 x,i8 y){ //求最大公约数
        i8 t;
        while(y){
          t=y;
          y=x%y;
          x=t;
        }
        return x;
      }
      void calc(i8 &v1,i8 &v2,i8 u1,i8 u2){ //求v1/v2+u1/u2
        v1=v1*u2+v2*u1;
        v2=v2*u2;
        i8 g=gcd(v1,v2);
        v1/=g;
        v2/=g;
      }
      void topo(){ //拓扑排序
        queue<int> q;
        for(int i=1;i<=m;++i){
          fz[i]=1;   //分子初值
          q.push(i); //m个入口入队
        }  
        while(!q.empty()){
          int u=q.front(); q.pop();
          for(int i=h[u];i;i=ne[i]){
            int v=to[i];
            calc(fz[v],fm[v],fz[u],fm[u]*cd[u]);
            if(--rd[v]==0) q.push(v);
          }
        } 
      }
      void write(i8 x){ //大数字输出
        if(x>9) write(x/10);
        putchar(x%10+'0');
      }
      int main(){
        n=read(),m=read();
        for(int i=1,d;i<=n;++i){
          d=read();
          fm[i]=1;  //分母初值
          cd[i]+=d; //出度
          for(int j=1,x;j<=d;++j){
            x=read();
            add(i,x);
            rd[x]++; //入度
          }
        }
        
        topo();
        for(int i=1;i<=n;++i)if(!cd[i]){ //枚举出口
          write(fz[i]); putchar(' ');
          write(fm[i]); putchar('\n');
        }
        return 0;
      }
      
      • 1

      D154 拓扑排序[NOIP2020] 排水系统

      信息

      ID
      2000
      时间
      1000ms
      内存
      512MiB
      难度
      9
      标签
      递交数
      13
      已通过
      3
      上传者