题解：P13368 [GCJ 2011 #1B]

#include<iostream>
#include<iomanip>
using namespace std;
int main(){
	int T;
	cin>>T;
	for(int t=1;t<=T;t++){
		int N;
		cin>>N;
		char g[105][105];
		for(int i=0;i<N;i++)cin>>g[i];
		int w[105]={0},c[105]={0};
		for(int i=0;i<N;i++){
			for(int j=0;j<N;j++){
				if(g[i][j]!='.'){
					c[i]++;
					if(g[i][j]=='1')w[i]++;
				}
			}
		}
		double a[105],b[105],d[105],r[105];
		for(int i=0;i<N;i++)a[i]=(double)w[i]/c[i];
		for(int i=0;i<N;i++){
			int k=0;
			double p=0.0;
			for(int j=0;j<N;j++){
				if(g[i][j]!='.'){
					int x=w[j];
					int y=c[j]-1;
					if(g[j][i]=='1')x--;
					p+=(double)x/y;
					k++;
				}
			}
			b[i]=p/k;
		}
		for(int i=0;i<N;i++){
			int k=0;
			double p=0.0;
			for(int j=0;j<N;j++){
				if(g[i][j]!='.'){
					p+=b[j];
					k++;
				}
			}
			d[i]=p/k;
		}
		for(int i=0;i<N;i++)r[i]=0.25*a[i]+0.5*b[i]+0.25*d[i];
		cout<<"Case #"<<t<<":"<<endl;
		for(int i=0;i<N;i++)cout<<fixed<<setprecision(12)<<r[i]<<endl;
	}
	return 0;
}