bzoj1002 轮状病毒 暴力打标找规律/基尔霍夫矩阵+高斯消元

//暴力代码
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<set>

using namespace std;

typedef long long ll;
typedef long long LL;
typedef pair<int,int> pii;
const int inf = 0x3f3f3f3f;
const int maxn = 200+10;
const ll mod = 1e9+9;
int ans,cnt,n;
int head[maxn],fa[maxn],vis[maxn];
struct Table{
    int to,next;
    bool ban;
}table[maxn<<2];
void add(int x,int y){
    table[cnt].to=y;
    table[cnt].next=head[x];
    head[x]=cnt++;
}
bool bfs()
{
    memset(fa,-1,sizeof(fa));
    memset(vis,0,sizeof(vis));
    queue<int>q;q.push(0);vis[0]=true;int tot=1;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u];i!=-1;i=table[i].next){
            if(!table[i].ban){
                int v=table[i].to;
                if(v==fa[u]) continue;
                if(vis[v]) return 0;
                fa[v]=u;
                vis[v]=1;
                q.push(v);
                tot++;
            }
        }
    }
    if(tot<=n) return 0;
    else return 1;
}
void dfs(int x){
    if((x<<1)>=cnt){
        if(bfs()){
            ++ans;
        }
        return;
    }
    table[x<<1].ban=table[x<<1|1].ban=0;
    dfs(x+1);
    table[x<<1].ban=table[x<<1|1].ban=1;
    dfs(x+1);
}
int main(){
    int num;
    scanf("%d",&num);
    for(int i=1;i<=num;i++){
        ans=0;cnt=0;n=i;
        memset(head,-1,sizeof(head));
        for(int j=1;j<=i;j++){
            add(0,j);add(j,0);add(j,j%n+1);add(j%n+1,j);
        }
        dfs(0);
        cout<<ans<<" ";
    }
    return 0;
}
//大数规律代码
#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;

typedef long long ll;
const int inf = 0x3f3f3f3f;
const int maxn = 100+10;

struct BigInteger{
    int arr[maxn];
    int len;
    BigInteger(){
        memset(arr,0,sizeof(arr));
        len=0;
    }
};
BigInteger Mul(BigInteger a,int k){
    for(int i=1;i<=a.len;i++){
        a.arr[i]*=k;
    }
    for(int i=1;i<=a.len;i++){
        a.arr[i+1]+=a.arr[i]/10;
        a.arr[i]%=10;
    }
    if(a.arr[a.len+1]!=0) a.len++;
    return a;
}
BigInteger Sub(BigInteger a,BigInteger b){
    a.arr[1]+=2;
    int j=1;
    while(a.arr[j]>=10){
        a.arr[j+1]++;
        a.arr[j]%=10;
        j++;
    }
    //！！！
    a.len=max(a.len,j);
    for(int i=1;i<=a.len;i++){
        a.arr[i]-=b.arr[i];
        if(a.arr[i]<0){
            a.arr[i]+=10;
            a.arr[i+1]--;
        }
    }
    while(a.arr[a.len]==0) a.len--;
    return a;
}
int main(){
    BigInteger f[maxn];
    f[1].arr[1]=1;f[1].len=1;
    f[2].arr[1]=5;f[2].len=1;
    int n;
    scanf("%d",&n);
    for(int i=3;i<=n;i++){
        f[i]=Sub(Mul(f[i-1],3),f[i-2]);
    }
    for(int i=f[n].len;i>0;i--){
        printf("%d",f[n].arr[i]);
    }
    printf("\n");
    return 0;
}