[数学题] pku 3472 Holey Square Tiling

 

首先介绍本题的解题思想:分类思想,化繁为简。善于从复杂问题中抽象出自己熟悉的模型进行解决。

http://www.math.unam.mx/EMIS/journals/JIS/VOL7/Tauraso/tauraso3.pdf的解题报告对我的作用很大。刚开始我按照poj monthly contest解题报告的思想，列出了本题的递推关系，并自己进行程序的编写，发现很是麻烦，并且提交给了我一个超时的警告！无奈之下我只能另寻求其他的方法，想了好久没有什么好的进展，最后在师兄的指导下，我参照了上文的解题报告，利用了论文中关于famous Fibonacci identities的性质将我的递推公式进行了简化      // 4*( 2 * f(n-1)^2 +(-1)^(n-1) )^2   //多n分奇或偶讨论

   但是运行的结果仍然是错误的，调试了好久也没有什么好的结果，经过长达几个小时的折磨，无意中我将short改为int发现结果就正常了，原来都是越界惹的祸啊！提交成功，8000多ms，显然一般比赛的时候这样的程序很难能过，于是我开始着手多高精度算法进行改进，的确，程序运行的时间有了很大的提高。在原来10进制的基础上我用了1000进制，程序运行的时间就减少为1000多ms,的确算法的力量是无穷的！

   多余1000进制在输出的时候应该注意，第一个用printf("%d"),以后的必须用printf("%03d");这里防止出现错误。

   以下是我的高精度算法改进的部分：

const int M=1000;

/高精度乘法
void mult(int *a,int *b,int *c)
{//c=a*b;
  int i,j,m,n,k;
  m=a[0];n=b[0];
  k=n+m-1;
  for(i=0;i<=k;i++)c[i]=0;
  for(i=1;i<=m;i++) {
   if(a[i])for(j=1;j<=n;j++)if(b[j])c[i+j-1]+=a[i]*b[j];
  }
  for(i=k;i>=1;i--)if(c[i]>=M){  c[i-1]+=c[i]/M;c[i]%=M;} 
  if(c[0]>0)
  {k++; for(i=k;i>=1;i--)c[i]=c[i-1];   }c[0]=k;
}

//高精度加法

void add(int *a,int *b,int *c)
{ // c=a+b;
 int m,n,i,j,k;
  m=a[0]; n=b[0];
  if(m<n)k=n; else k=m;
  for(i=0;i<=k;i++) c[i]=0;
  j=k;
  for(i=n;i>=1;i--,j--) c[j]=b[i]; 
  j=k;
  for(i=m;i>=1;i--,j--) c[j]+=a[i];   
  for(i=k;i>=1;i--)
  { if(c[i]>=M) { c[i-1]++; c[i]-=M; }    //
  }
  if(c[0])
  {k++; for(i=k;i>=1;i--) c[i]=c[i-1]; }    
  c[0]=k; 
}