HDU 2276 Kiki & Little Kiki 2( 矩阵快速幂 + 循环同构矩阵 )
蒟蒻的我还需深入学习
**链接:****传送门 **
题意:给出一个长度为 n,n 不超过100的 01 串 s ,每当一个数字左侧为 1 时( 0的左侧是 n-1 ),这个数字就会发生改变,整个串改变一次需要 1s ,询问 M s 后此串变为什么样子,例如 0101111 ,1s 后变为 1111000
思路:
根据题意可以得到这样一个规律 s[ i ] = ( s[ i - 1 ] + s[ i ] ) % 2,特别的 s[ 0 ] = ( s[ n-1 ] + s[ 0 ] ) % 2 ( s[ ] 不再考虑为char ),构造矩阵a、b、ans
- 假设 01 串长为 5,下图为矩阵a
1 0 0 0 1 1 1 0 0 0 0 1 1 0 0 0 0 1 1 0 0 0 0 1 1 - 矩阵b
s0 0 0 0 0 s1 0 0 0 0 s2 0 0 0 0 s3 0 0 0 0 s4 0 0 0 0 - 矩阵 ans
s0 0 0 0 0 s1 0 0 0 0 s2 0 0 0 0 s3 0 0 0 0 s4 0 0 0 0 - 矩阵 ans = a * b,假设经过 Ms 后矩阵 ans = pow( a, M ) * b,此时的 b 中的 s[ ] 为起始的 01 串,所以通过矩阵快速幂就可以解决了
解决此题的三种姿势
1.普通矩阵快速幂,经过上面分析很容易写出该方法对应的程序
2.循环同构矩阵优化矩阵快速幂,经上面分析可以看出矩阵 a 为典型的循环同构矩阵,( 什么是循环矩阵?**戳这里! ** )对于循环同构矩阵可以先算出来第 0 行然后递推下面的剩余行,从而将计算 pow( a , n ) 的复杂度 O( n ^ 3 ) 降为 O( n ^ 2 ),需要注意的时,这种方法应该只能用来优化同构矩阵的计算,如果是一个同构矩阵 × 其他矩阵,递推剩下行数是错误的!这道题不会卡复杂度......
3.循环同构矩阵优化矩阵快速幂 + 位运算,可以把矩阵乘法中的 * 换成 & , % 2 换成 ^(XOR)
姿势1:
/*************************************************************************
> File Name: hdu2276.cpp
> Author: WArobot
> Blog: http://www.cnblogs.com/WArobot/
> Created Time: 2017年05月04日 星期四 14时32分45秒
************************************************************************/
#include<bits/stdc++.h>
using namespace std;
int n,len;
string s;
const int MOD = 2;
const int maxn = 110;
#define ll long long
#define mod(x) ((x)%MOD)
#define cls(x) memset(x,0,sizeof(x));
struct mat{
int m[maxn][maxn];
}unit;
void init_unit(){
for(int i=0;i<maxn;i++) unit.m[i][i] = 1;
}
mat operator *(mat a,mat b){
mat ret;
ll x;
for(int i=0;i<len;i++){
for(int j=0;j<len;j++){
x = 0;
for(int k=0;k<len;k++)
x += mod( a.m[i][k]*b.m[k][j] );
ret.m[i][j] = mod(x);
}
}
return ret;
}
mat pow_mat(mat a,int x){
mat ret = unit;
while(x){
if(x&1) ret = ret*a;
a = a*a;
x >>= 1;
}
return ret;
}
mat a,b;
void init(){
cls(a.m);
a.m[0][0] = a.m[0][len-1] = 1;
for(int i=1;i<len;i++) a.m[i][i-1] = a.m[i][i] = 1;
}
int main(){
init_unit();
while(cin >> n >> s){
len = s.size();
init();
mat ans = pow_mat(a,n);
cls(b.m);
for(int i=0;i<len;i++)
b.m[i][0] = s[i]-'0';
ans = ans*b;
for(int i=0;i<len;i++) printf("%d",ans.m[i][0]);
printf("\n");
}
return 0;
}
姿势2:
/*************************************************************************
> File Name: hdu2276t2.cpp
> Author: WArobot
> Blog: http://www.cnblogs.com/WArobot/
> Created Time: 2017年05月04日 星期四 14时57分05秒
************************************************************************/
#include<bits/stdc++.h>
using namespace std;
int n,len;
string s;
const int MOD = 2;
const int maxn = 110;
#define ll long long
#define mod(x) ((x)%MOD)
#define cls(x) memset(x,0,sizeof(x));
struct mat{
int m[maxn][maxn];
}unit;
void init_unit(){
for(int i=0;i<maxn;i++) unit.m[i][i] = 1;
}
// 根据循环矩阵可以优化矩阵乘法
// 因为循环矩阵a[i][i] = a[i-1][i-1] , 所以只需要计算出第0行然后递推剩下的其他行就ok了
mat operator *(mat a,mat b){
mat ret;
cls(ret.m);
for(int i=0;i<len;i++)
for(int j=0;j<len;j++)
ret.m[0][i] = mod( ret.m[0][i] + mod(a.m[0][j]&b.m[j][i]) );
for(int i=1;i<len;i++)
for(int j=0;j<len;j++)
ret.m[i][j] = ret.m[i-1][ (j-1+len)%len ];
return ret;
}
mat pow_mat(mat a,int x){
mat ret = unit;
while(x){
if(x&1) ret = ret*a;
a = a*a;
x >>= 1;
}
return ret;
}
mat a,b;
void init(){
cls(a.m);
a.m[0][0] = a.m[0][len-1] = 1;
for(int i=1;i<len;i++) a.m[i][i-1] = a.m[i][i] = 1;
}
int main(){
init_unit();
while(cin >> n >> s){
len = s.size();
init();
// 使用循环矩阵只能计算矩阵a^n,因为只有矩阵a有循环矩阵的特点
mat ans = pow_mat(a,n);
for(int i=0;i<len;i++){
ll tmp = 0;
for(int j=0;j<len;j++){
tmp = mod( tmp + mod(ans.m[i][j]*(s[j]-'0')) );
}
cout<<tmp;
}
cout<<endl;
}
return 0;
}
姿势3:
/*************************************************************************
> File Name: hdu2276t3.cpp
> Author: WArobot
> Blog: http://www.cnblogs.com/WArobot/
> Created Time: 2017年05月04日 星期四 15时41分45秒
************************************************************************/
#include<bits/stdc++.h>
using namespace std;
int n,len;
string s;
const int MOD = 2;
const int maxn = 110;
#define ll long long
#define mod(x) ((x)%MOD)
#define cls(x) memset(x,0,sizeof(x));
struct mat{
int m[maxn][maxn];
}unit;
void init_unit(){
for(int i=0;i<maxn;i++) unit.m[i][i] = 1;
}
// 根据循环矩阵可以优化矩阵乘法
// 因为循环矩阵a[i][i] = a[i-1][i-1] , 所以只需要计算出第0行然后递推剩下的其他行就ok了
mat operator *(mat a,mat b){
mat ret;
cls(ret.m);
for(int i=0;i<len;i++)
for(int j=0;j<len;j++)
ret.m[0][i] ^= (a.m[0][j] & b.m[j][i]);
for(int i=1;i<len;i++)
for(int j=0;j<len;j++)
ret.m[i][j] = ret.m[i-1][ (j-1+len)%len ];
return ret;
}
mat pow_mat(mat a,int x){
mat ret = unit;
while(x){
if(x&1) ret = ret*a;
a = a*a;
x >>= 1;
}
return ret;
}
mat a,b;
void init(){
cls(a.m);
a.m[0][0] = a.m[0][len-1] = 1;
for(int i=1;i<len;i++) a.m[i][i-1] = a.m[i][i] = 1;
}
int main(){
init_unit();
while(cin >> n >> s){
len = s.size();
init();
// 使用循环矩阵只能计算矩阵a^n,因为只有矩阵a有循环矩阵的特点
mat ans = pow_mat(a,n);
for(int i=0;i<len;i++){
ll tmp = 0;
for(int j=0;j<len;j++)
tmp ^= (ans.m[i][j] & (s[j]-'0'));
printf("%lld",tmp);
}
printf("\n");
}
return 0;
}
HDU 2276 Kiki & Little Kiki 2( 矩阵快速幂 + 循环同构矩阵 )的更多相关文章
- HDU1757-A Simple Math Problem,矩阵快速幂,构造矩阵水过
A Simple Math Problem 一个矩阵快速幂水题,关键在于如何构造矩阵.做过一些很裸的矩阵快速幂,比如斐波那契的变形,这个题就类似那种构造.比赛的时候手残把矩阵相乘的一个j写成了i,调试 ...
- HDU6395-Sequence 矩阵快速幂+除法分块 矩阵快速幂模板
目录 Catalog Solution: (有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 Catalog Problem:Portal传送门 原题目描述在最下面. Solution ...
- HDU_4965 Fast Matrix Calculation 2014多校9 矩阵快速幂+机智的矩阵结合律
一开始看这个题目以为是个裸的矩阵快速幂的题目, 后来发现会超时,超就超在 M = C^(N*N). 这个操作,而C本身是个N*N的矩阵,N最大为1000. 但是这里有个巧妙的地方就是 C的来源其实 ...
- LA 3704 (矩阵快速幂 循环矩阵) Cellular Automaton
将这n个格子看做一个向量,每次操作都是一次线性组合,即vn+1 = Avn,所求答案为Akv0 A是一个n*n的矩阵,比如当n=5,d=1的时候: 不难发现,A是个循环矩阵,也就是将某一行所有元素统一 ...
- 经典矩阵快速幂之一-----poj3233(矩阵套矩阵
题意:给你一个矩阵A,求S=A+A^2+A^3+...+A^k. 其实这个当时我看着毫无头绪,看了他们给的矩阵发现好!精!妙! 我们这样看 是不是有点思路! 没错!就是右上角,我们以此类推可以得到A+ ...
- hdu 2604 Queuing(矩阵快速幂乘法)
Problem Description Queues and Priority Queues are data structures which are known to most computer ...
- (hdu 6030) Happy Necklace 找规律+矩阵快速幂
题目链接 :http://acm.hdu.edu.cn/showproblem.php?pid=6030 Problem Description Little Q wants to buy a nec ...
- HDU 6470 【矩阵快速幂】
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6470 写这道题是为了让自己不要忘记矩阵快速幂如何推出矩阵式子的. 注意 代码是TLE的!! #incl ...
- HDU 5318——The Goddess Of The Moon——————【矩阵快速幂】
The Goddess Of The Moon Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/ ...
随机推荐
- 数据结构(5) 第五天 快速排序、归并排序、堆排序、高级数据结构介绍:平衡二叉树、红黑树、B/B+树
01 上次课程回顾 希尔排序 又叫减少增量排序 increasement = increasement / 3 + 1 02 快速排序思想 思想: 分治法 + 挖坑填数 分治法: 大问题分解成各个小问 ...
- 【JavaScript框架封装】实现一个类似于JQuery的内容框架的封装
// 内容框架 (function (xframe) { // 需要参与链式访问的(必须使用prototype的方式来给对象扩充方法) xframe.extend({ /** * .html()用为读 ...
- npm run build 打包项目,图片等资源使用相对路径会出现路径错误的问题
在build下的utils.js中,3使用 ‘vue-style-loader’ 依赖的地方添加 publicPath: '../../' , 如图:
- [网络流24题] 方格取数问题/骑士共存问题 (最大流->最大权闭合图)
洛谷传送门 LOJ传送门 和太空飞行计划问题一样,这依然是一道最大权闭合图问题 “骑士共存问题”是“方格取数问题”的弱化版,本题解不再赘述“骑士共存问题”的做法 分析题目,如果我们能把所有方格的数都给 ...
- Django入门--模板路径配置及渲染
模板就是前端的页面,Django把html源码写到模板文件中,然后通过特定方法渲染后交给客户端. 模板路径设置方法有两种,分别是在项目目录下设置以及在应用目录下设置. 模板查找顺序:优先在DIRS设置 ...
- Oracle学习总结(6)—— SQL注入技术
不管用什么语言编写的Web应用,它们都用一个共同点,具有交互性并且多数是数据库驱动.在网络中,数据库驱动的Web应用随处可见,由此而存在的SQL注入是影响企业运营且最具破坏性的漏洞之一. SQL注入基 ...
- OA项目总结3
struts2自定义标签中 使用in 判断当前值 是否在某个集合中 该属性一方面可以获取前端页面传递过来的参数 另外一个作用就是在数据回显时把用户已经拥有的权限id存入该集合中 放在栈顶 ...
- js日期原型扩展
当初做统计业务需要处理时间 周报:本周 上周 下周 近一周 月报上月 本月 等 需要使用时间处理 所以扩展了这些方法 <!DOCTYPE html> <html xmlns=&quo ...
- POJ 1948
这道题我记得是携程比赛上的一道. 开始时想直接设面积,但发现不可以,改设能否构成三角形.设dp[i][j][k]为前i根木棍构成边长为j和k的三角形,那么转移可以为dp[i][j][k]=dp[i-1 ...
- UI组件之AdapterView及其子类(三)Spinner控件具体解释
Spinner提供了从一个数据集合中高速选择一项值的办法. 默认情况下Spinner显示的是当前选择的值.点击Spinner会弹出一个包括全部可选值的dropdown菜单或者一个dialog对话框,从 ...