「 洛谷 」P2768 珍珠项链
珍珠项链
题目限制
- 内存限制:125.00MB
- 时间限制:1.00s
- 标准输入输出
题目知识点
- 动态规划 \(dp\)
- 矩阵
- 矩阵乘法
- 矩阵加速
- 矩阵快速幂
题目来源
为了方便大家阅读通畅,题目可能略有改动,保证不会造成影响
题目
题目背景
小\(L\) 通过泥萌的帮助,成功解决了牛栏的修建问题
奶牛们觉得主人非常厉害,就再也不偷懒了:母牛们奋力挤奶、生娃
小\(L\) 也因此成为了一代富豪
但是一直困扰 小\(L\) 的就是单身问题
小\(L\) 经过长久的寻觅,小\(L\) 终于找到了一个心仪的漂亮妹子
于是,小\(L\) 打算在 \(5.20\) 那天给妹子一个惊喜
题目描述
小\(L\) 决定用 \(K\) 种珍珠为妹子做一串举世无双的珍珠垂饰
珍珠垂饰是由珍珠连接而成的,其长度可以认为就是珍珠垂饰上珍珠的个数
小\(L\) 现在腰缠万贯,每种珍珠他都拥有 \(N\) 颗
根据将珍珠垂饰打开后珍珠不同的排列顺序可以区别出不同种类的项链
现在,小\(L\) 好奇自己可以组成多少种长度为 \(1\) 至 \(N\) 且不同的珍珠垂饰
当然,为显富有,每串珍珠垂饰都要必须由 \(K\) 种珍珠连接而成
由于答案可能太大,只需要取模 1234567891
这一定难不倒聪明的你吧
如果你能帮 小\(L\) 解决这个问题,他会把最后的资产分给你一半哦
格式
输入格式
第一行输入一个整数 \(T\) ,表示测试数据的个数
- 对于每组数据:
- 输入占一行,包含两个整数 \(N\) 和 \(K\),用一个空格隔开
输出格式
- 对于每组数据
- 输出占一行,包含一个整数,表示项链的种类数
样例
样例输入
2
2 1
3 2
样例输出
2
8
提示
时间限制
对于 \(70\% \sim 100\%\) 的数据:时间限制 \(10ms\)
数据范围
对于 \(40\%\) 的数据:\(1<= N<= 1e5\),\(0 \leq K \leq 30\)
对于 \(100\%\) 的数据:\(T \leq 10\),\(1 \leq N \leq 1e9\),\(0 \leq K \leq 30\)
思路
题目所说的 珍珠项链长度为 \(i\) 的种类数 其实就是 用不同的方式连接出长度为 \(i\) 的珍珠项链方案数
我们可以先思考 \(dp\):\(dp[i][j]\) 表示 用 \(j\) 种珍珠 连接出 长度为 \(i\) 的珍珠项链的方案数
- 很多读者在这里可能会有一个疑问:\(dp\) 的状态 \(dp[i][j]\),似乎没有考虑每种珍珠的数量
- 因为每种珍珠的个数是 \(N\) 颗,所需要求的最大长度也是 \(N\);也就是说,可以只用 \(1\) 种珍珠连接出所有的项链,每种珍珠的数量都是足够的;因此,珍珠的数量是不会造成影响的,可以忽略不计
我们可以进一步推出 \(dp\) 的状态转移方程:\(dp[i][j] = dp[i - 1][j - 1] * (K - (j - 1)) + dp[i - 1][j] * j\)
- \(dp[i - 1][j - 1] * (k - (j - 1))\):当 第 \(i\) 颗珍珠 与 前面 \(i - 1\) 颗用过的珍珠种类不同 的时候,所连接出长度为 \(i\) 的珍珠的方案数
- \(dp[i - 1][j - 1]\):用 \(j - 1\) 种珍珠 连接出 长度为 \(i - 1\) 的项链 的方案数
- \(k - (j - 1)\):前面的 \(i - 1\) 颗珍珠 已经用了 \(j - 1\) 种珍珠,而 第 \(i\) 颗使用没有用过的珍珠,就会有 剩下的珍珠种数 种情况,即 \(k - (j - 1)\)
- \(dp[i - 1][j] * j\):当 第 \(i\) 颗的珍珠的种类 是之前用过的时候,所连接出来的方案数
- \(dp[i - 1][j]\):用 \(j\) 种珍珠 连接出 长度为 \(i - 1\) 的项链 的方案数
- \(j\):前面 \(i - 1\) 颗珍珠用了 \(j\) 种,第 \(i\) 颗珍珠 还是使用 这 \(j\) 种中的某一种,就有 \(j\) 种情况
最后的答案就是 \(dp[1][K] + dp[2][K] + dp[3][K] + ··· + dp[N][K] = \sum_{i = 1}^{N} dp[i][K]\)
分析
如果这道题的数据较小的话,用 \(dp\) 就可以 \(AC\) 了,可惜 \(N\) 的范围太大了
这时候我们就要思考如何优化 \(dp\) 的时间复杂度了
\(dp\) 是在进行递推转移,不妨可以把 \(dp\) 放在矩阵里求解
我们设 \(ans[i] = dp[1][K] + dp[2][K] + ··· + dp[i][K]\),我们所求的答案就是 \(ans[N]\)
我们可以构造一个 原始矩阵:
& dp[1][1] & dp[1][2] & dp[1][3] & ··· & dp[1][K] & ans[0] & \\
\end{bmatrix}
\]
为了状态转移,我们需要将它变成:
& dp[2][1] & dp[2][2] & dp[2][3] & ··· & dp[2][K] & ans[1] & \\
\end{bmatrix}
\]
继续进化成:
& dp[3][1] & dp[3][2] & dp[3][3] & ··· & dp[3][K] & ans[2] & \\
\end{bmatrix}
\]
以此类推,直到求出 \(ans[N]\) 为止
从 构造的原始矩阵 和 \(dp\) 的状态转移方程 中,我们可以推出 加速矩阵 (\(k\) 行 \(k\) 列):
& 1 & k - 1 & 0 &0 & ··· & 0 & 0 & \\
& 0 & 2 & k - 2 &0 & ··· & 0 & 0 & \\
& 0 & 0 & 3 & k - 3 & ··· & 0 & 0 & \\
& ··· & ··· & ··· & ··· & ··· & 0 & 0 & \\
& ··· & ··· & ··· & ··· & ··· & 1 & 0 & \\
& 0 &0 &0 &0 &0 &k &1 & \\
& 0 &0 &0 &0 &0 &0 &1 & \\
\end{bmatrix}
\]
刚开始是 \(ans[0]\),答案是 \(ans[N]\),所以要乘以 \(N\) 个加速矩阵,可是仍然要超时
这时候就需要矩阵快速幂了:根据矩阵乘法的结合律,先把 \(N\) 个加速矩阵乘起来,再用 原始矩阵 乘以 这个得到的矩阵,就可以得到最终的答案了
(注意:矩阵1 \(*\) 矩阵2 不一定等于 矩阵2 \(*\) 矩阵1,所以不能乘反了)
代码
#include <cstdio>
#include <cstring>
const long long MOD = 1234567891LL;
int T, N, K;
const int MAXK = 30;
struct Matrix
{
long long Mat[MAXK + 5][MAXK + 5];
int R, C;
Matrix()
{
for (int i = 1; i <= MAXK + 1; i++)
for (int j = 1; j <= MAXK + 1; j++)
Mat[i][j] = 0LL;
}
void Read() const // 矩阵的输入
{
if (R < 1 || C < 1) return ;
for (int i = 1; i <= R; i++)
for (int j = 1; j <= C; j++)
scanf("%lld", &Mat[i][j]);
}
void Write() const // 矩阵的输出
{
if (R < 1 || C < 1) return ;
for (int i = 1; i <= R; i++) {
for (int j = 1; j < C; i++)
printf("%lld ", Mat[i][j]);
printf("%lld\n", Mat[i][C]);
}
}
Matrix operator * (const Matrix One) const // 重载矩阵的乘号
{
Matrix Res;
Res.R = R, Res.C = C;
for (int i = 1; i <= Res.R; i++)
for (int j = 1; j <= Res.C; j++)
for (int k = 1; k <= One.R; k++)
Res.Mat[i][j] = (Res.Mat[i][j] + Mat[i][k] * One.Mat[k][j]) % MOD;
return Res;
}
}A, B;
Matrix Pow(Matrix One, long long k) // 矩阵快速幂
{
Matrix Res, cnt = One;
Res.R = K + 1, Res.C = K + 1;
for (int i = 1; i <= K + 1; i++)
Res.Mat[i][i] = 1LL; // 单位矩阵
for (int i = k; i >= 1; i >>= 1)
{
if (i & 1) Res = Res * cnt;
cnt = cnt * cnt;
}
return Res;
}
void Init(int k) // 初始化 原始矩阵 和 加速矩阵
{
A.R = 1, A.C = k + 1;
A.Mat[1][1] = (long long)k;
B.R = B.C = k + 1;
B.Mat[1][1] = 1LL;
for (int i = 2; i <= k; i++)
{
B.Mat[i][i] = (long long)i;
B.Mat[i - 1][i] = (long long)k + 1LL - i;
}
B.Mat[k][k + 1] = B.Mat[k + 1][k + 1] = 1LL;
}
int main()
{
scanf("%d", &T);
while (T--)
{
scanf("%d %d", &N, &K);
Init(K);
Matrix ans = A * Pow(B, N);
printf("%lld\n", ans.Mat[1][K + 1]);
}
return 0;
}
「 洛谷 」P2768 珍珠项链的更多相关文章
- 「 洛谷 」P4539 [SCOI2006]zh_tree
小兔的话 推荐 小兔的CSDN [SCOI2006]zh_tree 题目限制 内存限制:250.00MB 时间限制:1.00s 标准输入输出 题目知识点 思维 动态规划 \(dp\) 区间\(dp\) ...
- 「 洛谷 」P2151 [SDOI2009]HH去散步
小兔的话 欢迎大家在评论区留言哦~ HH去散步 题目限制 内存限制:125.00MB 时间限制:1.00s 标准输入 标准输出 题目知识点 动态规划 \(dp\) 矩阵 矩阵乘法 矩阵加速 矩阵快速幂 ...
- 「区间DP」「洛谷P1043」数字游戏
「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...
- 「P4994」「洛谷11月月赛」 终于结束的起点(枚举
题目背景 终于结束的起点终于写下句点终于我们告别终于我们又回到原点…… 一个个 OIer 的竞赛生涯总是从一场 NOIp 开始,大多也在一场 NOIp 中结束,好似一次次轮回在不断上演.如果这次 NO ...
- 「洛谷4197」「BZOJ3545」peak【线段树合并】
题目链接 [洛谷] [BZOJ]没有权限号嘤嘤嘤.题号:3545 题解 窝不会克鲁斯卡尔重构树怎么办??? 可以离线乱搞. 我们将所有的操作全都存下来. 为了解决小于等于\(x\)的操作,那么我们按照 ...
- 「洛谷3338」「ZJOI2014」力【FFT】
题目链接 [BZOJ] [洛谷] 题解 首先我们需要对这个式子进行化简,否则对着这么大一坨东西只能暴力... \[F_i=\sum_{j<i} \frac{q_iq_j}{(i-j)^2}-\s ...
- 「BZOJ2733」「洛谷3224」「HNOI2012」永无乡【线段树合并】
题目链接 [洛谷] 题解 很明显是要用线段树合并的. 对于当前的每一个连通块都建立一个权值线段树. 权值线段树处理操作中的\(k\)大的问题. 如果需要合并,那么就线段树暴力合并,时间复杂度是\(nl ...
- 「洛谷3870」「TJOI2009」开关【线段树】
题目链接 [洛谷] 题解 来做一下水题来掩饰ZJOI2019考炸的心情QwQ. 很明显可以线段树. 维护两个值,\(Lazy\)懒标记表示当前区间是否需要翻转,\(s\)表示区间还有多少灯是亮着的. ...
- 「洛谷5300」「GXOI/GZOI2019」与或和【单调栈+二进制转化】
题目链接 [洛谷传送门] 题解 按位处理. 把每一位对应的图都处理出来 然后单调栈处理一下就好了. \(and\)操作处理全\(1\). \(or\)操作处理全\(0\). 代码 #include & ...
随机推荐
- 2w+长文带你剖析ConcurrentHashMap~!
并发编程实践中,ConcurrentHashMap是一个经常被使用的数据结构,相比于Hashtable以及Collections.synchronizedMap(),ConcurrentHashMap ...
- HttpClient4.5X使用-集成微服务
HttpClient4.5X使用-集成微服务 1.什么是HttpClient HTTP 协议可能是现在 Internet 上使用得最多.最重要的协议了,越来越多的 Java 应用程序需要直 ...
- http://www.etymon.cn/yingyucigen/3093.html
import requests import lxml.etree as etree import xml.etree.ElementTree as ET # 详情页 # 3093-148 # htt ...
- 「NOIP2009」最优贸易 题解
「NOIP2009」最优贸易 题解 题目TP门 题目描述 \(C\)国有\(n\)个大城市和\(m\)条道路,每条道路连接这\(n\)个城市中的某两个城市.任意两个城市之间最多只有一条道路直接相连.这 ...
- [Kafka][1][初识Kafka]
目录 第1章 初识Kafka 1.1 发布与订阅消息系统 1.1.1 如何开始 1.1.2 独立的队列系统 1.2 Kafka登场 1.2.1 消息和批次(Message and batch) 1.2 ...
- 厉害!这份阿里面试官 甩出的Spring源码笔记,GitHub上已经爆火
前言 时至今日,Spring 在 Java 生态系统与就业市场上,面试出镜率之高,投产规模之广,无出其右.随着技术的发展,Spring 从往日的 IoC 框架,已发展成 Cloud Native 基础 ...
- Camtasia中对录制视频进行编辑——旁白
相信很多人都遇见过想要录制视频,但是不知道在电脑上用哪一款软件比较好,害怕自己录的视频导出来之后会有水印,或者在录制的过程中遇到麻烦,更或者下载一款带有病毒的软件.那么今天我便给大家推荐一款专业录制屏 ...
- 如何卸载MathType 7?
作为好用的公式编辑器,一般情况下是不会将其从电脑上卸载的,但是当电脑负荷过多,导致电脑运行缓慢时,就需要考虑卸载一些软件,本节就来学习卸载MathType 7的方法. 具体操作步骤如下: 1.打开控制 ...
- 如何修改Docker的默认网段?
背景 公司的服务器全部都是购买的阿里云ECS主机,默认内网网段为172.16.0.0/12,测试服上的服务基本都是使用Docker部署的. 之前的同事在该机器上安装的Docker全部都是默认配置,网段 ...
- 走进flex布局
简介:flex 是一个CSS的display 属性中新添加一个值. 随着inline-flex的使用,它将使它适用的元素成为一个flex container(伸缩容器),而这个元素的每个子元素将成为 ...