蒟蒻的我还需深入学习

**链接:****传送门 **

题意:给出一个长度为 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( 矩阵快速幂 + 循环同构矩阵 )的更多相关文章

  1. HDU1757-A Simple Math Problem,矩阵快速幂,构造矩阵水过

    A Simple Math Problem 一个矩阵快速幂水题,关键在于如何构造矩阵.做过一些很裸的矩阵快速幂,比如斐波那契的变形,这个题就类似那种构造.比赛的时候手残把矩阵相乘的一个j写成了i,调试 ...

  2. HDU6395-Sequence 矩阵快速幂+除法分块 矩阵快速幂模板

    目录 Catalog Solution: (有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 Catalog Problem:Portal传送门  原题目描述在最下面. Solution ...

  3. HDU_4965 Fast Matrix Calculation 2014多校9 矩阵快速幂+机智的矩阵结合律

    一开始看这个题目以为是个裸的矩阵快速幂的题目, 后来发现会超时,超就超在  M = C^(N*N). 这个操作,而C本身是个N*N的矩阵,N最大为1000. 但是这里有个巧妙的地方就是 C的来源其实 ...

  4. LA 3704 (矩阵快速幂 循环矩阵) Cellular Automaton

    将这n个格子看做一个向量,每次操作都是一次线性组合,即vn+1 = Avn,所求答案为Akv0 A是一个n*n的矩阵,比如当n=5,d=1的时候: 不难发现,A是个循环矩阵,也就是将某一行所有元素统一 ...

  5. 经典矩阵快速幂之一-----poj3233(矩阵套矩阵

    题意:给你一个矩阵A,求S=A+A^2+A^3+...+A^k. 其实这个当时我看着毫无头绪,看了他们给的矩阵发现好!精!妙! 我们这样看 是不是有点思路! 没错!就是右上角,我们以此类推可以得到A+ ...

  6. hdu 2604 Queuing(矩阵快速幂乘法)

    Problem Description Queues and Priority Queues are data structures which are known to most computer ...

  7. (hdu 6030) Happy Necklace 找规律+矩阵快速幂

    题目链接 :http://acm.hdu.edu.cn/showproblem.php?pid=6030 Problem Description Little Q wants to buy a nec ...

  8. HDU 6470 【矩阵快速幂】

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6470 写这道题是为了让自己不要忘记矩阵快速幂如何推出矩阵式子的. 注意 代码是TLE的!! #incl ...

  9. 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/ ...

随机推荐

  1. 数据结构(5) 第五天 快速排序、归并排序、堆排序、高级数据结构介绍:平衡二叉树、红黑树、B/B+树

    01 上次课程回顾 希尔排序 又叫减少增量排序 increasement = increasement / 3 + 1 02 快速排序思想 思想: 分治法 + 挖坑填数 分治法: 大问题分解成各个小问 ...

  2. 【JavaScript框架封装】实现一个类似于JQuery的内容框架的封装

    // 内容框架 (function (xframe) { // 需要参与链式访问的(必须使用prototype的方式来给对象扩充方法) xframe.extend({ /** * .html()用为读 ...

  3. npm run build 打包项目,图片等资源使用相对路径会出现路径错误的问题

    在build下的utils.js中,3使用 ‘vue-style-loader’ 依赖的地方添加 publicPath: '../../' , 如图:

  4. [网络流24题] 方格取数问题/骑士共存问题 (最大流->最大权闭合图)

    洛谷传送门 LOJ传送门 和太空飞行计划问题一样,这依然是一道最大权闭合图问题 “骑士共存问题”是“方格取数问题”的弱化版,本题解不再赘述“骑士共存问题”的做法 分析题目,如果我们能把所有方格的数都给 ...

  5. Django入门--模板路径配置及渲染

    模板就是前端的页面,Django把html源码写到模板文件中,然后通过特定方法渲染后交给客户端. 模板路径设置方法有两种,分别是在项目目录下设置以及在应用目录下设置. 模板查找顺序:优先在DIRS设置 ...

  6. Oracle学习总结(6)—— SQL注入技术

    不管用什么语言编写的Web应用,它们都用一个共同点,具有交互性并且多数是数据库驱动.在网络中,数据库驱动的Web应用随处可见,由此而存在的SQL注入是影响企业运营且最具破坏性的漏洞之一. SQL注入基 ...

  7. OA项目总结3

    struts2自定义标签中  使用in  判断当前值 是否在某个集合中 该属性一方面可以获取前端页面传递过来的参数   另外一个作用就是在数据回显时把用户已经拥有的权限id存入该集合中   放在栈顶 ...

  8. js日期原型扩展

    当初做统计业务需要处理时间 周报:本周 上周 下周 近一周 月报上月 本月 等 需要使用时间处理 所以扩展了这些方法 <!DOCTYPE html> <html xmlns=&quo ...

  9. POJ 1948

    这道题我记得是携程比赛上的一道. 开始时想直接设面积,但发现不可以,改设能否构成三角形.设dp[i][j][k]为前i根木棍构成边长为j和k的三角形,那么转移可以为dp[i][j][k]=dp[i-1 ...

  10. UI组件之AdapterView及其子类(三)Spinner控件具体解释

    Spinner提供了从一个数据集合中高速选择一项值的办法. 默认情况下Spinner显示的是当前选择的值.点击Spinner会弹出一个包括全部可选值的dropdown菜单或者一个dialog对话框,从 ...