题目大意:

一个由自然数组成的数列按下式定义:

对于i <= k:ai = bi

对于i > k: ai = c1ai-1 + c2ai-2 + ... + ckai-k

其中bj和 cj (1<=j<=k)是给定的自然数。写一个程序,给定自然数m <= n, 计算am + am+1 + am+2 + ... + an, 并输出它除以给定自然数p的余数的值。

题解

首先显然我们可以构造一个矩阵递推\(a_i\)。

如果直接从m递推到n,会超时(我也没写过,不知道),我们在矩阵中加一维,记录\(s_i\),具体可以见代码

注意一个问题:减法取模时应该加成正数。因为这个WA了好几次。

代码

#include <bits/stdc++.h>
#define ll long long
#define ull unsigned long long
using namespace std;
const ll maxn = 20;
ll k, B[maxn], C[maxn], s[maxn];
ll m, n, p;
struct M {
ll n, m;
ll a[maxn][maxn];
} a, b;
M operator*(M a, M b) {
M c;
c.n = a.n;
c.m = b.m;
memset(c.a, 0, sizeof(c.a));
for (ll i = 1; i <= c.n; i++) {
for (ll j = 1; j <= c.m; j++) {
for (ll k = 1; k <= a.m; k++)
c.a[i][j] = (c.a[i][j] + (ull)(a.a[i][k] * b.a[k][j]) % p) % p;
}
}
return c;
}
M pow(M a, ll b) {
M ret;
ret.n = a.n;
ret.m = a.m;
memset(ret.a, 0, sizeof(ret.a));
for (ll i = 1; i <= ret.n; i++)
ret.a[i][i] = 1;
while (b) {
if (b & 1)
ret = ret * a;
a = a * a;
b >>= 1;
}
return ret;
}
void print(M x) {
for (ll i = 1; i <= x.n; i++) {
for (ll j = 1; j <= x.m; j++)
cout << x.a[i][j] << ' ';
cout << endl;
}
}
ll calc(ll x) {
M y = pow(a, x + 1);
// prll (y);
y = y * b;
return y.a[1][1];
}
int main() {
// freopen("input", "r", stdin);
scanf("%lld", &k);
a.n = k + 1;
a.m = k + 1;
b.n = k + 1;
b.m = 1;
s[0] = 0;
for (ll i = 1; i <= k; i++) {
scanf("%lld", &B[i]);
}
for (ll i = 1; i <= k; i++)
scanf("%lld", &C[i]);
scanf("%lld %lld %lld", &m, &n, &p);
for (ll i = 1; i <= k; i++)
s[i] = (B[i] + s[i - 1]) % p;
b.a[1][1] = s[k - 1];
for (ll i = 1; i <= k; i++)
b.a[i + 1][1] = B[k - i + 1];
a.a[1][1] = a.a[1][2] = 1;
for (ll i = 1; i <= k; i++)
a.a[2][i + 1] = C[i];
for (ll i = 1; i < k; i++)
a.a[i + 2][i + 1] = 1;
// a = a * b;
ll ans1, ans2;
// calc(1);
if (m - 1 > k)
ans1 = calc(m - 1 - k);
else
ans1 = s[m - 1];
if (n > k)
ans2 = calc(n - k);
else
ans2 = s[n];
printf("%lld\n", (ans2 - ans1 + p) % p);
return 0;
}

[bzoj3231][SDOI2008]递归数列——矩阵乘法的更多相关文章

  1. bzoj 3231 [Sdoi2008]递归数列——矩阵乘法

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3231 矩阵乘法裸题. 1018是10^18.别忘了开long long. #include& ...

  2. 【bzoj3231】[Sdoi2008]递归数列 矩阵乘法+快速幂

    题目描述 一个由自然数组成的数列按下式定义: 对于i <= k:ai = bi 对于i > k: ai = c1ai-1 + c2ai-2 + ... + ckai-k 其中bj和 cj  ...

  3. [luogu2461 SDOI2008] 递归数列 (矩阵乘法)

    传送门 Description 一个由自然数组成的数列按下式定义: 对于i <= k:ai = bi 对于i > k: ai = c1ai-1 + c2ai-2 + ... + ckai- ...

  4. P2461 [SDOI2008]递归数列 矩阵乘法+构造

    还好$QwQ$ 思路:矩阵快速幂 提交:1次 题解: 如图: 注意$n,m$如果小于$k$就不要快速幂了,直接算就行... #include<cstdio> #include<ios ...

  5. BZOJ3231: [Sdoi2008]递归数列

    BZOJ3231: [Sdoi2008]递归数列 Description 一个由自然数组成的数列按下式定义: 对于i <= k:ai = bi 对于i > k: ai = c1ai-1 + ...

  6. BZOJ 3231: [Sdoi2008]递归数列( 矩阵快速幂 )

    矩阵乘法裸题..差分一下然后用矩阵乘法+快速幂就可以了. ----------------------------------------------------------------------- ...

  7. bzoj 3231 [ Sdoi 2008 ] 递归数列 —— 矩阵乘法

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3231 裸矩阵乘法. 代码如下: #include<iostream> #incl ...

  8. BZOJ-3231 [SDOI2008]递归数列

    转成矩阵连乘后,矩阵快速幂加速解决. 一开始没把需要longlong的变量补全..而且没初始化2333 #include <cstdlib> #include <cstdio> ...

  9. BZOJ_3231_[Sdoi2008]递归数列_矩阵乘法

    BZOJ_3231_[Sdoi2008]递归数列_矩阵乘法 Description 一个由自然数组成的数列按下式定义: 对于i <= k:ai = bi 对于i > k: ai = c1a ...

随机推荐

  1. 实现一个简单版的express

    express应该算是早期最优秀的一个node框架了,刚开始学node做后端语言就是用的express,它的cli可以帮我们搭建好项目目录,就像现在的vue,react一样.express本身没有做太 ...

  2. 使用JDK自带的keytool工具生成证书

    一.keytool 简介 keytool 是java用于管理密钥和证书的工具,它使用户能够管理自己的公钥/私钥对及相关证书,用于(通过数字签名)自我认证(用户向别的用户/服务认证自己)或数据完整性以及 ...

  3. 51单片机数码管字符H自右向左移动

    #include <reg51.h> #define uint unsigned int #define uchar unsigned char sfr P0M0 = 0x94; sfr ...

  4. Laxcus大数据管理系统2.0(3)- 第一章 基础概述 1.2 产品特点

    1.2 产品特点 Laxcus大数据管理系统运行在计算机集群上,特别强调软件对分布资源可随机增减的适应性.这种运行过程中数据动态波动和需要瞬时感知的特点,完全不同与传统的集中处理模式.这个特性衍生出一 ...

  5. Sersync实时备份服务部署实践

  6. 0.爬虫 urlib库讲解 urlopen()与Request()

    # 注意一下 是import urllib.request 还是 form urllib import request 0. urlopen() 语法:urllib.request.urlopen(u ...

  7. 无法注册DLL/OCX:regsvr32 failed with exit code 0x5

    最近在服务器上装一个等值线控件,确报错: RegSvr32 failed with exit code 0x5   尝试解决方法1: 在cmd下 运行for %1 in (%windir%\syste ...

  8. URAL 1736 Chinese Hockey(网络最大流)

    Description Sergey and Denis closely followed the Chinese Football Championship, which has just come ...

  9. NO1——线段树

    /* 数组存储 */ /* 预处理 */ #include <iostream> #include <cstdio> #include <algorithm> #i ...

  10. sqoop将oracle数据导入hdfs集群

    使用sqoop将oracle数据导入hdfs集群 集群环境: hadoop1.0.0 hbase0.92.1 zookeeper3.4.3 hive0.8.1 sqoop-1.4.1-incubati ...