HDU 4049

好题mark

1. 图论方法（最大权闭合子图）

#include <cstdio>
#include <cstring>
#define min(a,b) ((a)<(b)?(a):(b))
 
const int N=1110;
const int E=5000;
const int oo=1000000000;
 
int node,src,dest,ne;
int head[N],work[N],Q[N],dist[N];
int pnt[E],nxt[E],flow[E];
int val[N][N],bonus[N][N],cost[N];
 
inline void init(int _node,int _src,int _dest)
{
    node=_node;
    src=_src;
    dest=_dest;
    for(int i=ne=0;i<node;i++)   head[i]=-1;
}
inline void add(int u,int v,int c1,int c2)
{
    pnt[ne]=v,flow[ne]=c1,nxt[ne]=head[u],head[u]=ne++;
    pnt[ne]=u,flow[ne]=c2,nxt[ne]=head[v],head[v]=ne++;
}
bool dinic_bfs(void)
{
    int i,u,v,l,r=0;
    for(i=0;i<node;i++)  dist[i]=-1;
    dist[Q[r++]=src]=0;
    for(l=0;l<r;l++) for(i=head[u=Q[l]];~i;i=nxt[i])
        if(flow[i] && dist[v=pnt[i]]<0)
        {
            dist[Q[r++]=v]=dist[u]+1;
            if(v==dest) return 1;
        }
    return 0;
}
int dinic_dfs(int u,int exp)
{
    if(u==dest) return exp;
    for(int &i=work[u],v,tmp;~i;i=nxt[i])
        if(flow[i] && dist[v=pnt[i]]==dist[u]+1 && (tmp=dinic_dfs(v,min(flow[i],exp)))>0)
        {
            flow[i]-=tmp;
            flow[i^1]+=tmp;
            return tmp;
        }
    return 0;
}
int dinic_flow(void)
{
    int i,res=0,delta;
    while(dinic_bfs())
    {
        for(i=0;i<node;i++)  work[i]=head[i];
        while(delta=dinic_dfs(src,oo))  res+=delta;
    }
    return res;
}
int main(void)
{
    int n,m;
    while(scanf("%d%d",&n,&m)==2 && (n||m))
    {
        init(n*m+n*n*m+2,0,n*m+n*n*m+1);
        for(int i=1;i<=m;i++)    scanf("%d",cost+i);
        for(int i=1;i<=n;i++)    for(int j=1;j<=m;j++)
        {
            scanf("%d",&val[i][j]);
            val[i][j]-=cost[j];
        }
        for(int i=1;i<=n;i++)    for(int j=1;j<=n;j++)    scanf("%d",&bonus[i][j]);
        int sum=0;
        for(int i=1;i<=n;i++)    for(int j=1;j<=m;j++)
        {
            if(val[i][j]>0)
            {
                sum+=val[i][j];
                add(src,(i-1)*m+j,val[i][j],0);
            }
            else    add((i-1)*m+j,dest,-val[i][j],0);
            if(j>1)  add((i-1)*m+j,(i-1)*m+j-1,oo,0);
        }
        for(int i=1;i<=n;i++)    for(int j=i+1;j<=n;j++)  for(int k=1;k<=m;k++)
        {
            sum+=bonus[i][j];
            add(src,n*m+(((i-1)*n+j)-1)*m+k,bonus[i][j],0);
            add(n*m+(((i-1)*n+j)-1)*m+k,(i-1)*m+k,oo,0);
            add(n*m+(((i-1)*n+j)-1)*m+k,(j-1)*m+k,oo,0);
        }
        int res=sum-dinic_flow();
        if(res) printf("%d\n",res);
        else    puts("STAY HOME");
    }
    return 0;
}

2. DP方法（状态压缩）

#include<string.h>
#include<iostream>
#include <stdio.h>
#include <memory.h>
using namespace std;
//typedef long long LL;
typedef __int64 LL;
const int maxn=1<<11;
LL A[11][11],B[11],g[11][11],dp[12][maxn],w[12][maxn],ww[maxn][maxn],m,n;

LL Cal()
{
    LL i1,j1,k1,ii,kk,i,j,k,a[15];

    for(i=0;i<m;i++)
        for(j=0;j<1<<n;j++)
            dp[i][j]=-100000000;
    memset(w,0,sizeof(w));



    for(j=0;j<1<<n;j++)
    {
        i=j,ii=0,k=0;
        while(i)
        {
            if(i&1)
                a[ii++]=k;
            k++;
            i=i>>1;
        }
        for(kk=0;kk<m;kk++)
        {
            for(i=0;i<ii;i++)
                w[kk][j]+=A[a[i]][kk]-B[kk];
            for(i=0;i<ii;i++)
                for(k=i+1;k<ii;k++)
                    w[kk][j]+=g[a[i]][a[k]];
        }
    }//这里计算在ww景点，状态j游玩此景点的收益



    for(j=0;j<1<<n;j++)
        dp[0][j]=w[0][j];//初始化dp，



    for(i=0;i<1<<n;i++)
    {
        for(k=1,j=0;j<=i;j++)
            if((j|i)==i)
                ww[i][k++]=j;
        ww[i][0]=k;

    }//这里for循环做的是寻找如果状态i游玩了上一个景点，那么游玩下一个景点的状态可能有ww[i][0]个，分别是从ww[i][1]到ww[i][ww[i][0]]



    for(i1=1;i1<m;i1++)//在i1点，在i1-1点剩的人的二进制表示为k1
    {
        for(j1=0;j1<1<<n;j1++)
        {
            for(k1=1;k1<ww[j1][0];k1++)
            {
                ii=ww[j1][k1];
                if(dp[i1][ii]<dp[i1-1][j1]+w[i1][ii])
                    dp[i1][ii]=dp[i1-1][j1]+w[i1][ii];
            }
        }
    }//dp[i][j]表示状态j在i游玩，则总的最大收益是这个

    for(kk=0,i=0;i<1<<n;i++)
        if(dp[m-1][i]>kk) kk=dp[m-1][i];
    return kk;
}
void Init()
{
    LL i,j;
    for(i=0;i<m;i++)//每个人进所需的钱
        scanf("%I64d",&B[i]);
    for(i=0;i<n;i++)//i人到j地方的兴趣
        for(j=0;j<m;j++)
            scanf("%I64d",&A[i][j]);
    for(i=0;i<n;i++)//i与j同时进再加的兴趣
        for(j=0;j<n;j++)
            scanf("%I64d",&g[i][j]);
}

int main()
{
    LL ii;
    while(scanf("%I64d%I64d",&n,&m)!=EOF)
    {
        if(n==0&&m==0) break;
        Init();
        ii=Cal();
        if(ii==0) printf("STAY HOME\n");
        else printf("%I64d\n",ii);
    }
    return 0;
}