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

广告:

#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 105
using 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;
}

你可能感兴趣的:(动态规划,time,拓扑图,Meeting,BZOJ3890)