HDU 5351 MZL's Border (规律，大数)

[HDU 5351 MZL's Border]
题意
　　定义字符串$f_1=b，f_2=a，f_i=f_{i-1}f_{i-2}$。

　　对$f_n$的长度为$m$的前缀$s$，

　　求最大的$k$满足$s[1]=s[m-k+1],s[2]=s[m-k+2]...s[k]=s[m]$

---

Solution

　　先列出前几个字符串

　　$f_3=ab,f_4=aba,f_5=abaab,f_6=abaababa,f_7=abaababaabaab$。

　　对于n，m。

　　假设$f_{n-1}$的长度大于等于$m$，那么相当于求n-1，m。

　　假设$f_{n-1}$的小于$m$

　　由斐波那契的性质，答案要么是$m-length(f_{n-1})$,

　　要么是$m-length(f_{n-2})$

　　通过对前几个字符串的分析可以发现，当$m+1<length(f_n)$时答案是后者。

　　于是只要先把前1000个字符串的长度求出来。需要使用高精度。

　　然后直接按照上面的规则判断就行了。

#include <bits/stdc++.h>
#define lson x << 1
#define rson x << 1 | 1
#define LL long long
#define se second
#define fi first
#define LL long long
using namespace std;

const int N = ;

string  fib[N], m, ans;

inline string add ( string a, string b )
{
    if ( a.size() < b.size() ) swap ( a, b );
    int la = a.size(), lb = b.size(), x = , i = ;
    for ( i = ; i < lb; ++i ) {
        x += ( a[i] - '' + b[i] - '' );
        a[i] = '' + x % ;
        x /= ;
    }
    while ( x && i < la ) {
        x += ( a[i] - '' );
        a[i] = '' + x % ;
        x /= ;
        ++i;
    }
    if ( x ) a += '' + x;
    return a;
}

inline string sub ( string a, string b )
{
    int la = a.size(), lb = b.size(), x = , i = ;
    for ( i = ; i < lb; ++i ) {
        if ( a[i] - x >= b[i] ) {
            a[i] = a[i] - x - b[i] + '';
            x = ;
        } else {
            a[i] =  + '' + a[i] - b[i] - x;
            x = ;
        }
    }
    while ( x && i < la ) {
        if ( a[i] != '' ) {
            a[i] = a[i] - ;
            x = ;
        } else {
            a[i] = '';
            ++i;
        }
    }
    i = la - ;
    while ( a[i] == '' && i !=  ) {
        a.erase ( a.begin() + i );
        --i;
    }
    return a;
}

int T, n;

inline bool pd ( string a, string b )
{
    if ( a.size() == b.size() )  {
        for ( int i = a.size() - ; i >= ; --i ) {
            if ( a[i] != b[i] ) return a[i] > b[i];
        }
        return ;
    }
    return a.size() > b.size();
}

void make ( int n, string m )
{
    if ( n <  ) {
        ans = "";
        return ;
    }
    if ( pd ( m, fib[n - ] ) ) {
        string s1 = sub ( fib[n], "" );
        if ( !pd ( m, s1 ) ) {
            ans = sub ( m, fib[n - ] );
        } else {
            ans = sub ( m, fib[n - ] );;
        }
        return ;
    } else make ( n - , m );
}

void MOD ( string s, LL k )
{
    LL tem = ;
    for ( int i = s.size() - ; i >= ; --i ) {
        tem = ( tem *  + s[i] - '' ) % k;
    }
    cout << tem << endl;
}
int main()
{
    fib[] = "";
    fib[] = "";
    for ( int i = ; i <= ; ++i ) {
        fib[i] = add ( fib[i - ] , fib[i - ] );
    }
    cin >> T;
    while ( T-- ) {
        cin >> n >> m;
        reverse ( m.begin(), m.end() );
        make ( n, m );
        MOD ( ans,  );
    }
}