【BZOJ3890】【Usaco2015 Jan】Meeting Time 拓扑图简单DP

来源:互联网 发布:网络信息点测试 编辑:程序博客网 时间:2024/05/16 10:36

广告:

#include <stdio.h>int main(){    puts("转载请注明出处[vmurder]谢谢");    puts("网址:blog.csdn.net/vmurder/article/details/43971435");}

题意:

给出一个n个点m条边的有向无环图,每条边两个边权。
n<=100,没有重边。
然后要求两条长度相同且尽量短的路径,
路径1采用第一种边权,路径2采用第二种边权。
没有则输出”IMPOSSIBLE”

题解:

简单拓扑图DP。
bool型数组f[i][j]、g[i][j]表示第i个点是否有1到此点权值为j的路径1、2。

代码:

#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>#define N 105using namespace std;bool f[N][N*N],g[N][N*N];int map1[N][N],map2[N][N];int n,m,stk[N],top,d[N];int main(){//  freopen("test.in","r",stdin);    int i,j,k;    int a,b,c;    memset(map1,-1,sizeof map1);    memset(map2,-1,sizeof map2);    scanf("%d%d",&n,&m);    while(m--)    {        scanf("%d%d%d%d",&a,&b,&c,&k);        map1[a][b]=c,map2[a][b]=k,d[b]++;    }    for(i=1;i<=n;i++)if(!d[i])stk[++top]=i;    f[1][0]=g[1][0]=1;    while(top)    {        a=stk[top--];        for(i=1;i<=n;i++)if(map1[a][i]+1)        {            if(--d[i]==0)stk[++top]=i;            for(j=map1[a][i];j<=10000;j++)f[i][j]|=f[a][j-map1[a][i]];            for(j=map2[a][i];j<=10000;j++)g[i][j]|=g[a][j-map2[a][i]];        }    }    for(i=0;i<=10000;i++)    {        if(f[n][i]&g[n][i])        {            printf("%d\n",i);            return 0;        }    }    puts("IMPOSSIBLE");    return 0;}
0 0
原创粉丝点击