D10：洪水&&[COCI2006-2007#1] Slikar+邦德I

1235. 洪水 (Standard IO)

Time Limits: 1000 ms  Memory Limits: 65536 KB

Description

一天， 一个画家在森林里写生，突然爆发了山洪，他需要尽快返回住所中，那里是安
全的。
森林的地图由R行C列组成，空白区域用点“.”表示，洪水的区域用“*”表示，而
岩石用“X”表示，另画家的住所用“D”表示，画家用“S”表示。
有以下几点需要说明：
1、 每一分钟画家能向四个方向移动一格（上、下、左、右）
2、 每一分钟洪水能蔓延到四个方向的相邻格子（空白区域）
3、 洪水和画家都不能通过岩石区域
4、 画家不能通过洪水区域（同时也不行，即画家不能移到某个格子，该格子在画家达到的同时被洪水蔓延到了，这也是不允许的）
5、 洪水蔓不到画家的住所。
给你森林的地图，编写程序输出最少需要花费多长时间才能从开始的位置赶回家中。

Input

输入第一行包含两个整数R和C(R,C<=50)。
接下来R行每行包含C个字符(“.”、“*”、“X”、“D”或“S”)。地图保证只有一个“D”和一个“S”。

Output

输出画家最快安全到达住所所需的时间，如果画家不可能安全回家则输出“KAKTUS”。

Sample Input

输入1：
3 3 
D.* 
... 
.S. 

输入2：
3 3 
D.* 
...
..S

输入3：
3 6 
D...*. 
.X.X.. 
....S. 

Sample Output

输出1：
3

输出2：
KAKTUS 

输出3：
6

 思路

 

代码

 

1236. 邦德I【推荐】 (Standard IO)

Time Limits: 2000 ms  Memory Limits: 65536 KB

Description

每个人都知道詹姆斯邦德，著名的007，但很少有人知道很多任务都不是他亲自完成的，而是由他的堂弟们吉米邦德完成（他有很多堂弟），詹姆斯已经厌倦了把一个个任务分配给一个个吉米，他向你求助。
每个月，詹姆斯都会收到一些任务，根据他以前执行任务的经验，他计算出了每个吉米完成每个任务的成功率，要求每个任务必须分配给不同的人去完成，每个人只能完成一个任务。
请你编写程序找到一个分配方案使得所有任务都成功完成的概率。

Input

输入第一行包含一个整数N,表示吉米邦德的数量以及任务的数量（正好相等,1<=N<=20）。
接下来N行，每行包含N个0到100之间整数，第i行的第j个数Aij表示吉米邦德i完成任务j成功的概率为Aij%

Output

输出所有任务成功完成最大的概率，结果保留6位小数。

Sample Input

输入1：
2 
100 100 
50 50 

输入2：
2 
0 50 
50 0 

输入3：
3 
25 60 100 
13 0 50 
12 70 90 

Sample Output

输出1：
50.000000

输出2：
25.000000

输出3：
9.100000

思路

比赛时用暴力+O2优化过了55

#pragma GCC optimize(2)
#include<cstdio>
int n,x;
double a[25][25];
double ans=0;
bool f[25];
void dfs(int i,int x,double res)
{
	res*=a[i][x]*1.0;
	if(a[i][x]==0)return;
	if(i==n)
	{
		if(res>ans)ans=res;
		return;
	}
	f[x]=1;
	for(int k=1;k<=n;k++)
		if(!f[k])dfs(i+1,k,res);
	f[x]=0;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			scanf("%d",&x);
			a[i][j]=x*1.0/100;
		}
	for(int i=1;i<=n;i++)
		dfs(1,i,1);
	printf("%.6f",ans*100);
}

后来知道正解是状压dp

设f[i]为n个吉米中用1表示接了任务，0表示没接任务的最大概率（i表示为二进制数）

状态转移方程式：f[i]=max{f[i]-2^j *a[j]}

代码

#pragma GCC optimize(2)
#include<cstdio>
#include<algorithm>
using namespace std;
int n;
double a[25][25];
double f[2097152];//2的21次方 
int pow2[25],b[25];
int main()
{
	scanf("%d",&n);
	pow2[1]=1;
	for(int i=1;i<=n;i++)
	{
		pow2[i+1]=pow2[i]*2;
		for(int j=1;j<=n;j++)
		{
			scanf("%lf",&a[i][j]);
			a[i][j]/=100;
		}
	}
	f[0]=100;
	int cnt;
	for(int i=1;i<=pow2[n+1]-1;i++)
	{
		cnt=0;
		for(int j=1;j<=n;j++)
			if(i&pow2[j])b[++cnt]=j;
		for(int j=1;j<=cnt;j++)
			f[i]=max(f[i],f[i-pow2[b[j]]]*a[b[j]][cnt]);
	}
	printf("%.6f",f[pow2[n+1]-1]);
}