bzoj 2326

矩阵乘法。。

这个矩阵的构造有点麻烦。。关键就从n转移到n-1（其实是我题做得少。。

分位数统计即可

#include<bits/stdc++.h>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return x*f;
}
struct mat{
    ll a[4][4];
}tmp,t;
ll n,m,k,inf;
mat operator*(const mat&x,const mat&y){
    mat s;mem(s.a);
    inc(i,1,3)
    inc(j,1,3)
    inc(k,1,3)s.a[i][j]=(x.a[i][k]*y.a[k][j]%inf+s.a[i][j]+inf)%inf;
    return s;
}
void out(mat x){
    inc(i,1,3){inc(j,1,3)printf("%lld ",x.a[i][j]);putchar('\n');
    }putchar('\n');
}
int main(){
    freopen("data.in","r",stdin);
    inc(i,1,3)tmp.a[i][i]=1;
    t.a[2][2]=t.a[1][3]=t.a[3][3]=t.a[3][2]=1;
    m=n=read();inf=read();m%=inf;
    for(k=1;k<=n&&k;k*=10){
        ll _t=min(k*10-1,n);
        t.a[1][1]=10*k%inf;ll v=_t-k+1;
//        out(t);
        for(mat s=t;v;v>>=1,s=s*s)if(v&1)tmp=s*tmp;
//        out(tmp);
    }
    printf("%lld\n",(tmp.a[1][3]+tmp.a[1][2]+inf)%inf);
    return 0;
}