poj 3735 大数量反复操作问题(矩阵高速幂)
题意:一个一维数组,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<iostream>
#include<cstring>
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<a.y;k++)
for(int i=0;i<a.x;i++)
if(a.bat[i][k])
{
for(int j=0;j<b.y;j++)
{
c.bat[i][j]+=(a.bat[i][k]*b.bat[k][j]);
}
}
return c;
}
juz quickf(juz a,int k)
{
juz c=a;
for(int i=0;i<a.x;i++)
for(int j=0;j<a.x;j++)
c.bat[i][j]=(i==j);
while(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<k;i++)
{
cin>>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<<c.bat[0][i]<<" ";
else cout<<c.bat[0][i]<<endl;
}
return 0;
}
poj 3735 大数量反复操作问题(矩阵高速幂)的更多相关文章
- poj 2778 AC自己主动机 + 矩阵高速幂
// poj 2778 AC自己主动机 + 矩阵高速幂 // // 题目链接: // // http://poj.org/problem?id=2778 // // 解题思路: // // 建立AC自 ...
- POJ 3613 Cow Relays (floyd + 矩阵高速幂)
题目大意: 求刚好经过K条路的最短路 我们知道假设一个矩阵A[i][j] 表示表示 i-j 是否可达 那么 A*A=B B[i][j] 就表示 i-j 刚好走过两条路的方法数 那么同理 我们把 ...
- poj 3735 Training little cats(矩阵快速幂,模版更权威,这题数据很坑)
题目 矩阵快速幂,这里的模版就是计算A^n的,A为矩阵. 之前的矩阵快速幂貌似还是个更通用一些. 下面的题目解释来自 我只想做一个努力的人 @@@请注意 ,单位矩阵最初构造 行和列都要是(猫咪数+1) ...
- POJ 3735 Training little cats(矩阵快速幂)
Training little cats Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 11787 Accepted: 2892 ...
- HDU2842-Chinese Rings(递推+矩阵高速幂)
pid=2842">题目链接 题意:求出最少步骤解出九连环. 取出第k个的条件是,k-2个已被取出,k-1个仍在支架上. 思路:想必九连环都玩过吧,事实上最少步骤就是从最后一个环開始. ...
- [POJ 3735] Training little cats (结构矩阵、矩阵高速功率)
Training little cats Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 9613 Accepted: 2 ...
- poj 3233(矩阵高速幂)
题目链接:http://poj.org/problem?id=3233. 题意:给出一个公式求这个式子模m的解: 分析:本题就是给的矩阵,所以非常显然是矩阵高速幂,但有一点.本题k的值非常大.所以要用 ...
- [POJ 3150] Cellular Automaton (矩阵高速幂 + 矩阵乘法优化)
Cellular Automaton Time Limit: 12000MS Memory Limit: 65536K Total Submissions: 3048 Accepted: 12 ...
- hdu 3221 Brute-force Algorithm(高速幂取模,矩阵高速幂求fib)
http://acm.hdu.edu.cn/showproblem.php?pid=3221 一晚上搞出来这么一道题..Mark. 给出这么一个程序.问funny函数调用了多少次. 我们定义数组为所求 ...
随机推荐
- Winform控件Enable=false显示优化
在B/S开发中(ASP.NET),往往可以css样式表来让页面控件更加美观,但是在C/S中(Winform)里面,我们则需要通过其他取巧的 方式来实现.例如:当你因为某个需求需要将控件设置为Reado ...
- ThinkPHP接入支付宝支付功能
最近做系统,需要实现在线支付功能,毫不犹豫,选择的是支付宝的接口支付功能.这里我用的是即时到帐的接口,具体实现的步骤如下: 一.下载支付宝接口包 下载地址:https://b.alipay.com/o ...
- 在Linux终端执行clear或top命令时出现:'xterm': unknown terminal type
在Linux终端执行clear或top命令时出现:'xterm': unknown terminal type的错误. 例如: [root@localhost phpmyadmin]# clear ' ...
- 如何在Android Studio中使用Gradle发布项目至Jcenter仓库
简述 目前非常流行将开源库上传至Jcenter仓库中,使用起来非常方便且易于维护,特别是在Android Studio环境中,只需几步配置就可以轻松实现上传和发布. Library的转换和引用 博主的 ...
- UIview lianxi
// 创建一个和屏幕大小相同的window,记住[UIScreen mainScreen].bounds 是获取当前屏幕大小 self.window = [[[UIWindow alloc] init ...
- MIT教授将网页开发整合为完整独立的程式语言Ur/Web
MIT 的软体技术教授 Adam Chlipala 设计了新的 Ur/Web 程式语言,这是一个整合 HTML.CSS.XML.SQL 及 JavaScript 等网路标准的“完整独立”语言,强调快速 ...
- 通用GPIO模拟串口,提供源代码,本人经过测试OK(第一版)
--------------------------serial.h------------------------------------------ #ifndef _SERIAL_H_ #def ...
- 解决VS2010中产生的ipch文件夹和sdf文件
纯粹参考他人博文,转载时注明参考链接!!! 参考链接: http://blog.sina.com.cn/s/blog_707a4cbf0100r5jq.html 在VS2010开发的项目中,ipch文 ...
- Entropy (huffman) 优先队列)
http://acm.hdu.edu.cn/showproblem.php?pid=1053 Huffman问题利用STL中的priority_queue解决: #include<stdio.h ...
- BZOJ 1029 [JSOI2007]建筑抢修 已更新
1029: [JSOI2007]建筑抢修 Time Limit: 4 Sec Memory Limit: 162 MBSubmit: 2748 Solved: 1213[Submit][Statu ...