bzoj 4294 [PA2015]Fibonacci 数学

本文介绍了一种使用矩阵快速幂的方法来解决斐波那契数列取模的问题,并通过深度优先搜索结合剪枝策略来枚举每一位数,最终找到满足特定条件的斐波那契数列项。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

斐波那契取模10x的循环节是610x(别问我怎么证。。。)
求第x项取模可以矩乘。
因此可以深搜枚举每一位然后验证减枝。

#include <bits/stdc++.h>
using namespace std;
#define ull unsigned long long
char s[21];
int n;
ull ten[21],v,mod;
ull mul(ull x,ull y)
{
    ull ret=0;
    while(y)
    {
        if(y&1)ret=(ret+x)%mod;
        x=(x+x)%mod;y>>=1;
    }
    return ret;
}
struct matrix 
{
    ull w[2][2];
    friend matrix operator * (const matrix &r1,const matrix &r2)
    {
        matrix ret;
        memset(&ret,0,sizeof(ret));
        for(int i=0;i<2;i++)
            for(int j=0;j<2;j++)
                for(int k=0;k<2;k++)
                    (ret.w[i][j]+=mul(r1.w[i][k],r2.w[k][j]))%=mod;
        return ret;
    }
}mat;
matrix qpow(matrix x,ull y)
{
    matrix ret;
    ret.w[0][0]=1;ret.w[1][1]=1;
    ret.w[0][1]=ret.w[1][0]=0;
    while(y)
    {
        if(y&1)ret=ret*x;
        x=x*x;y>>=1;
    }
    return ret;
}
ull fib(ull x)
{
    mat.w[0][0]=0;mat.w[0][1]=1;
    mat.w[1][0]=1;mat.w[1][1]=1;
    mat=qpow(mat,x);
    return mat.w[0][1];
}
void dfs(ull x,int y)
{
    mod=ten[y];
    if(fib(x)!=v%ten[y])return;
    if(y==n){printf("%llu\n",x+6*ten[18]);exit(0);}
    for(int i=0;i<10;i++)dfs(x+6*ten[y]*i,y+1);
}
int main()
{
    //freopen("tt.in","r",stdin);
    scanf("%s",s+1);n=strlen(s+1);
    for(int i=1;i<=n;i++)v=v*10+s[i]-'0';
    ten[0]=1;
    for(int i=1;i<=20;i++)ten[i]=ten[i-1]*10;
    for(int i=0;i<60;i++)
    {
        dfs(i,1);
    }
    puts("NIE");
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值