首页 > 代码库 > hdu4780 费用流 (机器任务工作不中断问题)

hdu4780 费用流 (机器任务工作不中断问题)

机器,任务 ,每个任务有有时间,不可中断。

题意:m个机器,n个糖果要加工,给出每个糖果的工作时间(s,t),以及糖果之间、机器预备时间以及费用,求最小费用。

这题开始受原来可以时间中断那题影响,开始用时间建图,巨麻烦,后来学习了,才觉悟时间只是计算费用的,没有帮毛钱关系, s-->机器-》糖果-》t;

因为要每个糖果都加工一次,糖果拆点,必经过(-inf)。有一个没过,则无解。

网上的添加一组新点法也可以。我的使用释放压力+无穷碧流法也好。

跪点:注意!用inf=0x3f3f3f3f,最后加上n*inf,爆int! 这种以后要注意!顺变 e[][3]也要longlong !

本来这样我就A了!但是竟然因为一个变量的命名搞混了WA一天!!!SB了!  自己搞变量不要搞这么相似啊 !!!!!!!!!!


#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
const long long inf=0x3f3f3f3f;
const int maxn=505,maxe=80000;
int n,m,k;
int ss,tt;
int head[maxn];long long e[maxe][4];int nume=0;
void inline adde(int i, int j, int c,long long w)
{
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;
    e[nume][2]=c;e[nume++][3]=w;
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;
    e[nume][2]=0;e[nume++][3]=-w;
}
int inq[maxn];long long d[maxn];int pre[maxn];int prv[maxn];
bool spfa(long long &sum)
{
    for(int i=0;i<=tt;i++)
       {
           inq[i]=0;d[i]=inf;
       }
    queue<int>q;
    q.push(ss);
    inq[ss]=1;d[ss]=0;
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
        inq[cur]=0;
        for(int j=head[cur];j!=-1;j=e[j][1])
        {
            int v=e[j][0];
            if(d[v]>d[cur]+e[j][3]&&e[j][2]>0)
            {
                d[v]=d[cur]+e[j][3];
                pre[v]=j;
                prv[v]=cur;
                if(!inq[v])
                {
                    q.push(v);
                    inq[v]=1;
                }
            }
        }
    }
    if(d[tt]==inf)return 0;
   int minf=inf;
    int cur=tt;
    while(cur!=ss)
    {
        minf=minf<e[pre[cur]][2]?minf:e[pre[cur]][2];
        cur=prv[cur];
    }
     cur=tt;
    while(cur!=ss)
    {
       e[pre[cur]][2]-=minf;
       e[pre[cur]^1][2]+=minf;
        cur=prv[cur];
    }
    sum+=d[tt]*minf;
    return 1;
}
long long mincost()
{
    long long sum=0;
    while(spfa(sum));
    return sum;
}
int mac[105][2];int cmtime[105][105];int cmony[105][105];
int cgtime[105][105];int cgmony[105][105];
void build()
{
    for(int i=1;i<=m;i++)
    {
        adde(ss,i,1,0);
        adde(i,tt,1,0);
       for(int j=1;j<=n;j++)
       {
           if(cmtime[j][i]<mac[j][1])
           {
               if(cmtime[j][i]<=mac[j][0])
               adde(i,j+m,1,cmony[j][i]);
               else
               adde(i,j+m,1,cmony[j][i]+k*(cmtime[j][i]-mac[j][0]));
           }
       }
    }
    for(int i=1;i<=n;i++)
     {
         adde(i+m,i+m+n,1,-inf);
         adde(i+m+n,tt,1,0);
        for(int j=1;j<=n;j++)
        {
         if(i!=j&&mac[j][1]>mac[i][1]+cgtime[i][j])
         {
             if(mac[i][1]+cgtime[i][j]<=mac[j][0])
             adde(i+m+n,j+m,1,cgmony[i][j]);
             else
             adde(i+n+m,j+m,1,cgmony[i][j]+k*(mac[i][1]+cgtime[i][j]-mac[j][0]));
         }
        }
     }
}
void init()
{
    nume=0;
    memset(head,-1,sizeof(head));
    ss=0,tt=m+n+n+1;
}
int main()
{

    while(~scanf("%d%d%d",&n,&m,&k)&&(n||m||k))
    {
        init();
        for(int i=1;i<=n;i++)
             scanf("%d%d",&mac[i][0],&mac[i][1]);

        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
         {
             scanf("%d",&cmtime[i][j]);
         }
         for(int i=1;i<=n;i++)
           for(int j=1;j<=m;j++)
         {
             scanf("%d",&cmony[i][j]);
         }
         for(int i=1;i<=n;i++)
           for(int j=1;j<=n;j++)
         {
             scanf("%d",&cgtime[i][j]);
         }
        for(int i=1;i<=n;i++)
           for(int j=1;j<=n;j++)
         {
             scanf("%d",&cgmony[i][j]);
         }
        build();
        long long ans=mincost()+n*inf;
        if(ans>=inf)
         printf("-1\n");
        else
          printf("%I64d\n",ans);
    }
}



hdu4780 费用流 (机器任务工作不中断问题)