【gmoj】【高斯消元】 折射伤害

【gmoj】【高斯消元】 折射伤害

题目

---

解题思路

题意
每当受到伤害时就转移给别人
不停的转移

暴力纯模拟直到修改的微乎其微

根据样例可以知道每人最终受到的伤害
由自己原来的和从别人那转移的组成
f1=1+½f2
f2=0+½f1
f3=2
组成了多元一次的方程组
用高斯消元解决

---

代码

#include<iostream>
#include<cstdio>
#include<cmath> 
using namespace std;
double z,p[220],f[220][220];
int n,m,x,y,a[220];
int main()
{
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		f[i][n+1]=a[i];
		f[i][i]=-1;
	}
	scanf("%d",&m);
	for (int i=1;i<=m;i++)
	{
		scanf("%d%d%lf",&x,&y,&z);
		p[x]+=z;
		f[y][x]+=z;
	}
	for (int i=1;i<=n;i++)  //枚举列
	{
		int w=i; 
	    for (int j=i+1;j<=n;j++)
	        if (fabs(f[j][i])>fabs(f[w][i])) w=j;  //找到最大的
	    for (int j=i;j<=n;j++)  //交换，保证最大系数都在f[i][i]上
	         swap(f[i][j],f[w][j]);
	    for (int j=1;j<=n;j++)  //加减消元
	    if (i!=j)
	    {
	        double tmp=1.0*f[j][i]/f[i][i];  //求一个数，让第i行减第j行能消掉第i列
	        for (int k=1;k<=n+1;k++)  
	            f[j][k]-=tmp*f[i][k];  
		} 
	}
	for (int i=1;i<=n;i++)
	    printf("%.6lf\n",abs(f[i][n+1]/f[i][i]*(1-p[i])));
}
//5x+y+6z=25
//3x+2y+z=10
//tmp=0.6
//3x+2y+z-3x-0.6y-3.6z=10-15
//1.4y-2.6z=-5
//消第一项的模拟