求电阻

电阻
【问题描述】
物理课上老师让小 w 测一个电路板的电阻，小 w 用欧姆表很快就得到了答案。但是，
欧姆表巨大的误差是令小 w 无法容忍的，于是他决定用数学方法准确算出总电阻。
电路板一共有 N 个连接点，M 个电阻，小 w 需要算出 1 号连接点和 N 号连接点的总电
阻。
【输入格式】
多组数据，输入直到文件结束
每组数据第一行两个整数 N，M
接下来 M 行，每行三个非负整数 X，Y，R，表示电阻连接的两接点和阻值（X 与 Y 可能
相等）。
【输出格式】
每组数据输出一行，一个实数，四舍五入到小数点后两位。
【输入输出样例】
输入：

2 1

1 2 1
输出：
1.00
【数据范围】
对于 20%的数据 N=2
对于另外 30%的数据 M=N-1 且任意两个连接点间的电阻不是无穷大。
对于 100%的数据保证 1<N<=100

数据保证不出现电阻无穷大的情况。

分析：

因为电脑没有欧姆表，我们只能用伏安法解决问题。

我们设点1电势为1V，点n电势为0V，然后根据每个节点流入的电流等于流出的电流和欧姆定律，列出方程组，然后高斯消元即可，求出每个点的电势，然后欧姆定律。

不知道电势也没有关系，其实电势差值就是电压，这样就明白了吧。

其实高斯消元并不是非常复杂。只需自己画一个梯形矩阵就明白了。

参考程序：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#include<cmath>
using namespace std;
const int maxn=150;
bool vis[maxn];
int INF,n,m;
double f[maxn][maxn],s[maxn][maxn];
double a[maxn][maxn];
bool flag[maxn];
queue<int>Q;
void bfs(int x){
	memset(vis,0,sizeof(vis));
	Q.push(x);
	vis[x]=true;
	while (!Q.empty()){
		int u=Q.front();Q.pop();
		for (int i=1;i<=n;i++)
			if (i!=u && f[u][i]!=INF){
				if (u==1)flag[i]=true;
				a[i][u]=(double)1/f[u][i];
				a[u][i]=(double)-1/f[u][i];
				if (!vis[i]){
					vis[i]=true;
				    Q.push(i);
				}
			}
	}
}
int main(){
	freopen("resistor.in","r",stdin);
	freopen("resistor.out","w",stdout);
	while (scanf("%d%d",&n,&m)==2){
	memset(f,0x3f,sizeof(f));
	memset(s,0,sizeof(s));
	INF=f[0][0];
	for (int i=0;i<m;i++){
		int u,v,w;
		scanf("%d%d%d",&u,&v,&w);
		s[u][v]=s[u][v]+(double)1/w;
		if (u==v)s[u][v]=0;
		s[v][u]=s[u][v];
	}
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n;j++)
			f[i][j]=1/s[i][j];
	memset(a,0,sizeof(a));
	bfs(1);
	for (int i=1;i<=n;i++){
		for (int j=1;j<=n;j++)
			if (j!=i){
				if (a[i][j]<0)a[i][j]=-a[i][j];
					a[i][i]-=a[i][j];
		    }
	}
	for (int i=1;i<n;i++){
		a[i][n]=0;
		a[i][n+1]=-a[i][1];
		a[i][1]=0;
	}
	/*for (int i=1;i<=n;i++){
		for (int j=1;j<=n+1;j++)
			if (a[i][j]!=0)printf("a[%d,%d]= %.2f ",i,j,a[i][j]);
		printf("\n");
	}*/
	for (int i=2;i<=n-1;i++){
		int r=i;
		for (int j=i+1;j<=n-1;j++)
			if (fabs(a[j][i])>fabs(a[r][i]))r=j;
		if (a[r][i]==0)continue;
		if (r!=i)for (int j=2;j<=n+1;j++)swap(a[r][j],a[i][j]);
		for (int k=i+1;k<=n-1;k++){
			double f=a[k][i]/a[i][i];
			for (int j=i;j<=n+1;j++)
				a[k][j]-=f*a[i][j];
		}	
	}
	/*printf("Finish------->\n");
	for (int i=1;i<=n;i++){
		for (int j=1;j<=n;j++)
			if (a[i][j]!=0)printf("a[%d,%d]= %.2f ",i,j,a[i][j]);
		printf("\n");
	}*/
	for (int i=n-1;i>=2;i--){
		//printf("%d %.3f\n",i,a[i][i]);
		for (int j=i+1;j<=n-1;j++)
			a[i][n+1]-=a[j][n+1]*a[i][j];
		a[i][n+1]/=a[i][i];
	}
	double I=0;
	for (int i=2;i<=n;i++){
		//printf("%d %d %.3f\n",i,flag[i],a[i][n+1]);
		if (flag[i])I+=(double)((double)1-a[i][n+1])/f[1][i];
	}
	printf("%.2f\n",1/I);
    }
	return 0;
}