随机程序(数学期望,高斯约当消元法)

const int N = 105;
double A[N][N];
int n;
void gauss_jordan()
{
	f(i, 0, n - 1)
	{
		int r = i;
		f(j, i + 1, n - 1)
			if (fabs(A[j][i]) > fabs(A[r][i]))r = j;
		if (fabs(A[r][i]) < eps)continue;
		if (r != i)f(j, 0, n)swap(A[r][j], A[i][j]);
		f(k, 0, n - 1)
			if (k != i)ff(j, n, i)A[k][j] -= A[k][i] / A[i][i] * A[i][j];
	}
}
vector<int> G[N];
int d[N];
int if_inf[N];
int main()
{
	//freopen("in.txt", "r", stdin);
	int cas = 0;
	while (scanf("%d",&n)==1&&n)
	{
		memset(d, 0, sizeof d);
		f(i, 0, n)G[i].clear();
		int x, y;
		while (scanf("%d%d", &x, &y) == 2 && x)
		{
			x--, y--;
			G[y].emplace_back(x);//记录前驱来计算期望
			d[x]++;
		}
		memset(A, 0, sizeof A);
		f(i, 0, n - 1)
		{
			if (i == 0)A[0][n] = 1;
			A[i][i] = 1;
			for (auto I : G[i])
			{
				int pre_son = d[I];//前驱点的孩子
				A[i][I] -= 1.0 / pre_son;
			}
		}
		gauss_jordan();
		memset(if_inf, 0, sizeof if_inf);
		ff(i, n - 1, 0)
		{
			if (fabs(A[i][i]) < eps&&fabs(A[i][n]) > eps)if_inf[i] = 1;
			f(j, i + 1, n - 1)
				if (fabs(A[i][j]) > eps&&if_inf[j] == 1)if_inf[i] = 1;
			//和无穷变量扯上关系的变量也是无穷的
		}
		int q = in();
		printf("Case #%d:\n", ++cas);
		while (q--)
		{
			int x = in();x--;
			if (if_inf[x] == 1)puts("infinity");
			else printf("%.3lf\n", A[x][x]<eps?0.0:A[x][n] / A[x][x]);
		}
	}
	return 0;
}	
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值