UVa 104 Arbitrage (DP+最短路+Floyd)

http://www.cnblogs.com/scau20110726/archive/2012/12/26/2834674.html

/*
DP构建（仿照Floyd）

这题说难挺难的，现在想想好像也不太难，是因为想复杂了。其实这题的要求不算多，它仅仅是要求最后转换回到出发时的那种硬币并且获利大于0.01，而转换次数不能超过n次，转换次数越少越好，然后输出最短的路径。注意：题目并没有说要你获利最大，只是是获利>=0.01即可，另外题目对路径并没有什么特殊要求，仅仅是要求不能超过n，在路径长度相同的情况任意一条都可以的。

所以我们就抓住一点：转换次数！！！！！！

也就是说每转换一次，我们就去找所有的货币，看他们回到自己的时候那些货币的获利已经大于0.01了，如果有货币的获利已经大于0.01了，那么我们的算法就结束了，就输出这种货币的转换路径（就算有多种货币的获利都大于0.01都好，输出哪种货币都行，因为题目也没有要求）

我们应该写一个最外层的循环，是表示第几次转换的，这个循环不能大于n，因为最多只能转换n次不管成功还是失败都要结束

然后，用转换次数来做最外层的循环这个策略是正确的，为什么呢？假设当前是第i+1次转换，那么是需要用到第i次转化的信息的！！！

为什么呢？我们可以回顾一下Floyd的本质。

原本的Floyd是一个三维数组   dp[k][i][j]=max{ dp[k-1][i][j]  , dp[k-1][i][k]+dp[k-1][k][j] }

我们平时写的Floyd都是二维的，其实是三维的压缩回来的，没什么能压缩掉最高的那一维呢？单单看状态转移方程都可能看出来，当前要构建的dp[k][i][j]其实只和之前的dp[k-1][][]有关，也就是现在的信息完全是基于前一次的信息推导出来的。这就像我们的滚动数组一样

 

所以我们来定义个状态dp[i][j][s]  i和j就是第i种和第j种货币，s是转换次数，s从1开始到n。也就是说，在第s次转换中，试图从i转化内j

那么状态转移方程是什么呢？ dp[i][j][s]=max{ dp[i][j][s-1] , dp[i][k][s-1]*g[k][j]}

g[i][j]是输入的数据也就是两两货币之间的转化率，其实可以保存在dp[i][j][0]里面

那么这个状态转移方程是什么意思呢？就是说在前一次转化后，现在想从i转换到j，试图通过k，看能不能让值更大，能的话就变大（变大一定不吃亏，因为题目也没要求货币变成多少，只是要求要大于0.01，所以当前要是能变大的话那当然变大）
*/
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<climits>
#include<cstring>
#include<string>
#include<set>
#include<map>
#include<queue>
#include<vector>
using namespace std;

const int maxn=20+10;
int n;
double d[maxn][maxn][maxn];
int p[maxn][maxn][maxn];

void print_path(int i,int j,int s)
{
	if(s==0)
		printf("%d",i);
	else
	{
		print_path(i,p[i][j][s],s-1);
		printf(" %d",j);
	}
}

void DP()
{
	int s,m;
	for(s=2;s<=n;s++)
	{
		for(int k=1;k<=n;k++)
			for(int i=1;i<=n;i++)
				for(int j=1;j<=n;j++)
					if(d[i][j][s]<d[i][k][s-1]*d[k][j][1])
					{
						d[i][j][s]=d[i][k][s-1]*d[k][j][1];
						p[i][j][s]=k;
					}
		int i;
		for(i=1;i<=n;i++)
		{
			if(d[i][i][s]>1.01)
			{
				m=i;
				break;
			}
		}
		if(i>n)continue;
		else break;
	}
	if(s>n)printf("no arbitrage sequence exists\n");
	else
	{
		print_path(m,m,s);
		printf("\n");
	}
}

int main()
{
	while(~scanf("%d",&n))
	{
		memset(d,0,sizeof(d));
		memset(p,0,sizeof(p));
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
			{
				if(i==j)d[i][j][1]=1.0;
				else scanf("%lf",&d[i][j][1]);
				p[i][j][1]=j;
			}
        DP();
	}
}