1369 - Pell数列 --- 高精度算法

** 1369 - Pell数列 --- 高精度算法
** 来源：东方博宜oj   oj.czos.cn
#include<bits/stdc++.h>
using namespace std;
string add(string s1,string s2)
{
    int a[1000],b[1000],c[1000];
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    for(int i=0;i<s1.size();i++)  a[i]=s1[s1.size()-1-i]-'0';
    for(int i=0;i<s2.size();i++)  b[i]=s2[s2.size()-1-i]-'0';
    int len=max(s1.size(),s2.size());
    for(int i=0;i<len;i++)
    {
        c[i]=c[i]+a[i]+b[i];
        if(c[i]>=10)
        {
            c[i+1]=c[i+1]+c[i]/10;
            c[i]=c[i]%10;
        }
    }
    while(c[len]==0&&len>=1) len--;
    string s="";
    for(int i=0;i<=len;i++) s=char(c[i]+'0')+s;
    return s;
}
string mul(string s)
{
    int a[1000];
    memset(a,0,sizeof(a));
    for(int i=0;i<s.size();i++) a[i]=s[s.size()-1-i]-'0';
    int len=s.size()+1;
    for(int i=0;i<s.size();i++)  a[i]=a[i]*2;
    for(int i=0;i<s.size();i++)
    {
        if(a[i]>=10)
        {
            a[i+1]=a[i+1]+a[i]/10;
            a[i]=a[i]%10;
        }
    }
    while(a[len]==0&&len>1) len--;
    string s3;
    for(int i=0;i<=len;i++) s3=to_string(a[i])+s3;
    return s3;
}
int main()
{
    int n;
    string x,y,z;
    cin>>n;
    if(n<=2)
    {
        cout<< n;
        return 0;
    }
    x="1";y="2";
    for(int i=3;i<=n;i++)
    {
        z=add(mul(y),x);
        x=y;
        y=z;
    }
    cout<< z;
    return 0;
}