bzoj 1002 [FJOI2007]轮状病毒——打表找规律

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1002

看 Zinn 的博客：https://www.cnblogs.com/Zinn/p/9252831.html

时隔六个月，自己终于也 A 了这道题。

#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 ;
}