Wiki 1935(美食节-费用流动态加点/边)

最新推荐文章于 2021-01-04 19:16:04 发布

原创最新推荐文章于 2021-01-04 19:16:04 发布 · 1.4k 阅读

0 ·

CC 4.0 BY-SA版权

费用流专栏收录该内容

7 篇文章

订阅专栏

本文介绍了一个关于菜品分配的算法问题，旨在通过合理的菜品制作顺序安排，使得顾客的总等待时间达到最小。文章提供了完整的代码实现，并通过具体示例解释了算法的工作原理。

CZ市为了欢迎全国各地的同学，特地举办了一场盛大的美食节。作为一个喜欢尝鲜的美食客，小M自然不愿意错过这场盛宴。他很快就尝遍了美食节所有的美食。然而，尝鲜的欲望是难以满足的。尽管所有的菜品都很可口，厨师做菜的速度也很快，小M仍然觉得自己桌上没有已经摆在别人餐桌上的美食是一件无法忍受的事情。于是小M开始研究起了做菜顺序的问题，即安排一个做菜的顺序使得同学们的等待时间最短。小M发现，美食节共有n种不同的菜品。每次点餐，每个同学可以选择其中的一个菜品。总共有m个厨师来制作这些菜品。当所有的同学点餐结束后，菜品的制作任务就会分配给每个厨师。然后每个厨师就会同时开始做菜。厨师们会按照要求的顺序进行制作，并且每次只能制作一人份。此外，小M还发现了另一件有意思的事情: 虽然这m个厨师都会制作全部的n种菜品，但对于同一菜品，不同厨师的制作时间未必相同。他将菜品用1, 2, ..., n依次编号，厨师用1, 2, ..., m依次编号，将第j个厨师制作第i种菜品的时间记为 ti,j 。小M认为：每个同学的等待时间为所有厨师开始做菜起，到自己那份菜品完成为止的时间总长度。换句话说，如果一个同学点的菜是某个厨师做的第k道菜，则他的等待时间就是这个厨师制作前k道菜的时间之和。而总等待时间为所有同学的等待时间之和。现在，小M找到了所有同学的点菜信息: 有 pi 个同学点了第i种菜品（i=1, 2, ..., n）。他想知道的是最小的总等待时间是多少。

输入文件的第1行包含两个正整数n和m，表示菜品的种数和厨师的数量。第2行包含n个正整数，其中第i个数为pi，表示点第i种菜品的人数。接下来有n行，每行包含m个非负整数，这n行中的第i行的第j个数为ti,j，表示第j个厨师制作第i种菜品所需的时间。输入文件中每行相邻的两个数之间均由一个空格隔开，行末均没有多余空格。

输出仅一行包含一个整数，为总等待时间的最小值。

3 2
3 1 1
5 7
3 6
8 9

47

【样例说明】
厨师1先制作1份菜品2，再制作2份菜品1。点这3道菜的3个同学的等待时间分别为3，3+5=8，3+5+5=13。
厨师2先制作1份菜品1，再制作1份菜品3。点这2道菜的2个同学的等待时间分别为7，7+9=16。
总等待时间为3+8+13+7+16=47。
虽然菜品1和菜品3由厨师1制作更快，如果这些菜品都由厨师1制作，总等待时间反而更长。如果按上述的做法，将1份菜品1和1份菜品3调整到厨师2制作，这样厨师2不会闲着，总等待时间更短。
可以证明，没有更优的点餐方案。

【数据规模及约定】
对于100%的数据，n <= 40, m <= 100, p <= 800, ti,j <= 1000（其中p = ∑pi，即点菜同学的总人数）。
每组数据的n、m和p值如下：
测试点编号 n m p
1 n = 5 m = 5 p = 10
2 n = 40 m = 1 p = 400
3 n = 40 m = 2 p = 300
4 n = 40 m = 40 p = 40
5 n = 5 m = 40 p = 100
6 n = 10 m = 50 p = 200
7 n = 20 m = 60 p = 400
8 n = 40 m = 80 p = 600
9 n = 40 m = 100 p = 800
10 n = 40 m = 100 p = 800

弱菜一开始不会建图，各种被虐。

然后总算会建了（一开始想着拆边，完全忽略还能拆点TNT），动态加边，继续被虐 60分 2 Hit

给SPFA+ 上了 bool 是否进队判断 90分 3 Hit

被虐被虐被虐被虐。。。。。

研究N次标程后,把SLF 删了 100分我X

研究结果:SPFA SLF小根 80分 SLF大根 90分没有SLF 100分

丧心病狂，有人故意卡么。。。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
#include<cmath>
#include<cctype>
#include<cassert>
#include<climits>
#include<queue>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define RepD(i,n) for(int i=n;i>=0;i--)
#define MEM(a) memset(a,0,sizeof(a))
#define MEMI(a) memset(a,127,sizeof(a))
#define MEMi(a) memset(a,128,sizeof(a))
#define INF (2139062143)
#define F (1000000009)
#define MAX_N (400+10)
#define MAX_M (100+10)
#define MAX_P (800+10)
#define MAXN (80400+10)
#define MAXM (720000+2600+10)
#define read(x) { \
                while (!isdigit(c=getchar())); \
                x=c-48; \
                while (isdigit(c=getchar())) x=x*10+c-48; \
                }
typedef long long ll;
char c;
int n1,n,s,t,n2,p[MAX_N],tim[MAX_N][MAX_M];
int edge[MAXM],next[MAXM]={0},pre[MAXN]={0},size=1;
ll weight[MAXM],cost[MAXM];
void addedge(int u,int v,int w,int c)
{
   edge[++size]=v;
   weight[size]=w;
   cost[size]=c;
   next[size]=pre[u];
   pre[u]=size;
}
void addedge2(int u,int v,int w,int c){addedge(u,v,w,c),addedge(v,u,0,-c);}
int pre2[MAXN],q[MAXN*8];
int d[MAXN];
int who_cook_num[MAXN][2];
bool b[MAXN]={0};
void SPFA(int s)
{
   memset(d,127,sizeof(d));
   d[s]=0;
   int head=1,tail=1;
   q[1]=s;pre2[s]=0;
   while (head<=tail)
   {
      int now=q[head];b[now]=0;head++;
      Forp(now)
      {
         int &v=edge[p];
         if (weight[p]>0&&edge[p]^s&&d[now]+cost[p]<d[v]) 
         {
            d[v]=d[now]+cost[p];
            pre2[v]=p;
            if (b[v]) continue;b[v]=1;
            /*if (head>1&&d[v]>d[head])
            {
               q[--head]=v;
            }else*/ q[++tail]=v;
         }
      }
                             
   }  
}
ll costflow(int s)
{
   ll totcost=0,totflow=0;
   while(1)
   {
      SPFA(s);if (d[t]==INF) break;
      ll fl=INF,nowcost=0;
      for(int i=t;i^s;i=edge[pre2[i]^1]) fl=min(fl,weight[pre2[i]]),nowcost+=cost[pre2[i]];
      for(int i=t;i^s;i=edge[pre2[i]^1]) weight[pre2[i]]-=fl,weight[pre2[i]^1]+=fl;
      int whos=edge[pre2[t]^1];
      who_cook_num[++n][0]=who_cook_num[whos][0],who_cook_num[n][1]=who_cook_num[whos][1]+1;
      For(i,n1) addedge2(i,n,1,tim[i][who_cook_num[n][0]]*who_cook_num[n][1]);
      addedge2(n,t,1,0);
      totcost+=fl*nowcost;totflow+=fl;      
   }
   return totcost;
}
int main()
{
// freopen("wiki1935.in","r",stdin);
   read(n1) read(n2)
   int S=0;
   For(i,n1)
   {
      read(p[i]);
      S+=p[i];
   }
   n=n1+n2+2,s=n-1,t=n;
   For(i,n1)
   {
      addedge2(s,i,p[i],0);
   }
   For(i,n2)
   {
      addedge2(n1+i,t,1,0);
      who_cook_num[n1+i][0]=i,who_cook_num[n1+i][1]=1;
   }
   
   For(i,n1)
      For(j,n2)
      {
         read(tim[i][j]);
         addedge2(i,j+n1,1,tim[i][j]); 
      }
   
   cout<<costflow(s)<<endl;
   
// while(1);   
   return 0;
}