BZOJ4562 [Haoi2016]食物链

最新推荐文章于 2019-07-25 15:21:00 发布

原创最新推荐文章于 2019-07-25 15:21:00 发布 · 1.1k 阅读

1 ·

CC 4.0 BY-SA版权

DP、递推专栏收录该内容

84 篇文章

订阅专栏

本文介绍了一种使用SPFA算法解决特定食物链计数问题的方法。通过定义f[x]表示以x为结尾的食物链数量，g[x]表示未完全确定的食物链数量，利用SPFA算法进行遍历计算。代码实现包括了图的构建、节点处理等关键步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

f[x]表示到x为结尾的食物链条数，g[x]表示还没给出边传过去的食物链条数，拿spfa跑一遍即可，其实直接记忆化搜索就行，当时想到了BZ上拿魔法或者物理打怪那题，就傻逼一样写了个用spfa跑的递推-_-

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<iomanip>
#include<ctime>
#include<cmath>
#include<map>
#include<bitset>
#include<stack>
#include<queue>
#include<set>
#include<vector>
using namespace std;
#define ll long long 
#define MAXN 200010
#define MAXM 200010
#define INF 2000000000
#define MOD 2000000007
#define QMOD 200000
struct vec{
    int to;
    int fro;
};
vec mp[MAXM];
int tai[MAXN],cnt;
int n,m;
int d[MAXN];
int q[MAXN*10],h,t;
int f[MAXN],g[MAXN];
bool iq[MAXN];
int ans;
inline void be(int x,int y){
    mp[++cnt].to=y;
    mp[cnt].fro=tai[x];
    tai[x]=cnt;
}
int main(){
    int i;
    int ix,iy;
    scanf("%d%d",&n,&m);
    for(i=1;i<=m;i++){
        cin>>ix>>iy;
        be(ix,iy);
        d[iy]++;
    }
    for(i=1;i<=n;i++){
        if(!d[i]){
            q[t++]=i;
            if(tai[i]){
                f[i]=g[i]=1;
            }
            iq[i]=true;
            t%=QMOD;
        }
    }
    int x,y;
    while(h!=t){
        x=q[h++];
        h%=QMOD;
        iq[x]=false;
        for(i=tai[x];i;i=mp[i].fro){
            y=mp[i].to;
            f[y]+=g[x];
            g[y]+=g[x];
            if(!iq[y]){
                q[t++]=y;
                t%=QMOD;
                iq[y]=true;
            }
        }
        g[x]=0;
    }
    for(i=1;i<=n;i++){
        if(!tai[i]){
            ans+=f[i];
        }
    }
    printf("%d\n",ans);
    return 0;
}
   
/*
4
a e
b e
c e
d e
   
   
   
   
*/