poj 3735Training little cats(矩阵快速幂 完成模拟 注意优化)
来源:互联网 发布:sas聚类分析共现矩阵 编辑:程序博客网 时间:2024/05/21 22:48
Description
Facer's pet cat just gave birth to a brood of little cats. Having considered the health of those lovely cats, Facer decides to make the cats to do some exercises. Facer has well designed a set of moves for his cats. He is now asking you to supervise the cats to do his exercises. Facer's great exercise for cats contains three different moves:
g i : Let the ith cat take a peanut.
e i : Let the ith cat eat all peanuts it have.
s i j : Let the ith cat and jth cat exchange their peanuts.
All the cats perform a sequence of these moves and must repeat it m times! Poor cats! Only Facer can come up with such embarrassing idea.
You have to determine the final number of peanuts each cat have, and directly give them the exact quantity in order to save them.
Input
The input file consists of multiple test cases, ending with three zeroes "0 0 0". For each test case, three integers n, m and k are given firstly, where n is the number of cats and k is the length of the move sequence. The following k lines describe the sequence.
(m≤1,000,000,000, n≤100, k≤100)
Output
For each test case, output n numbers in a single line, representing the numbers of peanuts the cats have.
Sample Input
3 1 6g 1g 2g 2s 1 2g 3e 20 0 0
Sample Output
2 0 1
Source
/*难点在与怎么将这三种操作转换为矩阵的操作,以题目上的测试数据为例,定义表示猫咪的花生的向量A:1 0 0 0A[id]表示第id只猫的初始花生数量,初始时猫咪没有花生,所以都是0。当然,首个元素固定是1,这是为了做乘法用的。定义单位初始方阵T:1 0 0 00 1 0 00 0 1 00 0 0 1定义 R = A * T,那么R[id]表示第id只猫的最终花生数量。上面向量A那个突兀的1就是为了取矩阵T第1行各个值的。那么g i := ++T[0][i]e i := set T[i][i] = 0 即将矩阵的元素(i,i)置为0s i j := swap(T[t][i], T[t][j]) (0<=t<=height) 即将矩阵的i列和j列互换。T执行m次,结果R = A * Tm 。轮到矩阵的幂出场了,将上次的模板改改即可:*/ #include<iostream>#include<string>#include<cstring>#include<cstdio>typedef long long ll;using namespace std;struct matrix{ll mat[150][150];matrix(){memset(mat,0,sizeof(mat));}};matrix ret,d;ll n,m,k;matrix mul(matrix a,matrix b){memset(ret.mat,0,sizeof(ret.mat));for(int i=0;i<=n;i++){for(int k=0;k<=n;k++)//优化 防止超时 if(a.mat[i][k]){for(int j=0;j<=n;j++){ret.mat[i][j]+=a.mat[i][k]*b.mat[k][j];}}}return ret;}matrix ks(matrix a,long long k){if(k==1)return a;matrix e;for(int i=0;i<=n;i++){e.mat[i][i]=1;}if(k==0)return e;while(k){if(k&1){e=mul(e,a);}a=mul(a,a);k>>=1;}return e;}int main(){while(~scanf("%lld%lld%lld",&n,&m,&k)){if(!n&&!m&&!k)break;ll a,b;memset(d.mat,0,sizeof(d.mat));for(int i=0;i<=n;i++)d.mat[i][i]=1;char op[5];while(k--){scanf("%s",op);if(op[0]=='g')// 第一行表示每只猫的豆豆数 当输入g的时候a猫的豆豆加一 当输入e时 把a猫的所有豆豆清空 交换时把两列的数据全部交换 {scanf("%lld",&a);d.mat[0][a]++;}else if(op[0]=='e'){scanf("%lld",&a);for(int i=0;i<=n;i++)d.mat[i][a]=0;}else {scanf("%lld%lld",&a,&b);for(int i=0;i<=n;i++){ll t=d.mat[i][a];d.mat[i][a]=d.mat[i][b];d.mat[i][b]=t;}}}matrix ans=ks(d,m);printf("%lld",ans.mat[0][1]);for(int i=2;i<=n;i++){printf(" %lld",ans.mat[0][i]);}printf("\n");}return 0;}
- poj 3735Training little cats(矩阵快速幂 完成模拟 注意优化)
- poj 3735 Training little cats(矩阵快速幂)
- poj 3735 Training little cats(矩阵构造,快速幂)
- POJ 3735 Training little cats(矩阵快速幂)
- POJ 3735 Training little cats(矩阵快速幂)
- POJ 3735 Training little cats(矩阵快速幂)
- POJ - 3735 Training little cats (矩阵快速幂)
- poj 3735Training little cats (矩阵快速幂)
- POJ 3735 Training little cats 解题报告(矩阵构造+快速幂优化)
- POJ 3745 : Training little cats(矩阵快速幂加优化)
- POJ-3735-Training little cats-构造矩阵+矩阵快速幂+稀疏矩阵乘法优化
- POJ 3735 Training little cats 矩阵二分幂 + 矩阵优化
- POJ--3735[Training little cats] 矩阵建模+快速幂
- poj 3735 Training little cats(构造矩阵快速幂)
- POJ 3735 Training little cats 矩阵快速幂应用
- POJ 3735 Training little cats (矩阵快速幂)
- poj-3735-Training little cats-矩阵快速幂
- POJ 3735 Training little cats (矩阵快速幂)
- 关于处理Android 异常NetworkOnMainThreadException解决方案
- Android工具类--日志管理类
- 人人皆懂得KMP算法
- UVALive - 4256 || LA 4256 Salesmen 商人(DP)
- SVM边学边总结系列——非线性可分情况
- poj 3735Training little cats(矩阵快速幂 完成模拟 注意优化)
- java基础知识点
- 为什么32位机器最大只能用到4GB内存
- SpringMVC之后台接收参数与前台传递数据
- hdfs 机架感知
- C#中的泛型
- LeetCode Learning 1
- 数据源,连接池,数据库连接之间的关系
- cocos2d-android - 1.地图动态加载原理(世界地图的实现)