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函数调用了多少次. 我们定义数组为所求 ...
随机推荐
- 《Python 二三事》——python学习必看(转载)
面向初学者介绍Python相关的一些工具,以及可能遇到的常见问题. 原文出处 原文作者:八八年出生的男性,互联网上常用id是 jagttt .目前正从事 IT 行业的工作.业余爱好是动漫游加电 ...
- C语言结构体(struct)使用方法
基本定义:结构体,通俗讲就像是打包封装,把一些变量有共同特征(比如同属于某一类事物的属性)的变量封装在内部,通过一定方法访问修改内部变量. 结构体定义: 第一种:只有结构体定义 struct stuf ...
- 我和Cpp的第一次正式约会
今天是我和Cpp的第一次正式的约会,大一的时候学校开了Cpp课,可是那时候玩性未收,没有好好学习,而如今我主动约Cpp,是真的想跟他进一步走下去^_^,正在学习<C++ Primer>,每 ...
- macOS 自动修改mac地址脚本
介于 某公众号提供了通过修改mac地址来链接BUPT_mobile 的推送,上网上查了一下咋写脚本,实现一键修改mac地址的功能 网上有自动修改mac地址的程序,但是很坑爹的要收费,所以不如自力更生写 ...
- Python Lib:pyzmq
http://git.oschina.net/gitlab/StartWithCoding/tree/master/example/python/pyzmq
- 用汇编语言研究C语言的全局变量、局部变量、参数、返回值放在哪里
前提知识 c0s调用main函数的地址: 11ah main函数的连接地址: 01fah 一.全局变量与局部变量 测试程序 int a1,a2,a3; void f(void); void g(voi ...
- 如何编写一个简单的makefile
一个规则的构成 目标:依赖1,依赖2······ 命令 例子: objs := init.o nand.o head.o main.o nand.bin : $(objs) arm-linux-ld ...
- 转:Win7 IIS7应用PHP Manager使用FastCGI通道快速部署PHP支持
原文来自于:http://www.jb51.net/os/windows/62390.html 正常情况下,我们在Windows系统中部署WEB服务器(iis)支持PHP是采用ISAPI通道.参照这篇 ...
- BZOJ 1048 分割矩阵
Description 将一个a*b的数字矩阵进行如下分割:将原矩阵沿某一条直线分割成两个矩阵,再将生成的两个矩阵继续如此分割(当然也可以只分割其中的一个),这样分割了(n-1)次后,原矩阵被分割成了 ...
- MongoDB入门简介
MongoDB入门简介 http://blog.csdn.net/lolinzhang/article/details/4353699 有关于MongoDB的资料现在较少,且大多为英文网站,以上内容大 ...