NOIp TG 2003 神经网络

stp.1 对于题目的理解

除第一级外，每一个单元都接受来自上一级的信号，同时根据数学公式处理收到的信号，决定发出的信号；
给定初始状态，要求出最终状态

stp.2 主体思路：拓扑排序

先把状态已给出的点加入队列，对于队列的元素，它们发出的边所到达的点累加权值，但减少入度。如果某点入度为零则进入队列。对于拓扑排序来说，队列实现显然要比栈实现方便

stp.3 Codes

#include <iostream>
#include <fstream>
#include <cstdio>
#include <algorithm>
#include <string.h>
using namespace std;
const int maxn = 101;
const int maxp=maxn*maxn;
short int n,en=0,head[maxn],ind[maxn],outd[maxn],x,y,l,qans[maxn],tans=0,q[maxn],h=0,t=0,vis[maxn];
int u[maxn],p,c[maxn];
struct edge {
    int fr,gt,len,nxt;
};
edge e [maxp];
int qread() {
    int x=0,f=1;
    char ch=getchar();
    while (ch>'9'||ch<'0') {
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9') {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
void adg(int x,int y,int len) {
    ++en;
    e[en].fr=x;
    e[en].gt=y;
    e[en].len=len;
    e[en].nxt=head[x];
    head[x]=en;
    ++outd[x];
    ++ind[y];
}
int main() {
   // freopen ("1308.in","r",stdin);
    n=qread();
    p=qread();
 //   memset(c,0,sizeof(c));      
    for (int i=1; i<=n; i++) {
        c[i]=qread();
        u[i]=qread();
    }


    memset(head,0,sizeof(head));
    memset(ind,0,sizeof(ind));
    memset(outd,0,sizeof(outd));
    memset(vis,0,sizeof(vis));
    for (int i=1; i<=p; i++) {
        x=qread();
        y=qread();
        l=qread();
        adg(x,y,l);
    }
    for (int i=1; i<=n; i++) {
        //  cout<<i<<' '<<ind[i]<<' '<<outd[i]<<endl;
        if(!outd[i])qans[++tans]=i;

        if(!ind[i])q[++t]=i;
    }
    while (h<t) {
        int k=q[++h];
        for (int i=head[k];i;i=e[i].nxt){
        //  cout<<c[k]<<endl;
            c[e[i].gt]+=e[i].len*c[k];
            --ind[e[i].gt];

            if(ind[e[i].gt])continue;
        //  cout<<c[e[i].gt]<<endl;    
            if(c[e[i].gt]>=u[e[i].gt]){
                q[++t]=e[i].gt;
                c[e[i].gt]-=u[e[i].gt];
            }
        }
    }
    int cnt=0;
    for (int i=1;i<=tans;i++){
        if(c[qans[i]]>0){
            ++cnt;
            cout<<qans[i]<<' '<<c[qans[i]]<<endl;
        }
    }
    if(!cnt)cout<<"NULL"<<endl;
    //for (int i=1;i<=n;i++)cout<<outd[i]<<' ';
    fclose (stdin);
    fclose (stdout);
    return 0;
}