『HDU 5855』Less Time, More profit

题目链接:acm.hdu.edu.cn/showproblem.php?pid=5855

个人感想:
当时看错题意了,尼玛一直往贪心想,然后就GG了,贪心也可以过,只是因为数据太水了… 在这道题中并没有贪心策略可以选择呢…

然后出题解之后,我内心很崩溃,签到题居然是网络流,简直疯了…
不过确实是个水题,完完全全是个模板题,套上去就可以了…但首先,你得明白什么是最大权闭合图.
这个就百度一下吧.
然后分享个连接:
http://wenku.baidu.com/view/d9c9b9220722192e4536f6e1.html?from=search&qq-pf-to=pcqq.discussion
这个就是要搞明白最大流最小割问题,从而求最大收益…我尼玛我以前只是看过一下…我没当回事,我感觉并不会出这种题啊…真是acmer是上知天文,下知地理…,然后我今天恶补了概念,并把太空飞行计划问题刷了一遍.这完完全全是套这个模板的.只要看懂了,这很好做了…我不需要过渡解释了,这篇博文通俗易懂,没有一大堆公式,最好能保存下来以后用.

因为变量太多,所以要把时间t来进行二分来贪心就好了…

我贴两份代码,这都是自己写的,一个是太空飞行计划问题代码,和现在题目的AC代码.好好看一遍,这学习到很多东西.网络流模板来自于**《挑战程序与设计竞赛》**.

分析： 二分+网络流

太空飞行计划问题代码:

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e6;
int n,m;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }

 int label[maxn];
 struct edge
 {
     int to;
     int cap;
     int rev;
 };
 vector<edge> G[maxn];
 int level[maxn];
 int iter[maxn];

 void add_edge(int from,int to,int cap)
 {
    // cout<<from<<"-->"<<to<<" "<<cap<<endl;
     G[from].push_back((edge){to,cap,G[to].size()});
     G[to].push_back((edge){from,0,G[from].size()-1});
 }

 void bfs(int s)
 {
     memset(level,-1,sizeof(level));
     queue<int> que;
     level[s]=0;
     que.push(s);
     while(!que.empty())
     {
         int v=que.front();
         que.pop();
         for(int i=0;i<G[v].size();i++)
         {
             edge &e=G[v][i];
             if(e.cap>0 &&level[e.to]<0)
             {
                 level[e.to]=level[v]+1;
                 que.push(e.to);
             }
         }
     }
 }
 int dfs(int v,int t,int f)
 {
     if(v==t) return f;
     for(int &i=iter[v];i<G[v].size();i++)
     {
         edge &e=G[v][i];
         if(e.cap>0 && level[v]<level[e.to])
         {
             int d=dfs(e.to,t,min(f,e.cap));
             if(d>0)
             {
                 e.cap-=d;
                 G[e.to][e.rev].cap +=d;
                 return d;
             }
         }
     }
     return 0;
 }
int max_flow(int s,int t)
{
    int flow=0;
    for (;;)
    {
        bfs(s);
        if(level[t]<0) return flow;
        memset(iter,0,sizeof(iter));
        int f;
        while((f=dfs(s,t,INF))>0)
        {
            flow+=f;
        }
    }
}
void dfs1(int x)
{
   // cout<<x<<endl;
    label[x]=1;
    for(int j=0;j<G[x].size();j++)
    {
        edge &e=G[x][j];
        //cout<<e.to<<endl;
        if(e.cap>0&&!label[e.to]) dfs1(e.to);
    }
}
int main()
{

    freopen("coco.txt","r",stdin);
    freopen("lala.txt","w",stdout);


     int s,t;
     int val,to;
     scanf("%d%d",&n,&m);
     s=n+m+1;
     t=n+m+2;
     int total=0;
     for(int i=1;i<=n;i++)
     {

         scanf("%d",&val);
         add_edge(s,i,val);
         total+=val;

         char c;
         while(c=getchar())
         {
             if(c=='\n') break;
             else if(isdigit(c))
             {
                 ungetc(c,stdin);
                 scanf("%d",&to);
                 add_edge(i,n+to,INF);
             }
         }
     }
     for(int i=1;i<=m;i++)
     {
        scanf("%d",&val);
        add_edge(n+i,t,val);

     }
     int ans=max_flow(s,t);
     memset(label,0,sizeof(label));

     dfs1(s);
     //cout<<ans<<endl;
     for(int i=1;i<=n;i++)if( (label[i]^label[s]) ==0) cout<<i<<' ';
     cout<<endl;
     for(int i=1;i<=m;i++)if( label[n+i]^label[t]) cout<<i<<' ';
     cout<<endl;

     cout<<total-ans<<endl;
    return 0;
}

AC代码（网络流用了dinic算法,我只是担心EK很慢…）:

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=405;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }


int T;
int n,m,k;

int N,s,t;
struct edge
{
    int to;
    int cap;
    int rev;
};

vector<edge> G[maxn];
int level[maxn];
int iter[maxn];
bool vis[maxn];

void add_edge(int from,int to,int cap)
{
    G[from].push_back((edge){to,cap,G[to].size()});
    G[to].push_back((edge){from,0,G[from].size()-1});
}
void bfs(int s)
{
    memset(level,-1,sizeof(level));
    queue<int> que;
    level[s]=0;
    que.push(s);
    while(!que.empty())
    {
        int v=que.front();
        que.pop();
        for(int i=0;i<G[v].size();i++)
        {
            edge &e=G[v][i];
            if(e.cap >0 && level[e.to]<0)
            {
                level[e.to]=level[v]+1;
                que.push(e.to);
            }
        }
    }
}
int dfs(int v,int t,int f)
{
    if(v == t) return f;
    for(int &i=iter[v];i<G[v].size();i++)
    {
        edge &e=G[v][i];
        if(e.cap>0 && level[v]< level[e.to])
        {
            int d=dfs(e.to,t,min(f,e.cap));
            if(d>0)
            {
                e.cap-=d;
                G[e.to][e.rev].cap += d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t)
{
    int flow=0;
    for(;;)
    {
        bfs(s);
        if(level[t]<0) return flow;
        memset(iter,0,sizeof(iter));
        int f;
        while((f= dfs(s,t,INF))>0)
        {
            flow+=f;
        }
    }
}





class Indestry
{
public:
    int pay;
    int t;
};
Indestry I[maxn];
class Shop
{
public:
    int id;
    int profit;
    int MAX;
    int Iid[205];
};
Shop S[maxn];






int judge(int ti)
{
    int total=0;
    memset(vis,0,sizeof(vis));

    for(int i=0;i<=N;i++) G[i].clear();

    for(int i=1;i<=m;i++)
    {
       if(S[i].MAX<=ti)
       {
           total+=S[i].profit;
           add_edge( S[i].id ,t,S[i].profit);
           for(int j=1;j<=S[i].Iid[0];j++)
           {
                int son=S[i].Iid[j];
                if(!vis[son]){add_edge(s,son,I[son].pay);}
                vis[son]=true;
                add_edge(son,S[i].id,INF);
           }
       }
    }


    int ans=max_flow(s,t);

    return total-ans;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("coco.txt","r",stdin);
    freopen("lala.txt","w",stdout);
#endif
    T=Scan();
    for(int cas=1;cas<=T;cas++)
    {

        int l=0;
        int r=0;

        int val=-INF;//枚举用的答案.

        n=Scan();
        m=Scan();
        k=Scan();

        printf("Case #%d: ",cas);
        N=n+m+2;

        s=N-1;
        t=N;

        for(int i=1;i<=n;i++)//工厂
        {
          I[i].pay=Scan();
          I[i].t=Scan();

          r=max(I[i].t,r);//枚举2分用的变量.
        }
        r+=1;
        for(int i=1;i<=m;i++)//商店,顺便把商店所需要的工厂建立起来.
        {
            int mx=0;

            S[i].profit=Scan();
            S[i].id=n+i;

            int x;

            S[i].Iid[0]=Scan();
            for(int j=1;j<=S[i].Iid[0];j++)
            {
                x=Scan();
                S[i].Iid[j]=x;
                mx=max(mx,I[x].t);
            }
            S[i].MAX=mx;
        }


        while(r-l>1)
        {
            int mid=(r+l)/2;
            int tval=judge(mid);
            if(tval>=k)
            {
                val=tval;
                r=mid;
            }
            else l=mid;
        }


        if(val==-INF) printf("impossible\n");
        else printf("%d %d\n",r,val);



    }


    return 0;
}