数学-HDU 4565 so easy

最新推荐文章于 2021-03-24 22:50:27 发布

跑调的叶子

最新推荐文章于 2021-03-24 22:50:27 发布

阅读量815

点赞数

分类专栏：题和早期用的模板文章标签：数学

本文链接：https://blog.csdn.net/sinianluoye/article/details/24196873

版权

题和早期用的模板专栏收录该内容

29 篇文章 0 订阅

订阅专栏

一道数学味很浓的题，如果不是看了题解真忘了那些数学知识了····

这是解析

http://blog.csdn.net/crazy______/article/details/9021169

看明白了直接矩阵快速幂就可以了（原谅我英语不好···）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include <algorithm>
#define ll long long
#define eps 1e-8
#define ms(x,y) (memset(x,y,sizeof(x)))
#define fr(i,x,y) for(ll i=x;i<=y;i++)
using namespace std;
struct juzhen
{
    ll p[3][3];
    juzhen(){};
};
juzhen mul(juzhen a,juzhen b)
{
    juzhen c;
    c.p[1][1]=a.p[1][1]*b.p[1][1]+a.p[1][2]*b.p[2][1];
    c.p[1][2]=a.p[1][1]*b.p[1][2]+a.p[1][2]*b.p[2][2];
    c.p[2][1]=a.p[2][1]*b.p[1][1]+a.p[2][2]*b.p[2][1];
    c.p[2][2]=a.p[2][1]*b.p[1][2]+a.p[2][2]*b.p[2][2];
    return c;
}
juzhen fast_mod(juzhen a,ll n,ll m)
{
    juzhen ret;
    ret.p[1][1]=ret.p[2][2]=1;
    ret.p[1][2]=ret.p[2][1]=0;
    while(n)
    {
        if(n&1)
        {
            ret=mul(ret,a);
            fr(i,1,2)fr(j,1,2)
            {
                ret.p[i][j]%=m;
                if(ret.p[i][j]<0)ret.p[i][j]+=m;
            }
        }
        a=mul(a,a);
        fr(i,1,2)fr(j,1,2)
        {
            a.p[i][j]%=m;
            if(a.p[i][j]<0)a.p[i][j]+=m;
        }
        n/=2;
    }
    return ret;
}
int main()
{
    ll a,b,n,m;
    while(~scanf("%I64d%I64d%I64d%I64d",&a,&b,&n,&m))
    {
        juzhen ans;
        ans.p[1][1]=2*a;
        ans.p[1][2]=-(a*a-b);
        ans.p[2][1]=1;
        ans.p[2][2]=0;
        ans=fast_mod(ans,n,m);
        ll t=ans.p[2][1]*2*a+2*ans.p[2][2];
        cout<<t%m<<endl;
    }
}