uva8455 Twice Equation（pell方程+高精度）

题目好神啊orz
给sbw巨佬跪烂
求$n(n+1)=2m(m+1)$，使得有正整数解m的最小的n，要求n>=l。
我们化一下这个等式：
$4n^2+4n+2=8m^2+8m+2$
$(2n+1)^2+1=2(2m+1)^2$
设$x=2n+1,y=2m+1$
$x^2-2y^2=-1$
这就是个pell方程的形式2了。
我们有一组基本解$x_0=1,y_0=1$，可以递推得到$x_n=6x_{n-1}-x_{n-2},x0=1,x1=7$
需要一系列高精。
我的这份代码在计蒜客上a了。然而uva t了。我也不知道为什么qaq

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 200
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,tst,tot=250;
char s[N];
struct bigint{
    int a[N],n;
    bigint(){memset(a,0,sizeof(a));n=0;}
    friend bigint operator*(bigint a,int x){
        bigint res;res.n=a.n;
        for(int i=1;i<=res.n;++i) res.a[i]=a.a[i]*x;
        for(int i=1;i<=res.n;++i) res.a[i+1]+=res.a[i]/10,res.a[i]%=10;
        while(res.a[res.n+1]) res.n++,res.a[res.n+1]+=res.a[res.n]/10,res.a[res.n]%=10;
        return res;
    }friend bigint operator-(bigint a,bigint b){
        bigint res;res.n=a.n;
        for(int i=1;i<=res.n;++i) res.a[i]=a.a[i]-b.a[i];
        for(int i=1;i<=res.n;++i) if(res.a[i]<0) res.a[i]+=10,res.a[i+1]--;
        while(!res.a[res.n]) --res.n;return res;
    }friend bigint operator/(bigint a,int x){
        bigint res;res.n=a.n;int tmp=0;
        for(int i=res.n;i>=1;--i){
            tmp=tmp*10+a.a[i];res.a[i]=tmp/x;tmp%=x;
        }while(!res.a[res.n]) --res.n;return res;
    }friend bool operator<(bigint a,bigint b){
        if(a.n<b.n) return 1;
        if(a.n>b.n) return 0;
        for(int i=a.n;i>=1;--i){
            if(a.a[i]<b.a[i]) return 1;
            if(a.a[i]>b.a[i]) return 0;
        }return 0;
    }
}a[310],b;
int main(){
//  freopen("data.in","r",stdin);
//  freopen("a.out","w",stdout);
    a[0].n=1;a[0].a[1]=1;a[1].n=1;a[1].a[1]=7;
    for(int i=2;i<=tot;++i) a[i]=a[i-1]*6-a[i-2];
    for(int i=1;i<=tot;++i) a[i]=a[i]/2;
    while(~scanf("%d",&tst)){ 
        while(tst--){
            scanf("%s",s+1);
            n=strlen(s+1);b.n=n;
            for(int i=1;i<=n;++i) b.a[i]=s[n-i+1]-'0';
            for(int i=1;i<=tot;++i){
                if(a[i]<b) continue;
                for(int j=a[i].n;j>=1;--j) putchar(a[i].a[j]+'0');puts("");break;
            }
        }
    }return 0;
}