和斐波那契一个道理在最后加一个求和即可
 #include<cstdio>
#include<cstring>
#include<iostream>
//using namespace std;
const int maxn=;
const double eps=1e-;
long long modn;
long long n,l,r;
long long b[]={};
struct mat{
long long e[][];
mat(){ memset(e,,sizeof(e)); }
};
mat a;
mat Mul(mat x,mat y){
mat z;
for(int i=;i<=n+;i++){
for(int j=;j<=n+;j++){
for(int k=;k<=n+;k++){
z.e[i][j]+=x.e[i][k]*y.e[k][j];
z.e[i][j]%=modn;
}
}
}
return z;
}
mat Pow(mat x,long long k){
mat z;
for(int i=;i<=n+;i++){
z.e[i][i]=;
}
while(k>){
if(k&){
z=Mul(z,x);
}
x=Mul(x,x);
k/=;
}/*for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
std::cout<<z.e[i][j]<<' ';
}
std::cout<<std::endl;
}*/
return z;
}
long long doit(long long x){
if(x<n){
return b[x+];
}
mat z=Pow(a,x-n+);
long long ans=,s=,d=;
for(int i=;i<=n+;i++){
d+=z.e[n][i]*b[i];
s+=z.e[n+][i]*b[i];
d%=modn;s%=modn;
}
ans=(s+d)%modn;
ans%=modn;
return ans;
}
int main(){
scanf("%lld",&n);
n+=;
for(int i=;i<=n;i++){
scanf("%lld",&b[i]);
b[n+]+=b[i];
}
b[n+]-=b[n];
for(int i=n;i>;i--){
scanf("%lld",&a.e[n][i]);
}
long long l,r;
scanf("%lld%lld%lld",&l,&r,&modn);
for(int i=;i<=n;i++){
b[i]%=modn;
a.e[i-][i]=;a.e[n][i]%=modn;
}
a.e[n+][n+]=,a.e[n+][n]=;
/*for(int i=1;i<=n+1;i++){
for(int j=1;j<=n+1;j++){
std::cout<<a.e[i][j]<<' ';
}
std::cout<<std::endl;
}*/
long long ans=(doit(r)-doit(l-)+modn)%modn;
printf("%lld\n",ans);
return ;
}

最新文章

  1. c#变量缺少using引用,如何快速加上using,加Using的快捷键[bubuko.com]
  2. linux环境下配置solr5.3详细步骤
  3. python file 文件读写
  4. json:There is a cycle in the hierarchy!
  5. hdu 1020 Encoding
  6. MySQL数据库指定字符集
  7. 拾遗:『Linux Capability』
  8. Spark中加载本地(或者hdfs)文件以及SparkContext实例的textFile使用
  9. 九度OJ1207
  10. MSP430常见问题之指令系统类
  11. js多行字符串
  12. requirejs 合并方案
  13. BF算法(朴素的模式匹配算法)
  14. 深度分析DataTable如何筛选
  15. C#中Invoke的用法
  16. HDU3311Dig The Wells
  17. C#中添加log4net(日志文件)
  18. Nuget私有服务搭建实战
  19. emWin 界面切换注意事项
  20. 如何从angular2中的url获取查询参数?

热门文章

  1. 基本控件文档-UITextField属性---iOS-Apple苹果官方文档翻译
  2. windos下创建软链接,附Linux下创建软链接
  3. Java 注解全面解析
  4. c++刷题(6/100)最长上升子序列
  5. 一款已上市MMO手游地图同步方案总结
  6. Attention is all you need 论文详解(转)
  7. Machine Learning系列--隐马尔可夫模型的三大问题及求解方法
  8. Springboot swagger2 导出api文档
  9. Spring之IOC,DI,动态代理,反射
  10. [ python ] 购物系统