拯救莫莉斯(from Maple)

　　　　　　　　　　　　拯救莫莉斯

　　　　　　　　　　　　　　　　时间限制: 1 Sec  内存限制: 256 MB

问题描述

莫莉斯·乔是圣域里一个叱咤风云的人物，他凭借着自身超强的经济头脑，牢牢控制了圣域的石油市场。

圣域的地图可以看成是一个n*m的矩阵。每个整数坐标点(x , y)表示一座城市（1<=x<= n, 1<=y<=m）。两座城市间相邻的定义为：对于城市(Ax, Ay)和城市(Bx, By)，满足(Ax - Bx)2 + (Ay - By)2 = 1。

由于圣域的石油贸易总量很大，莫莉斯意识到不能让每笔石油订购单都从同一个油库里发货。为了提高效率，莫莉斯·乔决定在其中一些城市里建造油库，最终使得每一个城市X都满足下列条件之一：

1.该城市X内建有油库，

2.某城市Y内建有油库，且城市X与城市Y相邻。

与地球类似，圣域里不同城市间的地价可能也会有所不同，所以莫莉斯想让完成目标的总花费尽可能少。如果存在多组方案，为了方便管理，莫莉斯会选择建造较少的油库个数。

输入格式

第一行两个正整数n,m ( n * m <= 50 且m<=n)，表示矩阵的大小。

接下来一个n行m列的矩阵F，Fi, j表示在城市(i,j)建造油库的代价。

输出格式

输出两个数，建造方案的油库个数和方案的总代价。

 

输入样例：	输出样例：
3 3 6 5 4 1 2 3 7 8 9	3 6

数据范围

对于30%数据满足 n * m <= 25;

对于100%数据满足n * m <= 50; 0 <= Fi, j <= 100000

因为是考试题，一定要详细点：：

　　首先读到题时要认真审题，是n*m<=50!!当你真真正正的理解题了后，忽然发现，这就是个状压啊~

设f[i][j][k]数组存储的是最小的花费，第一维代表前i行，第二维代表上一行的状态，第三维代表这一行的状态。

预处理出cost[i][j]数组代表在第i行时选取状态j时的花费；

所以很容易得出状态转移方程：f[i][k][u]=min(f[i][j][k]+cost[i][u]);

经过严密的调查发现因为每一个城市至少要有一个油库与之相邻，所以这个状态能够转移过来的前提是(i | j | u | k<<1| k>>1)&(1<<m)-1==(1<<m)-1

k代表当前这一行。这里一定要注意优先级！！

因为转移到第n行时有两维状态不确定，所以我们转移到第n+1行；

在转移f数组时同时维护一个g数组代表最小个数就好了！！

#include<cstdio>
#include<iostream>
#include<cstring>
#define maxn 55
#define edge 8
using namespace std;
inline int read()
{   char c=getchar();int x=0,y=1;
    while(c<'0'||c>'9'){if(c=='-') y=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*y;
}
int ans2=0x7fffffff,n,m,w[maxn][maxn],f[maxn][1<<edge][1<<edge],g[maxn][1<<edge][1<<edge],co[maxn][1<<edge],ans=0x7fffffff,edg,sum[1<<edge];
inline int getn(int x){int tmp=0;while(x){if(x&1) tmp++;x>>=1;}return tmp;}
void dp()
{   for(int i=1;i<=n;i++)
        for(int j=0;j<=edg;j++)
            for(int k=0;k<=edg;k++)
            {   for(int l=0;l<=edg;l++)
                {   if(((l|j|k|(k<<1)|(k>>1))&edg)==edg)
                        if(f[i+1][k][l]>f[i][j][k]+co[i+1][l])
                        {   f[i+1][k][l]=f[i][j][k]+co[i+1][l];
                            g[i+1][k][l]=g[i][j][k]+sum[l];
                        }
                        else if(f[i+1][k][l]==f[i][j][k]+co[i+1][l])
                            g[i+1][k][l]=min(g[i+1][k][l],g[i][j][k]+sum[l]);
                }
            }
}
void init()
{   for(int i=1;i<=n;i++)
        for(int j=0;j<=edg;j++)
            for(int k=1;k<=m;k++)
                if(j&(1<<k-1)) co[i][j]+=w[i][k];
    for(int i=0;i<=edg;i++)
        sum[i]=getn(i);
}
int main()
{   //freopen("proj.in","r",stdin);
    //freopen("proj.out","w",stdout);
    n=read();m=read();
    memset(f,0x7f,sizeof(f));
    edg=(1<<m)-1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            w[i][j]=read();
    init();
    for(int i=0;i<=edg;i++) f[1][0][i]=co[1][i],g[1][0][i]=sum[i];
    dp();
    for(int i=0;i<=edg;i++)
        if(ans>f[n+1][i][0])
            ans=f[n+1][i][0],ans2=g[n+1][i][0];
        else if(ans==f[n+1][i][0]) ans2=min(ans2,g[n+1][i][0]);
    printf("%d %d",ans2,ans);
    return 0;
}
 
/**************************************************************
    Problem: 10425
    User: Anonymity
    Language: C++
    Result: 正确
    Time:72 ms
    Memory:29772 kb
****************************************************************/