【矩阵】[HNOI2011][HYSBZ/BZOJ2326][CQBZOJ2831]数学作业

本文介绍了一种使用矩阵快速幂的方法来高效计算从1到N的所有整数连接起来形成的数字串对M取模的问题。通过定义转移矩阵并利用快速幂算法,可以有效地解决该问题,尤其是在N较大时。

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

题目
分析:另f[n]为Concatenate(1..N) Mod M的值,那么f[n]=(f[n-1]*10^length[n]+n)%MOD。
所以

[f[n]n1]=[f[n1]n1]111011001

根据数字的长度分段做矩阵快速幂即可。

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define MAXMT 3
#define MOD m
typedef int matrix[MAXMT+1][MAXMT+1];
matrix a,b;
long long n;
int m,ans;
template<class T>
void Read(T &x){
    char c;
    while(c=getchar(),c!=EOF)
        if(c>='0'&&c<='9'){
            x=c-'0';
            while(c=getchar(),c>='0'&&c<='9')
                x=x*10+c-'0';
            ungetc(c,stdin);
            return;
        }
}
void mul(matrix a,matrix b,matrix c){
    matrix d;
    memset(d,0,sizeof d);
    int i,j,k;
    for(i=1;i<=MAXMT;i++)
        for(j=1;j<=MAXMT;j++)
            for(k=1;k<=MAXMT;k++)
                d[i][j]=(d[i][j]+1ll*a[i][k]*b[k][j])%MOD;
    memcpy(c,d,sizeof d);
}
void quick_pow(matrix a,long long b,matrix c){
    matrix d,t;
    memcpy(t,a,sizeof t);
    memset(d,0,sizeof d);
    for(int i=1;i<=MAXMT;i++)
        d[i][i]=1;
    while(b){
        if(b&1)
            mul(d,t,d);
        mul(t,t,t);
        b>>=1;
    }
    memcpy(c,d,sizeof d);
}
void solve(){
    a[1][1]=a[2][1]=a[2][2]=a[3][1]=a[3][2]=a[3][3]=1;
    long long t=1;
    while(n>=t){
        a[1][1]=1ll*a[1][1]*10%MOD;
        quick_pow(a,min(n,t*10-1)-t+1,b);
        ans=(1ll*ans*b[1][1]%MOD+(t-1)%MOD*b[2][1]%MOD+b[3][1])%MOD;
        t*=10;
    }
}
int main()
{
    Read(n),Read(m);
    solve();
    printf("%d\n",ans);
}
根据引用\[1\]和引用\[2\]的描述,题目中的影魔拥有n个灵魂,每个灵魂有一个战斗力ki。对于任意一对灵魂对i,j (i<j),如果不存在ks (i<s<j)大于ki或者kj,则会为影魔提供p1的攻击力。另一种情况是,如果存在一个位置k,满足ki<c<kj或者kj<c<ki,则会为影魔提供p2的攻击力。其他情况下的灵魂对不会为影魔提供攻击力。 根据引用\[3\]的描述,我们可以从左到右进行枚举。对于情况1,当扫到r\[i\]时,更新l\[i\]的贡献。对于情况2.1,当扫到l\[i\]时,更新区间\[i+1,r\[i\]-1\]的贡献。对于情况2.2,当扫到r\[i\]时,更新区间\[l\[i\]+1,i-1\]的贡献。 因此,对于给定的区间\[l,r\],我们可以根据上述方法计算出区间内所有下标二元组i,j (l<=i<j<=r)的贡献之和。 #### 引用[.reference_title] - *1* *3* [P3722 [AH2017/HNOI2017]影魔(树状数组)](https://blog.youkuaiyun.com/li_wen_zhuo/article/details/115446022)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [洛谷3722 AH2017/HNOI2017 影魔 线段树 单调栈](https://blog.youkuaiyun.com/forever_shi/article/details/119649910)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值