uva10828

https://vjudge.net/problem/UVA-10828

裸高斯消元。。。

但是要判无解和无穷解。

当出现一个环时会无解，环上每个点只有一个出边。

#include<bits/stdc++.h>
using namespace std;
const int N = ;
const double eps = 1e-;
int n, q;
double a[N][N], d[N];
vector<int> G[N];
int mark[N];
void build()
{
    a[][n + ] = -1.0;
    for(int i = ; i <= n; ++i)
    {
        a[i][i] = -1.0;
        for(int j = ; j < G[i].size(); ++j)
            a[i][G[i][j]] += 1.0 / d[G[i][j]];
    }
}
void gauss_jordan()
{
    for(int now = ; now <= n; ++now)
    {
        int x = now;
        for(int i = now + ; i <= n; ++i) if(fabs(a[i][now]) > fabs(a[x][now])) x = i;
        if(fabs(a[x][now]) < eps) continue;
        for(int i = ; i <= n + ; ++i) swap(a[now][i], a[x][i]);
        double t = a[now][now];
        for(int i = now; i <= n + ; ++i) a[now][i] /= t;
        for(int i = ; i <= n; ++i) if(i != now)
        {
            double t = a[i][now];
            for(int j = now; j <= n + ; ++j) a[i][j] -= a[now][j] * t;
        }
    }
}
int main()
{
    int T = ;
    while(scanf("%d", &n))
    {
        ++T;
        if(n == ) break;
        memset(d, , sizeof(d));
        memset(a, , sizeof(a));
        memset(mark, , sizeof(mark));
        while()
        {
            int a, b; scanf("%d%d", &a, &b);
            if(a ==  && b == ) break;
            G[b].push_back(a);
            d[a] += 1.0;
        }
        build();
        gauss_jordan();
        for(int i = n; i; --i)
        {
            if(fabs(a[i][i]) < eps && fabs(a[i][n + ]) > eps)
                mark[i] = ;
            for(int j = i + ; j <= n; ++j) if(fabs(a[i][j]) > eps && mark[j]) mark[i] = ;
        }
        scanf("%d", &q);
        printf("Case #%d:\n", T);
        while(q--)
        {
            int x; scanf("%d", &x);
            if(mark[x]) puts("infinity");
            else if(fabs(a[x][x]) < eps) puts("0.000");
            else printf("%.3f\n", fabs(a[x][n + ] / a[x][x]));
        }
        for(int i = ; i <= n; ++i) G[i].clear();
    }
    return ;
}