France '98

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=30506#problem/H

#include<map>
#include<set>
#include<list>
#include<cmath>
#include<ctime>
#include<deque>
#include<stack>
#include<bitset>
#include<cstdio>
#include<vector>
#include<cstdlib>
#include<cstring>
#include<iomanip>
#include<numeric>
#include<sstream>
#include<utility>
#include<iostream>
#include<algorithm>
#include<functional>

using namespace std ;
char str[ 20 ][ 50 ] ;
double val[ 20 ][ 20 ] ;
double dp[ 20 ][ 5 ] ;
 
int main()
{
	for( int i = 0 ; i < 16 ; ++i )
	{
		scanf( "%s" , str[ i ] ) ;
	}
	for( int i = 0 ; i < 16 ; ++i )
	{
		for( int j = 0 ; j < 16 ; ++j )
		{
			scanf( "%lf" , &val[ i ][ j ] ) ;
			val[ i ][ j ] /= 100 ;
		}
	}
	memset( dp , 0 , sizeof( dp ) ) ;
	for( int i = 0 ; i < 16 ; ++i )
	{
		dp[ i ][ 0 ] = 1 ;
	}
	int len = 1 ;
	for( int i = 1 ; i <= 4 ; ++i )
	{
	
		for( int j = 0 ; j < 16 ; ++j )
		{	
			int st ;
			if( ( j / len ) & 1 )
			{
				st = ( j / len - 1 ) * len ;
			}
			else
			{
				st = ( j / len + 1 ) * len ;
			}
			
			for( int k = st ; k < st + len ; ++k )
			{
				dp[ j ][ i ] += dp[ j ][ i - 1 ] * dp[ k ][ i - 1 ] * val[ j ][ k ] ;
			}
		}
		len = len * 2 ;
	}
	for( int i = 0 ; i < 16 ; ++i )
	{
		printf( "%-11sp=%.2lf%%\n" , str[ i ] , dp[ i ][ 4 ]* 100 ) ; 
	}
	//system( "pause" ) ;
    return 0;
}