矩阵快速幂,分1-9,10-99...看黄学长的代码理解。。。然而他直接把答案保存在最后一行(没有说明。。。好吧应该是我智障这都不知道。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define rep(i,n) for(int i=1;i<=n;i++)
#define clr(x,c) memset(x,c,sizeof(x))
#define ll long long
ll a[4][4],b[4][4],n,mod;
void mul(ll a[4][4],ll b[4][4],ll ans[4][4]){
ll tmp[4][4];clr(tmp,0);
rep(i,3) rep(j,3) rep(k,3) tmp[i][j]=(tmp[i][j]+a[i][k]*b[k][j])%mod;
rep(i,3) rep(j,3) ans[i][j]=tmp[i][j];
}
void cal(ll t,ll cnt){
clr(b,0);b[1][1]=t%mod,b[2][1]=b[2][2]=b[3][1]=b[3][2]=b[3][3]=1;
ll y=cnt-t/10+1;
while(y){
if(y&1) mul(a,b,a);
mul(b,b,b);y>>=1;
}
}
int main(){
scanf("%lld%lld",&n,&mod);
clr(a,0);a[3][3]=1;
ll t=10;
while(n>=t) cal(t,t-1),t*=10;
cal(t,n);printf("%lld\n",a[3][1]);
return 0;
}

  

2326: [HNOI2011]数学作业

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1711  Solved: 1018
[Submit][Status][Discuss]

Description

Input

 

Output

 

Sample Input

 

Sample Output

 

HINT

 

Source

 

[Submit][Status][Discuss]

最新文章

  1. 谈谈如何使用Netty开发实现高性能的RPC服务器
  2. 微信小程序-地图组件
  3. debuggap,移动端调试新方式
  4. tomcat 系统服务 outofmemory
  5. firebug使用
  6. [Redis源码阅读]redis持久化
  7. vue-cli 发布部署IIS
  8. 循环队列搜索 Search in Rotated Sorted Array
  9. nodejs+koa2微信app支付,小程序支付
  10. Android 使用 NYTimes Stores 缓存 network request
  11. @Html.Raw() 与Newtonsoft.Json.JsonConvert.SerializeObject()
  12. 人脸识别中的Procruster analysis应用
  13. Python3判断shell下进程是否存在&amp;&amp;启动&amp;&amp;邮件通知
  14. from setuptools import setup ImportError: No module named setuptools【转】
  15. LeetCode 1. Two Sum (JavaScript)
  16. NUMA架构
  17. R语言爬虫:爬取包含所有R包的名称及介绍
  18. Flutter入门之无状态组件
  19. RPi 2B QEMU 模拟树莓派
  20. nyoj 983 ——首尾相连数组的最大子数组和——————【最大子串和变形】

热门文章

  1. MVC的Ajax的异步请求
  2. [转载]C#时间函数
  3. CoreLocation MKMapView 地图
  4. MapInfo格式转arggis格式
  5. HDU 5629 Clarke and tree dp+prufer序列
  6. JNI中使用cl命令生成DLL文件
  7. 用hibernate自动创建mysql表,添加失败org.hibernate.exception.SQLGrammarException
  8. d3d11 effect state and default value tables
  9. HTML5 动画引擎 小记
  10. Javascript实现两张图片的延迟加载