UVA 11651 - Krypton Number System(DP+矩阵快速幂)

2024-06-01 19:58

本文主要是介绍UVA 11651 - Krypton Number System(DP+矩阵快速幂),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

UVA 11651 - Krypton Number System

题目链接

题意:给一个进制base,一个分数score求该进制下,有多少数满足一下条件:
1、没有连续数字
2、没有前导零
3、分数为score,分数的计算方式为相邻数字的平方差的和

思路:先从dp入手,dp[i][j]表示组成i,最后一个数字为j的种数,然后进行状态转移,推出前面一步能构成的状态,也就是到dp[(b - 1) * (b - 1)][x]。

然后可以发现后面的状态,都可以由前面这些状态统一转移出来,这样就可以利用矩阵快速幂进行优化,时间复杂度为O(n^3log(score))

这题矩阵快速幂姿势不够优美还会超时,矩阵相乘时候需要优化

代码:

#include <cstdio>
#include <cstring>typedef unsigned long long ll;
const int N = 155;
const ll MOD = (1ULL<<32);int n;struct Mat {ll v[N][N];Mat() {memset(v, 0, sizeof(v));}void init() {memset(v, 0, sizeof(v));}void init_one() {memset(v, 0, sizeof(v));for (int i = 0; i < n; i++)v[i][i] = 1;}Mat operator * (Mat c) {Mat ans;for (int k = 0; k < n; k++) {for (int i = 0; i < n; i++) {if (!v[i][k]) continue;for (int j = 0; j < n; j++)ans.v[i][j] = (ans.v[i][j] + v[i][k] * c.v[k][j] % MOD) % MOD;}}return ans;}
} B;Mat pow_mod(Mat x, int k) {Mat ans;ans.init_one();while (k) {if (k&1) ans = ans * x;x = x * x;k >>= 1;}return ans;
}ll dp[40][10], A[N];
int t, b, s, m;void init() {scanf("%d%d", &b, &s);m = (b - 1) * (b - 1);n = m  * b;memset(dp, 0, sizeof(dp));for (int i = 0; i < b; i++)dp[0][i] = 1;for (int i = 1; i <= m; i++) {for (int j = 0; j < b; j++) {for (int k = 0; k < b; k++) {if (k == j) continue;int tmp = (k - j) * (k - j);if (i - tmp < 0) continue;dp[i][j] = (dp[i][j] + dp[i - tmp][k]) % MOD;}}}
}void build() {B.init();for (int i = b; i < n; i++)B.v[i][i - b] = 1;for (int i = 1; i <= m; i++) {for (int j = 0; j < b; j++) {for (int k = 0; k < b; k++) {if (j == k) continue;if (i + (k - j) * (k - j) == m + 1) {B.v[(i - 1) * b + j][n - b + k] = 1;}}}}
}ll solve() {ll ans = 0;if (s <= m) {for (int i = 1; i < b; i++)ans = (ans + dp[s][i]) % MOD;return ans;}for (int i = 1; i <= m; i++) {for (int j = 0; j < b; j++)A[(i - 1) * b + j] = dp[i][j];}build();B = pow_mod(B, s - m);for (int i = n - b + 1; i < n; i++) {for (int j = 0; j < n; j++)ans = (ans + A[j] * B.v[j][i] % MOD) % MOD;}return ans;
}int main() {int cas = 0;scanf("%d", &t);while (t--) {init();printf("Case %d: %llu\n", ++cas, solve());}return 0;
}


这篇关于UVA 11651 - Krypton Number System(DP+矩阵快速幂)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1021976

相关文章

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s

C++快速排序超详细讲解

《C++快速排序超详细讲解》快速排序是一种高效的排序算法,通过分治法将数组划分为两部分,递归排序,直到整个数组有序,通过代码解析和示例,详细解释了快速排序的工作原理和实现过程,需要的朋友可以参考下... 目录一、快速排序原理二、快速排序标准代码三、代码解析四、使用while循环的快速排序1.代码代码1.由快

Win32下C++实现快速获取硬盘分区信息

《Win32下C++实现快速获取硬盘分区信息》这篇文章主要为大家详细介绍了Win32下C++如何实现快速获取硬盘分区信息,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 实现代码CDiskDriveUtils.h#pragma once #include <wtypesbase

Spring AI与DeepSeek实战一之快速打造智能对话应用

《SpringAI与DeepSeek实战一之快速打造智能对话应用》本文详细介绍了如何通过SpringAI框架集成DeepSeek大模型,实现普通对话和流式对话功能,步骤包括申请API-KEY、项目搭... 目录一、概述二、申请DeepSeek的API-KEY三、项目搭建3.1. 开发环境要求3.2. mav

Python如何快速下载依赖

《Python如何快速下载依赖》本文介绍了四种在Python中快速下载依赖的方法,包括使用国内镜像源、开启pip并发下载功能、使用pipreqs批量下载项目依赖以及使用conda管理依赖,通过这些方法... 目录python快速下载依赖1. 使用国内镜像源临时使用镜像源永久配置镜像源2. 使用 pip 的并

SpringBoot快速接入OpenAI大模型的方法(JDK8)

《SpringBoot快速接入OpenAI大模型的方法(JDK8)》本文介绍了如何使用AI4J快速接入OpenAI大模型,并展示了如何实现流式与非流式的输出,以及对函数调用的使用,AI4J支持JDK8... 目录使用AI4J快速接入OpenAI大模型介绍AI4J-github快速使用创建SpringBoot

使用Python快速实现链接转word文档

《使用Python快速实现链接转word文档》这篇文章主要为大家详细介绍了如何使用Python快速实现链接转word文档功能,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 演示代码展示from newspaper import Articlefrom docx import

Oracle数据库如何切换登录用户(system和sys)

《Oracle数据库如何切换登录用户(system和sys)》文章介绍了如何使用SQL*Plus工具登录Oracle数据库的system用户,包括打开登录入口、输入用户名和口令、以及切换到sys用户的... 目录打开登录入口登录system用户总结打开登录入口win+R打开运行对话框,输php入:sqlp

shell脚本快速检查192.168.1网段ip是否在用的方法

《shell脚本快速检查192.168.1网段ip是否在用的方法》该Shell脚本通过并发ping命令检查192.168.1网段中哪些IP地址正在使用,脚本定义了网络段、超时时间和并行扫描数量,并使用... 目录脚本:检查 192.168.1 网段 IP 是否在用脚本说明使用方法示例输出优化建议总结检查 1