#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cmath>

#define db double

using namespace std;

const int N=;

int n;db a[N][N];

db get()

{

  db ret=;

  for(int i=;i<=n;i++)

    {

      for(int j=i+;j<=n;j++)

    if(fabs(a[j][i])>fabs(a[i][i]))

      {

        for(int t=i;t<=n;t++)swap(a[i][t],a[j][t]);

      }

      if(!a[i][i])return ;

      for(int j=i+;j<=n;j++)

    {

      db sl=a[j][i]/a[i][i];

      for(int t=i;t<=n;t++)a[j][t]-=sl*a[i][t];

    }

      ret*=a[i][i];

    }

  return ret;

}

int main()

{

  int lm;scanf("%d",&lm);

  for(n=;n<=lm;n++)

    {

      for(int i=;i<=n;i++)for(int j=;j<=n;j++)a[i][j]=;

      for(int i=;i<=n;i++)

    {

      a[i][i]+=;

      i==?a[i][n]-=:a[i][i-]-=;

      i==n?a[i][]-=:a[i][i+]-=;

      a[i][n+]-=;a[n+][i]-=;

    }

      a[n+][n+]+=n;

      printf("%.0f ",get());

    }

  puts("");

  return ;

}

打表

#include<cstdio>

#include<cstring>

#include<algorithm>

using namespace std;

const int N=,base=1e8;

int n,a[N][N];

void mul(int *a,int *ret,int b)

{

  ret[]=a[];

  for(int i=;i<=a[];i++)

    ret[i]=a[i]*b;

  for(int i=;i<=ret[];i++)

    ret[i+]+=ret[i]/base,ret[i]%=base;

  while(ret[ret[]+])ret[]++,ret[ret[]+]+=ret[ret[]]/base,ret[ret[]]%=base;

}

void dec(int *a,int *b)

{

  for(int i=;i<=a[];i++)

    {

      a[i]-=b[i];

      if(a[i]<)a[i]+=base,a[i+]--;

    }

  while(a[]&&!a[a[]])a[]--;

}

void add(int *a,int b)

{

  a[]+=b;

  for(int i=;i<=a[]&&a[i]>=base;i++)

    a[i+]++,a[i]-=base;

  if(a[a[]+])a[]++;

}

void print(int *a)

{

  printf("%d",a[a[]]);

  for(int i=a[]-;i;i--)printf("%08d",a[i]);puts("");

}

int main()

{

  scanf("%d",&n);

  if(n<=){puts("");return ;}

  a[][]=;a[][]=;a[][]=;a[][]=;

  if(n<=){printf("%d\n",a[n][]);return ;}

  for(int i=;i<=n;i++)

    {

      mul(a[i-],a[i],);

      dec(a[i],a[i-]);add(a[i],);

    }

  print(a[n]);return ;

}