题意:一个一维数组,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 大数量反复操作问题(矩阵高速幂)的更多相关文章

  1. poj 2778 AC自己主动机 + 矩阵高速幂

    // poj 2778 AC自己主动机 + 矩阵高速幂 // // 题目链接: // // http://poj.org/problem?id=2778 // // 解题思路: // // 建立AC自 ...

  2. POJ 3613 Cow Relays (floyd + 矩阵高速幂)

    题目大意: 求刚好经过K条路的最短路 我们知道假设一个矩阵A[i][j] 表示表示 i-j 是否可达 那么 A*A=B  B[i][j]  就表示   i-j 刚好走过两条路的方法数 那么同理 我们把 ...

  3. poj 3735 Training little cats(矩阵快速幂,模版更权威,这题数据很坑)

    题目 矩阵快速幂,这里的模版就是计算A^n的,A为矩阵. 之前的矩阵快速幂貌似还是个更通用一些. 下面的题目解释来自 我只想做一个努力的人 @@@请注意 ,单位矩阵最初构造 行和列都要是(猫咪数+1) ...

  4. POJ 3735 Training little cats(矩阵快速幂)

    Training little cats Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 11787 Accepted: 2892 ...

  5. HDU2842-Chinese Rings(递推+矩阵高速幂)

    pid=2842">题目链接 题意:求出最少步骤解出九连环. 取出第k个的条件是,k-2个已被取出,k-1个仍在支架上. 思路:想必九连环都玩过吧,事实上最少步骤就是从最后一个环開始. ...

  6. [POJ 3735] Training little cats (结构矩阵、矩阵高速功率)

    Training little cats Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 9613   Accepted: 2 ...

  7. poj 3233(矩阵高速幂)

    题目链接:http://poj.org/problem?id=3233. 题意:给出一个公式求这个式子模m的解: 分析:本题就是给的矩阵,所以非常显然是矩阵高速幂,但有一点.本题k的值非常大.所以要用 ...

  8. [POJ 3150] Cellular Automaton (矩阵高速幂 + 矩阵乘法优化)

    Cellular Automaton Time Limit: 12000MS   Memory Limit: 65536K Total Submissions: 3048   Accepted: 12 ...

  9. hdu 3221 Brute-force Algorithm(高速幂取模,矩阵高速幂求fib)

    http://acm.hdu.edu.cn/showproblem.php?pid=3221 一晚上搞出来这么一道题..Mark. 给出这么一个程序.问funny函数调用了多少次. 我们定义数组为所求 ...

随机推荐

  1. 《Python 二三事》——python学习必看(转载)

        面向初学者介绍Python相关的一些工具,以及可能遇到的常见问题. 原文出处 原文作者:八八年出生的男性,互联网上常用id是 jagttt .目前正从事 IT 行业的工作.业余爱好是动漫游加电 ...

  2. C语言结构体(struct)使用方法

    基本定义:结构体,通俗讲就像是打包封装,把一些变量有共同特征(比如同属于某一类事物的属性)的变量封装在内部,通过一定方法访问修改内部变量. 结构体定义: 第一种:只有结构体定义 struct stuf ...

  3. 我和Cpp的第一次正式约会

    今天是我和Cpp的第一次正式的约会,大一的时候学校开了Cpp课,可是那时候玩性未收,没有好好学习,而如今我主动约Cpp,是真的想跟他进一步走下去^_^,正在学习<C++ Primer>,每 ...

  4. macOS 自动修改mac地址脚本

    介于 某公众号提供了通过修改mac地址来链接BUPT_mobile 的推送,上网上查了一下咋写脚本,实现一键修改mac地址的功能 网上有自动修改mac地址的程序,但是很坑爹的要收费,所以不如自力更生写 ...

  5. Python Lib:pyzmq

    http://git.oschina.net/gitlab/StartWithCoding/tree/master/example/python/pyzmq

  6. 用汇编语言研究C语言的全局变量、局部变量、参数、返回值放在哪里

    前提知识 c0s调用main函数的地址: 11ah main函数的连接地址: 01fah 一.全局变量与局部变量 测试程序 int a1,a2,a3; void f(void); void g(voi ...

  7. 如何编写一个简单的makefile

    一个规则的构成 目标:依赖1,依赖2······ 命令 例子: objs := init.o nand.o head.o main.o nand.bin : $(objs) arm-linux-ld ...

  8. 转:Win7 IIS7应用PHP Manager使用FastCGI通道快速部署PHP支持

    原文来自于:http://www.jb51.net/os/windows/62390.html 正常情况下,我们在Windows系统中部署WEB服务器(iis)支持PHP是采用ISAPI通道.参照这篇 ...

  9. BZOJ 1048 分割矩阵

    Description 将一个a*b的数字矩阵进行如下分割:将原矩阵沿某一条直线分割成两个矩阵,再将生成的两个矩阵继续如此分割(当然也可以只分割其中的一个),这样分割了(n-1)次后,原矩阵被分割成了 ...

  10. MongoDB入门简介

    MongoDB入门简介 http://blog.csdn.net/lolinzhang/article/details/4353699 有关于MongoDB的资料现在较少,且大多为英文网站,以上内容大 ...