A Simple Math Problem hdu 1757

30 篇文章 0 订阅
10 篇文章 0 订阅
#include <stdio.h>
#include <cstring>
#include <iostream>
using namespace std;
typedef long long LL;
LL tmp[11][11],a[11][11],b[11][11];
int mod;
void mul(LL a[][11],LL b[][11])
{
    for(int i=0; i<10; i++)
        for(int j=0; j<10; j++)
        {
            tmp[i][j]=0;
            for(int k=0; k<10; k++)
            {
                tmp[i][j]+=a[i][k]*b[k][j];
                if(tmp[i][j]>=mod) tmp[i][j]%=mod;
            }
        }
    memcpy(a,tmp,sizeof(tmp));
}
void pow(LL a[][11],LL b[][11],int n)
{
    memset(a,0,sizeof(a));
    for(int i=0; i<10; i++) a[i][i]=1;
    while(n)
    {
        if(n&1) mul(a,b);
        mul(b,b);
        n>>=1;
    }
}
int main()
{
    int k,x;
    while(scanf("%d%d",&k,&mod)==2)
    {
        memset(b,0,sizeof(b));
        memset(a,0,sizeof(a));
        for(int i=0; i<10; i++)
        {
            scanf("%d",&x);
            b[0][i]=x;
        }
        if(k<10)
        {
            printf("%d\n",k);
            continue;
        }
        for(int i=0; i<9; i++)
            b[i+1][i]=1;
        pow(a,b,k-9);
        int ans=0;
        for(int i=0; i<10; i++)
            ans+=a[0][i]*(9-i)%mod;
        printf("%d\n",ans%mod);
    }
    return 0;
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值