首师大附中科创教育平台 我的刷题记录 3120 LJX的校园：入学典礼

今天给大家献上“C”级题：LJX的校园：入学典礼！！

试题编号：3120

LJX的校园：入学典礼

难度级别：C； 运行时间限制：45ms； 运行空间限制：256000KB； 代码长度限制：2000000B

试题描述

LJX上小学啦！他与YSM,YSF,WHT,LTJ等人都是校友。今天，是他人生中“溺亡”的一天。今天，他要向同学们证明他的数学很“乐呵”。于是，刚学会简单的A+B问题的他，在课上，向冤家对头 斯沃琪 挑战 QAQ，斯沃琪 队有YZM,SJY,ZZQ等人。而LJX队有他的好朋（ji）友：YSM,YSF,WHT,LTJ,LZH等人，实力不弱小觑。有这么一道题：

给定正整数N，M，要求计算1,2，……,N连接起来（1234567891011……N）mod M的值。

“新兵”LJX想了想，要是M是3，或者9，他一定会。但是M什么都可以，只要小于INF。“预约”了各种方法以后，费劲脑筋想不出来。于是，右转向了（呵呵）他的（ZHU）队友们。可他们已经跑了 TAT。jue ruo LJX找到了你，他跪求你编一个程序帮他解决问题。否则，他将在毕业典礼上“锻炼”，并且被可怕的斯沃琪虐残 QAQ

输入

* 一行：正整数N，M。

输出

* 一行：按要求输出

输入示例

输入样例1
13 13
输入样例2
12345678910 1000000000

输出示例

输出样例1
4
输出样例2
345678910

其他说明

N<=10^18
M<=10^9
这数据是在坑LJX呀

用今天学的矩阵快速幂
哦，不，是分段矩阵快速幂

好的，以上就是LJX的校园：入学典礼的题目要求，现在献上代码！！！当当当！！！

#include<iostream>

#include<cstdio>

#include<cstdlib>

#include<cstring>

#include<algorithm>

#define MAXN 4

using namespace std;

typedef long long int LL;

int mod;

struct matrix

{

    LL p[MAXN][MAXN];

}ans,tmp;

matrix operator*(matrix a,matrix b)

{

    matrix c;

    for(int i=;i<=;i++)

        for(int j=;j<=;j++)

        {

            c.p[i][j]=;

            for(int k=;k<=;k++)

                c.p[i][j]=(c.p[i][j]+((a.p[i][k]%mod)*(b.p[k][j]%mod))%mod)%mod;

        }

    return c;

}

void cal(LL t,LL last) 

{

    memset(tmp.p,,sizeof(tmp.p));

    tmp.p[][]=t;

    tmp.p[][]=tmp.p[][]=tmp.p[][]=tmp.p[][]=tmp.p[][]=;

    LL y=last-t/+;

    while(y)

    {

        if(y&) ans=ans*tmp;

        tmp=tmp*tmp;

        y>>=;

    }

}

int main()

{

    for(int i=;i<=;i++)

        ans.p[i][i]=;

    LL n;

    scanf("%lld%lld",&n,&mod);

    LL t=;

    while(n>=t)

    {

        cal(t,t-);

        t*=;

    }

    cal(t,n);

    printf("%lld\n",ans.p[][]);

    return ;

}

LJX的校园：入学典礼！！！！！