P1251 餐巾计划问题

题目

题目

思路

啊哈，本题实质上是拆点+MCMF，我们把明天拆成早上和晚上，按照题意建边（不详细描述，见代码），跑MCMF即可。
code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
int n,s,t,x,y,tot=2,head[5001],dep[5001],fy;
bool book[5001];
long long w,ans,in[5001],ans2;
long long mn(long long x,long long y)
{
 return (x>y?y:x);
}
struct f{
 int to,net;
 long long w,fy;
} a[1000011];
int b[1000011];
void add(int x,int y,long long w,long long fy)
{
 a[tot].to=y,a[tot].w=w,a[tot].net=head[x],a[tot].fy=fy,head[x]=tot++;
 return;
}
bool spfa()
{
 memset(dep,0x3f,sizeof(dep));
 memset(book,0,sizeof(book));
 queue<int> u;
 u.push(s);
 book[s]=1,dep[s]=0,in[s]=0x7fffffff;
 while (u.size())
 {
  y=u.front();
  u.pop();
  book[y]=0;
  for (int j=head[y];j;j=a[j].net)
  {
   if (a[j].w!=0&&dep[y]+a[j].fy<dep[a[j].to])
   {
    dep[a[j].to]=dep[y]+a[j].fy;//记录单价和 
    in[a[j].to]=min(in[y],a[j].w);
    b[a[j].to]=j;
    if (!book[a[j].to]) book[a[j].to]=1,u.push(a[j].to);
   }
  }
 }
 return dep[t]!=1061109567;
}
void dfs()
{
 x=t;
 ans+=in[t];
 ans2+=dep[t]*in[t];
 while (x!=s)
 {
  a[b[x]].w-=in[t],a[b[x]^1].w+=in[t],x=a[b[x]^1].to; 
 }//清空增广路 
 return;
}
int main()
{
 scanf("%d",&n);
 s=4999,t=5000;
 for (int i=1;i<=n;i++)
 {
  scanf("%d",&w);
  add(s,i+n,w,0);
  add(i+n,s,0,0);
  add(i,t,w,0);
  add(t,i,0,0);
  if (i+1<=n)
  {
   add(i+n,i+n+1,1e30,0);
   add(i+n+1,i+n,0,0);
  }
 }
 int p,m,f,n2,s1;
 scanf("%d%d%d%d%d",&p,&m,&f,&n2,&s1);
 for (int i=1;i<=n;i++)
 {
  if (i+m<=n)
  {
   add(i+n,i+m,1e30,f);
   add(i+m,i+n,0,-f);
  }
  if (i+n2<=n)
  {
   add(i+n,i+n2,1e30,s1);
   add(i+n2,i+n,0,-s1);
  }
  add(s,i,1e30,p);
  add(i,s,0,-p);
 }
 while (spfa())
 {
  dfs();
 }
 cout<<ans2;
    return 0;
}