矩阵乘+高斯消元专题小结【高斯消元】

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=20851#overview

I-P

其它题太水了就不写了，这里说一下O和P。

【我太懒不想写】

OP这两个题是高斯消元求期望的题。

对于P题，求从@点到$点的步数期望，先从@点BFS一遍找到所有可以到达的点，对于可达点每个点是一个变量。

因为终点有任意多而起点只有一个，这里可以反过来求从任意一个终点走到起点的期望。

于是对于每一个终点的期望E=0。

对于其它任意一点x，走到这个点的期望 EX=sigma(ENEXT1+……+ENEXTK)/K+1

ENEXT是与X邻接并且可达的点。

这样建立方程组就可以求出起点到终点的期望。

对于O题，这题是走到头再反回来，于是可以把数轴翻转一下再接上去，比如0 1 2 …… n-1 n-2 …… 1

这样共有2*n-2个变量，于是只要向一个方向走就可以了，同样，终点有两个点t和n-t，还是用和上一题一样的方法列方程求解。

下面送上P题的代码

int n,m;
int vis[M][M],poi;
char s[M][M];
double a[N][N];
int dx[]={0,1,0,-1};
int dy[]={1,0,-1,0};
bool bfs()
{
	queue<pair<int,int> > Q;
	int _x,_y,flag=0;
	poi=1;
	for(int i=0;i<n;i++)
	for(int j=0;j<m;j++)
		if(s[i][j]=='@')_x=i,_y=j;
	Q.push(make_pair(_x,_y));
	vis[_x][_y]=1;
	while (! Q.empty())
	{
		int u=Q.front().first,v=Q.front().second;
		Q.pop();
		for(int i=0;i<4;i++)
		{
			int x=u+dx[i],y=v+dy[i];
			if(x>=0 && x<n && y>=0 && y<m)
			if(s[x][y]!='#')
			if(! vis[x][y])
			{
				Q.push(make_pair(x,y));
				vis[x][y]=++poi;
				if(s[x][y]=='$') flag=1;
			}
		}
	}
	return flag;
}
int dcmp(double x)
{
	if(x>EPS) return 1;
	if(x<-EPS) return -1;
	return 0;
}
int guass(double a[][N],int n,int m)
{
	int i,j,k,r;
	double tmp;
	for(i=r=1;i<m;i++)
	{
		for(j=r;j<=n && !dcmp(a[j][i]);j++);
		if(j>n) continue;
		if(j>r) for(k=1;k<=m;k++) swap(a[j][k],a[r][k]);
		for(j=1;j<=n;j++)
		if(j!=r && dcmp(a[j][i]))
		{
			tmp=a[j][i]/a[r][i];
			for(k=1;k<=m;k++)
				a[j][k]-=tmp*a[r][k];
		}
		r++;
	}
	for(i=n;i>0;i--)
	if(dcmp(a[i][m]))
	{
		for(j=1;j<m && !dcmp(a[i][j]);j++);
		if(j>=m) return -1;
	}
	return r-1;
}
int main()
{
	while (scanf("%d%d",&n,&m)!=-1)
	{
		for(int i=0;i<n;i++) scanf("%s",s[i]);
		memset(vis,0,sizeof(vis));
		if(! bfs())
		{
			printf("-1\n");
			continue;
		}
		memset(a,0,sizeof(a));
		for(int i=0;i<n;i++)
		for(int j=0;j<m;j++)
		if(vis[i][j])
		{
			if(s[i][j]=='$')
			{
				a[vis[i][j]][vis[i][j]]=1;
				continue;
			}
			int ss=0;
			for(int k=0;k<4;k++)
			{
				int x=i+dx[k],y=j+dy[k];
				if(x>=0 && x<n && y>=0 && y<m)
				if(vis[x][y])
				{
					ss++;
					a[vis[i][j]][vis[x][y]]=-1;
				}
			}
			a[vis[i][j]][vis[i][j]]=a[vis[i][j]][poi+1]=ss;
		}
		if(guass(a,poi,poi+1)==-1) printf("-1\n");
		else printf("%lf\n",a[1][poi+1]/a[1][1]);
	}
	return 0;
}