HDU 4049

来源：互联网发布：mac运行32位软件编辑：程序博客网时间：2024/06/06 01:59

好题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;}

0 0