博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
poj 3735 大数量反复操作问题(矩阵高速幂)
阅读量:4322 次
发布时间:2019-06-06

本文共 1764 字,大约阅读时间需要 5 分钟。

题意:一个一维数组,3种操作: a:  第i个数+1,b: 第i个数=0 ,c::交换某俩处的数。  由三种基本操作构成一组序列,反复该序列m次(m<10^9),问结果

属于一种综合操作反复型: 每次乘以一矩阵T,相当于做一次操作。关键是构造这个矩阵。

先构造最初矩阵A :  1*(n +1) ={1,0,0,0...} ,  第一个一时为了操作第一行数的,   

 T的构造:初始是N+1 * N+1单位阵 这样恰好操作第i个数, +1,就在第0行的第 i个加1;交换就相应列交换,清零就相应列清0.

ans= A*(T^m); 注意用;long long 

#include
#include
using namespace std;struct juz{ long long bat[105][105]; int x,y; //行 列 juz () { memset(bat,0,sizeof(bat)); x=0;y=0; }};juz mutp(juz a,juz b){ juz c; c.x=a.x;c.y=b.y; memset(c.bat,0,sizeof(c.bat)); for(int k=0;k
=1) { if(k%2) c=mutp(c,a); k=k/2; a=mutp(a,a); } return c;}int main(){ int n,m,k; while(cin>>n>>m>>k&&(n||m||k)) { juz a,b,c; a.x=1;a.y=n+1; b.x=n+1;b.y=n+1; for(int i=0;i<=n;i++) { a.bat[0][i]=0; b.bat[i][i]=1; } a.bat[0][0]=1; char tmp; int xx,yy; for(int i=0;i
>tmp; if(tmp=='g') { cin>>xx; b.bat[0][xx]++; } else if(tmp=='e') { cin>>xx; for(int i=0;i<=n;i++) b.bat[i][xx]=0; } else { cin>>xx>>yy; for(int i=0;i<=n;i++) { int tx=b.bat[i][xx]; b.bat[i][xx]=b.bat[i][yy]; b.bat[i][yy]=tx; } } } c=quickf(b,m); c=mutp(a,c); for(int i=1;i<=n;i++) if(i!=n)cout<
<<" "; else cout<
<

转载于:https://www.cnblogs.com/blfshiye/p/4021461.html

你可能感兴趣的文章
“==”运算符与equals()
查看>>
单工、半双工和全双工的定义
查看>>
Hdu【线段树】基础题.cpp
查看>>
时钟系统
查看>>
BiTree
查看>>
5个基于HTML5的加载动画推荐
查看>>
水平权限漏洞的修复方案
查看>>
静态链接与动态链接的区别
查看>>
Android 关于悬浮窗权限的问题
查看>>
如何使用mysql
查看>>
linux下wc命令详解
查看>>
敏捷开发中软件测试团队的职责和产出是什么?
查看>>
在mvc3中使用ffmpeg对上传视频进行截图和转换格式
查看>>
python的字符串内建函数
查看>>
Spring - DI
查看>>
微软自己的官网介绍 SSL 参数相关
查看>>
Composite UI Application Block (CAB) 概念和术语
查看>>
64位MATLAB和C混合编程以及联合调试
查看>>
原生js大总结二
查看>>
PHP基础
查看>>