状压DP 拯救莫莉斯
来源:互联网 发布:软件著作权申请加急 编辑:程序博客网 时间:2024/05/16 07:59
问题 C: 拯救莫莉斯
时间限制: 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
有点类似那个奶牛翻瓦片的问题传送门 ,因为m<n,所以m最大得7,状压每一行。某一行的状态由自己,上一行,下一行决定,下一行不可知,所以要保存上一行与这一行的选择状态。f[i][j][k]都表示啥我就不说了。OWO
f[i+1][j][k]=min(f[i][h][j]+val[i+1][k]); val表示第i+1行选k状态时耗费。
多维护一个g数组,就有了最小个数了。
那么啥时候能推到下一行呢?
f[i+1][j][k]和f[i][h][j] 则(h|k|k|k<<1|k>>1)&(2^m-1)==2^m-1就说明i这一行已经全部覆盖了。
最后一个。。一个问题(很严重)初始化和最后输出,我太菜初始化想了半天,旁边ltr大佬一遍想出。。%%%
初始化 f[1][0][i]=balabala()
输出 max(f[n+1][i][0]);
#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<ctime>#include<algorithm>using namespace std;int read(){int sum=0,f=1;char x=getchar();while(x<'0'||x>'9'){if(x=='-')f=-1;x=getchar();}while(x>='0'&&x<='9'){sum=sum*10+x-'0';x=getchar();}return sum*f;}int n,m,a[55][55],v[55][(1<<7)+2],ans=10000000,cnt=100;int f[55][(1<<7)+5][(1<<7)+5],g[55][(1<<7)+2][(1<<7)+2];int get(int x){int s=0;for(int i=0;i<m;i++) if(x&(1<<i))s++;return s;}int check(int x,int y,int z){return ((y|z|x|(x<<1)|(x>>1))&((1<<m)-1))==((1<<m)-1);}int main(){ // freopen("proj.in","r",stdin); // freopen("proj.out","w",stdout); n=read();m=read(); for(int i=1;i<=n;i++) { for(int j=1;j<=m;j++) a[i][j]=read(); for(int j=0;j<(1<<m);j++){for(int k=0;k<m;k++) if((1<<k)&j) v[i][j]+=a[i][k+1];v[i][0]=0;} } memset(f,12,sizeof(f)); memset(g,20,sizeof(g)); //for(int i=0;i<(1<<m);i++) for(int j=0;j<(1<<m);j++) f[1][0][j]=v[1][j],g[1][0][j]=get(j); for(int i=1;i<=n;i++) for(int j=0;j<(1<<m);j++) for(int k=0;k<(1<<m);k++) for(int h=0;h<(1<<m);h++) if(check(j,k,h)) { int d=get(k); if(f[i+1][j][k]>f[i][h][j]+v[i+1][k]) f[i+1][j][k]=f[i][h][j]+v[i+1][k],g[i+1][j][k]=g[i][h][j]+d; else if(f[i+1][j][k]==f[i][h][j]+v[i+1][k]) g[i+1][j][k]=min(g[i+1][j][k],g[i][h][j]+d); } for(int i=0;i<(1<<m);i++) if(ans>f[n+1][i][0]) ans=f[n+1][i][0],cnt=g[n+1][i][0]; else if(ans==f[n+1][i][0]) cnt=min(g[n+1][i][0],cnt); cout<<cnt<<" "<<ans;}
- 状压DP 拯救莫莉斯
- 拯救莫莉斯 状压dp
- 拯救莫莉斯
- 【bzoj3174】【TJOI2013】【拯救小矮人】【贪心+dp】
- 3643. 【GDOI2014】拯救莫莉斯
- 拯救莫莉斯(from Maple)
- 拯救
- 双向搜索+DP——拯救公主(二)
- BZOJ3174 Tjoi2013 拯救小矮人(贪心+DP)
- 【BZOJ3174】【codevs2544&&2075】拯救小矮人,DP+贪心
- [noip测试]拯救紫萱学姐(kmp+树形dp)
- 【NOIP 模拟题】[T2]拯救紫萱学姐(kmp+树形dp)
- BZOJ 3174: [Tjoi2013]拯救小矮人 贪心 dp
- 拯救拯救自己
- 拯救程序员
- 拯救wap
- 拯救。。20070422
- 拯救情感
- MySQL存储过程写法总结
- JavaScript作用域链和原型链
- 2018届秋招提前批阿里在线编程
- HTML5 总结
- Web 项目中 JSP 和 HTML 的取舍和资源引用时 "\" 根目录的寻找
- 状压DP 拯救莫莉斯
- Good Luck in CET-4 Everybody!(博弈)
- 创建型-单例模式(singleton)
- CountDownLatch监听初始化
- python str.format与制表符\t关于中文对齐的细节问题
- 【JVM】JVM执行引擎
- uva 10082 WERTYU
- 欢迎使用CSDN-markdown编辑器
- Next ESA SAR Toolbox (NEST)