• 【SPOJ 104】HIGH


    题目描述
    In some countries building highways takes a lot of time… Maybe that’s because there are many possiblities to construct a network of highways and engineers can’t make up their minds which one to choose. Suppose we have a list of cities that can be connected directly. Your task is to count how many ways there are to build such a network that between every two cities there exists exactly one path. Two networks differ if there are two cities that are connected directly in the first case and aren’t in the second case. At most one highway connects two cities. No highway connects a city to itself. Highways are two-way.

    输入输出格式
    输入格式:
    The input begins with the integer t, the number of test cases (equal to about 1000). Then t test cases follow. The first line of each test case contains two integers, the number of cities (1<=n<=12) and the number of direct connections between them. Each next line contains two integers a and b, which are numbers of cities that can be connected. Cities are numbered from 1 to n. Consecutive test cases are separated with one blank line.

    输出格式:
    The number of ways to build the network, for every test case in a separate line. Assume that when there is only one city, the answer should be 1. The answer will fit in a signed 64-bit integer.

    输入输出样例
    输入样例#1:
    4
    4 5
    3 4
    4 2
    2 3
    1 2
    1 3

    2 1
    2 1

    1 0

    3 3
    1 2
    2 3
    3 1
    输出样例#1:
    8
    1
    1
    3

    先用高斯消元转化为上三角,后求对角线之积算出行列式绝对值
    code:

    //Menteur_Hxy
    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #include<cmath>
    using namespace std;
    #define ll long long
    
    ll rd() {
        ll x=0,fla=1; char c=' ';
        while(c<'0' || c>'9') {c=getchar();if(c=='-') fla=-fla;}
        while(c>='0' && c<='9') x=x*10+c-'0',c=getchar();
        return x*fla;
    }
    
    const double eps=1e-12;
    int T,n,m;
    double a[3010][3010];
    
    void gauss() {
        for(int i=1;i<=n;i++) {
            int mx=i;
            for(int j=i+1;j<=n;j++) mx=(a[mx][i]-a[j][i]<-eps)?i:mx;
            if(mx!=i) swap(a[i],a[mx]);
            if(!a[i][i]) {printf("0
    ");return ;}
            for(int j=i+1;j<=n;j++) {
                double t=a[j][i]/a[i][i];
                for(int k=i;k<=n+1;k++)
                    a[j][k]-=t*a[i][k];
            }
        }
        double ans=1;
        for(int i=1;i<=n;i++) ans=ans*a[i][i];
        printf("%.0f
    ",fabs(ans));
    }
    
    int main() {
        T=rd();
        while(T--) {
            memset(a,0,sizeof a);
            n=rd()-1,m=rd();
            for(int i=1;i<=m;i++) {
                int u=rd(),v=rd();
                a[u][u]++,a[v][v]++;
                a[u][v]--,a[v][u]--;
            }
            gauss();
        }
        return 0;
    }
    版权声明:本文为博主原创文章,未经博主允许不得转载。 博主:https://www.cnblogs.com/Menteur-Hxy/
  • 相关阅读:
    语料和文本处理
    seq2seq+torch7聊天机器人bug处理
    unity3d inputfield标签控制台打印object
    多种语言tcp编程
    处理json中的空格
    安卓无法访问Azure服务器和微软API
    Xamarin/Unity3d无法访问Azure服务器或者微软API
    unity3d C# soket客户端接受失败
    unity3d之public变量引发错误
    unity3d,java,c#,python,rospy的socket通信测试
  • 原文地址:https://www.cnblogs.com/Menteur-Hxy/p/9247975.html
Copyright © 2020-2023  润新知