bzoj2467 生成树

传送门

生成树计数裸题，因为模数不是质数所以要用辗转相除的高斯消元。

de了很久的bug然后发现一个变量A赋值后下一行又申明了一个新的临时变量A（：

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<ctime>
const int mod=;
const int N=;
typedef long long LL;
using namespace std;
int T,n,g[N][N],ans;

template<typename T> void read(T &x) {
    T f=; x=; char ch=getchar();
    while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
    if(ch=='-') f=-,ch=getchar();
    for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
}

int guess(int n) {
    int res=,f=;
    for(int i=;i<=n;i++) {
        for(int j=i+;j<=n;j++) {
            int A=g[i][i],B=g[j][i];
            while(B) {
                int t=A/B; A=A%B; swap(A,B);
                for(int k=i;k<=n;k++)
                    g[i][k]=(g[i][k]-t*g[j][k]%mod+mod)%mod;
                for(int k=i;k<=n;k++)
                    swap(g[i][k],g[j][k]);
                f=-f;
            }
        }
        if(!g[i][i]) return ;
        res=(res*g[i][i])%mod;
    }
    if(f==-) res=(mod-res)%mod;
    return res;
}

int main() {
    read(T);
    while(T--) {
        read(n);
        if(n==) printf("4\n");
        else if(n==) printf("40\n");
        else {
            int tot=n;
            memset(g,,sizeof(g));
            for(int i=;i<=n;i++) {
                if(i==) g[i][n]--,g[i][]--;
                else if(i==n) g[i][]--,g[i][i-]--;
                else g[i][i-]--,g[i][i+]--;
                g[i][i]+=;
            }
            for(int i=;i<=tot;i++) {
                int x=n+,y=n+,z=n+,w=(i==n)?:i+;
                g[x][i]--; g[i][x]--;
                g[x][y]--; g[y][x]--;
                g[y][z]--; g[z][y]--;
                g[z][w]--; g[w][z]--;
                g[x][x]+=; g[y][y]+=; g[z][z]+=;
                n+=;
            }
            for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
                (g[i][j]+=mod)%=mod;
        }
        ans=guess(n-);
        printf("%d\n",ans);
        }
    }
    return ;
}