hdu 4309 状压+最大流

传送门

题目大意：有n个城市，每个城市有一个人数，城市与城市之间有三种通道(均为单向)：隧道，公路和桥。隧道和公路都可以通过无限制的人，并且隧道可以避难ai人。桥没有修理的时候只能通过一个人次，修理之后就和公路没有差别了，修理有费用。特别注意桥不超过12条。求通过最大人数，最少修理话费。

最多只有12个。。。。直接枚举可能就过了--4000+种可能情况而已。每次输入的时候记录有多少个桥需要修 cnt++ 状压枚举的时候一定是枚举修桥之后的情况

枚举方法(状压 1<<cnt

每个状态求一次流量，然后更新答案。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>

using namespace std;

const int maxn = 60007;
const int maxm = 65000;
const int oo = 99999999;

int idx;
int cur[maxn], pre[maxn];
int dis[maxn], gap[maxn];
int aug[maxn], head[maxn];

struct Node
{
    int u, v, w;
    int next;
}edge[maxm];

void addEdge(int u, int v, int w)
{
    edge[idx].u = u;
    edge[idx].v = v;
    edge[idx].w = w;
    edge[idx].next = head[u];
    head[u] = idx++;

    edge[idx].u = v;
    edge[idx].v = u;
    edge[idx].w = 0;
    edge[idx].next = head[v];
    head[v] = idx++;
}


int SAP(int s, int e, int n)
{
    int max_flow = 0, v, u = s;
    int id, mindis;
    aug[s] = oo;
    pre[s] = -1;
    memset(dis, 0, sizeof(dis));
    memset(gap, 0, sizeof(gap));
    gap[0] = n; // 我觉得这一句要不要都行，因为dis[e]始终为0
    for (int i = 0; i <= n; ++i)
    {   // 初始化当前弧为第一条弧
        cur[i] = head[i];
    }

    while (dis[s] < n)
    {
        bool flag = false;
        if (u == e)
        {
            max_flow += aug[e];
            for (v = pre[e]; v != -1; v = pre[v]) // 路径回溯更新残留网络
            {
                id = cur[v];
                edge[id].w -= aug[e];
                edge[id^1].w += aug[e];
                aug[v] -= aug[e]; // 修改可增广量，以后会用到
                if (edge[id].w == 0) u = v; // 不回退到源点，仅回退到容量为0的弧的弧尾
            }
        }
        for (id = cur[u]; id != -1; id = edge[id].next)
        {   // 从当前弧开始查找允许弧
            v = edge[id].v;
            if (edge[id].w > 0 && dis[u] == dis[v] + 1) // 找到允许弧
            {
                flag = true;
                pre[v] = u;
                cur[u] = id;
                aug[v] = min(aug[u], edge[id].w);
                u = v;
                break;
            }
        }
        if (flag == false)
        {
            if (--gap[dis[u]] == 0) break; /* gap优化，层次树出现断层则结束算法 */
            mindis = n;
            cur[u] = head[u];
            for (id = head[u]; id != -1; id = edge[id].next)
            {
                v = edge[id].v;
                if (edge[id].w > 0 && dis[v] < mindis)
                {
                    mindis = dis[v];
                    cur[u] = id; // 修改标号的同时修改当前弧
                }
            }
            dis[u] = mindis + 1;
            gap[dis[u]]++;
            if (u != s) u = pre[u]; // 回溯继续寻找允许弧
        }
    }
    return max_flow;
}
int val[1010];
struct node
{
    int u,v,w,p;
}Q[1010];
//int map[300][300];
int main()
{
    int t, n, m, pi, si, ei;
    int Max, sum, source, sink, vn;
    //scanf("%d", &t);
    t=1;
    for (int cas = 1; cas <= t; ++cas)
    {
        //scanf("%d%d",&n,&m);
        while(scanf("%d%d",&n,&m)!=EOF)
        {
            int flag=0;
            for(int i=1;i<=n;i++)
                scanf("%d",&val[i]);
          //  getchar();
            int cnt=0;
            for(int i=1;i<=m;i++)
            {
                scanf("%d%d%d%d",&Q[i].u,&Q[i].v,&Q[i].w,&Q[i].p);
                if(Q[i].p>0)
                    cnt++;
                if(Q[i].p<0&&Q[i].w!=0)
                {
                    flag=1;
                   // cout<<"*****"<<endl;
                }
               // cout<<Q[i].p<<" ";
            }
            if(!flag)
            {
                puts("Poor Heaven Empire");
                continue;
            }
            int mm=(1<<cnt);

            int maxnum=-oo,mincost=0;
            for(int k=0;k<mm;k++)
            {
                idx = 0;
                source = 0;sink = n + 1; vn = sink + 1;
                memset(head, -1, sizeof(head));
                sum = 0;  Max = 0;
           /* for (int i = 1; i <= n; ++i)
            {
                scanf("%d %d %d", &pi, &si, &ei);
                sum += pi;
                Max = max(Max, ei);
                addEdge(source, i, pi);
                for (int j = si; j <= ei; ++j)
                {
                    addEdge(i, n + j, 1);
                }
            }  */
                //cout<<endl;
                int j=0,cost=0;
                for(int i=1;i<=n;i++)
                    addEdge(0,i,val[i]);
                for(int i=1;i<=m;i++)
                {
                    //cout<<Q[i].u<<" "<<Q[i].v<<" "<<Q[i].w<<" "<<Q[i].p<<endl;
                    if(Q[i].p<0)
                    {
                        addEdge(Q[i].u,Q[i].v,oo);
                        addEdge(Q[i].u,sink,Q[i].w);
                    }
                    else if(Q[i].p==0)
                        addEdge(Q[i].u,Q[i].v,oo);
                    else
                    {
                        //cout<<Q[i].p<<"   ****"<<endl;
                        if(k&(1<<j))
                        {
                            addEdge(Q[i].u,Q[i].v,oo);
                            cost+=Q[i].w;
                           // cout<<cost<<endl;
                        }
                        else addEdge(Q[i].u,Q[i].v,1);
                        j++;
                    }
                }
                int ans=SAP(0,sink,vn);
                if(ans>maxnum||(ans==maxnum&&cost<mincost))
                {
                    maxnum=ans;
                    mincost=cost;
                }
            }
            if(maxnum>0)
                printf("%d %d\n",maxnum,mincost);
            else printf("Poor Heaven Empire\n");
            /*for (int i = 1; i <= Max; ++i)
            {
                addEdge(n + i, sink, m);
            }
            if (SAP(source, sink, vn) == sum)
                printf("Case %d: Yes\n\n", cas);
            else printf("Case %d: No\n\n", cas);  */
        }
    }
    return 0;
}