【bzoj4002】[JLOI2015]有意义的字符串

最新推荐文章于 2024-06-23 15:17:45 发布

愤怒的愣头青

最新推荐文章于 2024-06-23 15:17:45 发布

阅读量366

点赞数 1

分类专栏：矩阵乘法

本文链接：https://blog.csdn.net/w_yqts/article/details/78916458

版权

矩阵乘法专栏收录该内容

11 篇文章 0 订阅

订阅专栏

Description

B 君有两个好朋友，他们叫宁宁和冉冉。

有一天，冉冉遇到了一个有趣的题目：输入 b;d;n，求((b+sqrt(D)/2)^N的整数部分，请输出结果 Mod 7528443412579576937 之后的结果吧。
Input

一行三个整数 b;d;n

Output

一行一个数表示模 7528443412579576937 之后的结果。

Sample Input

1 5 9
Sample Output

76
HINT

#include<bits/stdc++.h>
#define ll unsigned long long
#define inf 1000000000
using namespace std;
#define mod 7528443412579576937UL
const int N=30005;
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;
}
ll ans[3][3],a[3][3],b,d,n;
ll cal(ll a,ll b)
{
    ll ans=0;
    while (b)
    {
        if (b&1) ans=(ans+a)%mod;
        b>>=1;
        a=(a+a)%mod;
    }
    return ans;
}
void mul(ll a[3][3],ll b[3][3],ll c[3][3])
{
    ll tmp[3][3];
    for (int i=1;i<=2;i++)
        for (int j=1;j<=2;j++)
        {
            tmp[i][j]=0LL;
            for (int k=1;k<=2;k++)
                tmp[i][j]=(tmp[i][j]+cal(a[i][k],b[k][j])%mod)%mod;
        }
    for (int i=1;i<=2;i++)
        for (int j=1;j<=2;j++)
            c[i][j]=tmp[i][j];
}
void MUL(ll n)
{
    while (n)
    {
        if (n&1) mul(ans,a,ans);
        mul(a,a,a);
        n/=2;
    }
}
int main()
{
    scanf("%lld%lld%lld",&b,&d,&n);if (n==0) return 0,puts("1");
    a[1][1]=b;a[1][2]=(d-cal(b,b))/4;a[2][1]=1;
    ans[1][1]=ans[2][2]=1;MUL(n-1);
    ll A=(cal(ans[1][1],b)%mod+cal(ans[1][2],2)%mod)%mod;
    if (n%2==0&&cal(b,b)!=d) A--;A=(A+mod)%mod;
    cout<<A;
    return 0;
}