老师开始帮我们查漏补缺啦!
我们的老师这两天给了我们一些我们没怎么学的函数和算法,比如STL的函数和拓扑排序之类的,这个题就是讲拓扑排序的。
先看题板:
题目背景 你知道食物链吗?Delia 生物考试的时候,数食物链条数的题目全都错了,因为她总是重复数了几条或漏掉了几条。于是她来就来求助你,然而你也不会啊!写一个程序来帮帮她吧。 题目描述 给你一个食物网,你要求出这个食物网中最大食物链的数量。 (这里的“最大食物链”,指的是生物学意义上的食物链,即最左端是不会捕食其他生物的生产者,最右端是不会被其他生物捕食的消费者。) Delia 非常急,所以你只有 11 秒的时间。 由于这个结果可能过大,你只需要输出总数模上 80112002 的结果。 输入格式 第一行,两个正整数 n、m,表示生物种类 n 和吃与被吃的关系数 m。 接下来 mm 行,每行两个正整数,表示被吃的生物A和吃A的生物B。 输出格式 一行一个整数,为最大食物链数量模上 80112002 的结果。 输入输出样例 输入 #1复制 5 7 1 2 1 3 2 3 3 5 2 5 4 5 3 4 输出 #1复制 5 说明/提示 各测试点满足以下约定:
n<=5000
嗯,看起来不错。直接暴力找出没有天敌的动物,然后操作就可以了。
我们先看辅助斯烤的关系图(仅限样例)
题目描述告诉我们,食物链的开始点是没有天敌的动物。我们就枚举一遍所有点的入度,看见有一个是0的,就把他插进队列。
然后我们从他吃的食物里找一个,把这个食物的天敌数量-1,把食物的路径变量+=自己的路径变量(一个节点的路径变量是指:以他为末尾,有多少种包含他的食物链),如果这个食物的天敌数变成0了,就说明他的路径变量完整了。这样就可以继续搜索了。
如果一个节点没有食物就饿死了就说明他是结尾了,而且他的天敌数如果同时是0的话,就表明已经可以统计食物链条数了。(已经知道以他为末尾,有多少种包含他的食物链,而且这个节点没有食物,那就结束了。)
好了,需要主要思想已经讲完了,看看代码吧:
#include<iostream> #include<cstdio> #include<queue> using namespace std; int a,b,n,m,bj[5005][5005],zshu,shu;//对了,用longlong的话空间会炸掉的哦(longlong是8个字节,比int多一倍) struct hehe { int as,bs,qz; }sz[5005];//sz[i]表示i的情况(比如天敌有几个,食物有几个,强制末尾是他的食物链有几个) queue<int>q; int main() { scanf("%d%d",&n,&m); for(int i=0;i<m;i++) { scanf("%d%d",&a,&b); sz[a].bs++;//bs是食物数 sz[b].as++;//as是天敌数 bj[a][b]=1;//bj[i][j]=1表示i会吃掉j } for(int i=1;i<=n;i++)//查找食物链的开始节点 { if(sz[i].as==0)//找到啦 { sz[i].qz=1; q.push(i); } } while(q.empty()!=true)//只要还有一个节点,while就休想退出! { shu=q.front(); for(int i=1;i<=n;i++) { if(bj[shu][i]==1)//研究表明:第shu号动物会吃掉第i号动物 { sz[i].as--;//把他的天敌-1(就是还需要收集多少个节点的路径数) sz[i].qz+=sz[shu].qz;//qz是路径总数的意思。 sz[i].qz%=80112002;//题目告诉我们要取模(这里不取会40分,别问我怎么知道的) if(sz[i].as==0)//符合要求,要不是结尾,要不是下一个。 { if(sz[i].bs==0)//他是结尾 { zshu+=sz[i].qz;//zshu是食物链的总数 zshu%=80112002;//正常操作,取模 continue;//他都结尾了加他干嘛? } q.push(i);//他不是结尾,加进去继续搜 } } } q.pop();//删除第一个,也就是现在这个 } cout<<zshu<<endl;//输出 return 0; }
这个题就很健康的讲完了,感觉还行。(感觉拓扑排序不难的样子,希望是这样)