ZOJ 3471 Most Powerful (状压DP)
Recently, researchers on Mars have discovered N powerful atoms. All of them are different. These atoms have some properties. When two of these atoms collide, one of them disappears and a lot of power is produced. Researchers know the way every two atoms perform when collided and the power every two atoms can produce.
You are to write a program to make it most powerful, which means that the sum of power produced during all the collides is maximal.
Input
There are multiple cases. The first line of each case has an integer N (2 <= N <= 10), which means there are N atoms: A1, A2, ... , AN. Then N lines follow. There are N integers in each line. The j-th integer on the i-th line is the power produced when Ai and Aj collide with Aj gone. All integers are positive and not larger than 10000.
The last case is followed by a 0 in one line.
There will be no more than 500 cases including no more than 50 large cases that N is 10.
Output
Output the maximal power these N atoms can produce in a line for each case.
Sample Input
2
0 4
1 0
3
0 20 1
12 0 1
1 10 0
0
Sample Output
4
22
Author: GAO, Yuan
Contest: ZOJ Monthly, February 2011
题目的意思是研究人员发现有N个原子两两组合会发生反应产生能量,且其中一个会被和谐掉,给出一个矩阵map,其中map[i][j]表示原子i与j发生反应且原子j被和谐掉所释放的能量,问给定的n个原子反应最多能产生多少能量。
显然直接枚举是不可行的,因为复杂度是O(n!),虽然10!在10^7左右,但是有50个测试数据的n=10,所以估计会超时的。
假设一个数,第i位表示第i个原子是否被灭掉,如果被灭掉则为1,没被灭掉为0,那么所有状态都可以用2^n范围内的数来表示。假设对于其中一个数(n=4)
1 0 1 0,(也就是十进制的10)则表示第二个原子和第四个原子都被灭掉,而第一、三个原子幸存,那么,假设知道这个状态下释放能量的最大值,那么我们可以拓展出1 1 1 0 和1 0 1 1 这两个状态通过1 0 1 0达到的值,那么对于每个状态进行这种拓展,那么就可以求出每个状态下释放能量的最大值,所以最后的最大值就求出来了。
#include<iostream> #include<cstdio> #include<cstring> using namespace std; int n,map[15][15]; int dp[1200]; int main(){ //freopen("input.txt","r",stdin); while(~scanf("%d",&n) && n){ for(int i=0;i<n;i++) for(int j=0;j<n;j++) scanf("%d",&map[i][j]); for(int i=0;i<(1<<n);i++) dp[i]=0; for(int i=0;i<(1<<n);i++) //对于每个状态进行拓展 for(int j=0;j<n;j++){ //枚举每一位,也就是说枚举每个原子将其和谐掉 if((i&(1<<j))!=0) //如果在状态i下该原子已经被和谐了,那么跳过看下一个原子 continue; for(int k=0;k<n;k++) //否则找一个没有被和谐的原子和谐掉原子 if((i&(1<<k))==0 && j!=k) dp[i|(1<<j)]=max(dp[i|(1<<j)],dp[i]+map[k][j]); } int ans=0; for(int i=0;i<(1<<n);i++) if(ans<dp[i]) ans=dp[i]; printf("%d\n",ans); } return 0; }
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· Spring AI + Ollama 实现 deepseek-r1 的API服务和调用
· 数据库服务器 SQL Server 版本升级公告
· C#/.NET/.NET Core技术前沿周刊 | 第 23 期(2025年1.20-1.26)
· 程序员常用高效实用工具推荐,办公效率提升利器!