acdream 1200 Jump! Jump! Jump! 概率dp+高斯消元

题意：有n个荷叶形成一个圈，有两只青蛙分别在两个荷叶上，初始位置相距d。每只青蛙在单位时间只能往左右相邻位的荷叶上

跳。不过他们不能控制自己，即往左往右分别有1／2的概率，问两只青蛙最后在跳到同一个荷叶上所需时间的期望。

思路：设dp[ i ]为两者从相距d到达目标状态(落在同一个荷叶)时间的期望。dp［i］=1/4 dp[ i + 2 ]+1/4 dp［i - 2］+1/2 dp［i］+1， 即

为转移方程。不过需要处理下每个距离值，因为有可能相距大于n/2,所以要取x,n-x的最小值。同时因为有的距离状态是无法到达的，

所以需要dfs对出现的距离状态重新编码建立方程，最后高斯消元得到答案，从起始状态无法到达距离为零的状态或高斯消元无解输

出Never Give Up。有一点需要注意的是：我们建方程的时候针对0也建立了方程，但是距离为0的这种状态已经无需转移，而是直接

dp[ 0 ]= 0即可。详见代码：

// file name: acdream1200.cpp //
// author: kereo //
// create time:  2014年09月08日 星期一 22时22分08秒 //
//***********************************//
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<stack>
#include<cmath>
#include<string>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN=100+100;
const int inf=0x3fffffff;
const int mod=1000000000+7;
const double eps=1e-9;
#define L(x) (x<<1)
#define R(x) (x<<1|1)
int n,d,equ,var,cnt;//方程数,var变量数
double a[MAXN][MAXN];
int num[MAXN];
int get_id(int x){
	x=(x%n+n)%n;
	if(x>n-x) x=n-x;
	return x;
}
int Gauss(){
    int i,col=0;
    for(i=0;i<equ && col<var;i++,col++){
        int r=i;
        for(int j=i+1;j<equ;j++)
            if(fabs(a[j][col])>fabs(a[r][col]))
                    r=j;
        if(fabs(a[r][col])<eps) return 0;
        if(r!=i){
            for(int j=col;j<=var;j++)
                swap(a[r][j],a[i][j]);
        }
        for(int j=col+1;j<=var;j++) a[i][j]/=a[i][col];
        a[i][col]=1;
        for(int k=0;k<equ;k++) if(k!=i){
            for(int j=col+1;j<=var;j++) a[k][j]-=a[i][j]*a[k][col];
            a[k][col]=0;
        }
    }
    return 1;
}
void dfs(int x){
	x=get_id(x);
	num[x]=cnt++;
	int x1=get_id(x+2);
	if(num[x1]==-1) dfs(x1);
	int x2=get_id(x-2);
	if(num[x2] == -1) dfs(x2);
}
int main()
{
	int kase=0;
	while(~scanf("%d%d",&n,&d)){
		printf("Case #%d: ",++kase);
		cnt=0;
		memset(num,-1,sizeof(num));
		memset(a,0,sizeof(a));
		d=get_id(d);
		dfs(d);
		if(num[0] == -1){
			printf("Never Give Up\n");
			continue;
		}
		for(int i=0;i<n;i++){
			if(num[i] == - 1) continue;
			int pos=num[i];
			a[pos][pos]=2;
			int x=get_id(i-2);
			a[pos][num[x]]-=1;
			x=get_id(i+2);
			a[pos][num[x]]-=1;
			a[pos][cnt]=4;
		}
		int x=num[0];
		memset(a[x],0,sizeof(a[x]));
		a[x][x]=1;
		equ=var=cnt;
		if(Gauss())
			printf("%.2lf\n",a[num[d]][cnt]);
		else 
			printf("Never Give Up\n");
	}
	return 0;
}