bzoj4666 小z的胡话

题目描述：

bz

题解：

乱搞好题哇。

众所周知斐波那契数列是有循环节的。

我们可以搞出在模$10^x$下与所给得数同余的集合，那么在模$10^{x+1}$下，同余集合一定是原集合及循环若干循环节的大集合的子集。

人话是，i have a xunhuanjie, i have a jihe.EN!another jihe and another xunhuanjie.

然后矩乘求值判定就好了。

代码：

#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
ll MOD;
void Mod(ll&x){if(x>=MOD)x-=MOD;}
ll fastadd(ll a,ll b)
{
    ll t = (long double)a*b/MOD;
    ll ret = a*b-t*MOD;
    while(ret<)ret+=MOD;
    return ret;
}
struct mt
{
    ll s[][];
    void reset(ll k)
    {
        s[][]=s[][]=;
        s[][]=s[][]=k;
    }
    void init()
    {
        s[][]=s[][]=s[][]=;
        s[][]=;
    }
    bool check()
    {
        for(int i=;i<=;i++)for(int j=;j<=;j++)
            if(s[i][j]!=(i==j))return ;
        return ;
    }
    mt operator * (const mt&a)const
    {
        mt ret;ret.reset();
        for(int i=;i<=;i++)
            for(int j=;j<=;j++)
                for(int k=;k<=;k++)
                    Mod(ret.s[i][j]+=fastadd(s[i][k],a.s[k][j]));
        return ret;
    }
};
mt operator ^ (mt&x,ll y)
{
    mt ret;ret.reset();
    while(y)
    {
        if(y&)ret=ret*x;
        x=x*x;y>>=;
    }
    return ret;
}
ll n,L[];
vector<ll>ve[];

int main()
{
    scanf("%lld",&n);
    L[] = ,ve[].push_back(),MOD = ;
    mt bas;bas.reset();
    for(int i=;i<=;i++)
    {
        MOD=MOD*;
        mt now;now.init();now=now^L[i&];
        L[!(i&)]=,ve[!(i&)].clear();
        mt nxt;nxt.reset();
        while(!L[!(i&)]||!nxt.check())
        {
            for(int j=,lim=(int)ve[i&].size();j<lim;j++)
            {
                mt tmp;tmp.init();tmp=tmp^(ve[i&][j]+L[!(i&)]);
                if(tmp.s[][]==n%MOD)ve[!(i&)].push_back(ve[i&][j]+L[!(i&)]);
            }
            L[!(i&)]+=L[i&];nxt=nxt*now;
        }
    }
    if(!ve[].size())puts("-1");
    else printf("%lld\n",ve[][]+);
    return ;
}