博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
poj 3735Training little cats(矩阵快速幂 完成模拟 注意优化)
阅读量:4048 次
发布时间:2019-05-25

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

Training little cats
Time Limit: 2000MS   Memory Limit: 65536K
Total Submissions: 12392   Accepted: 3049

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 nm 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

, Facer
/*难点在与怎么将这三种操作转换为矩阵的操作,以题目上的测试数据为例,定义表示猫咪的花生的向量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
#include
#include
#include
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;}

转载地址:http://vafci.baihongyu.com/

你可能感兴趣的文章
uboot start.s文件分析
查看>>
没有路由器的情况下,开发板,虚拟机Ubuntu,win10主机,三者也可以ping通
查看>>
本地服务方式搭建etcd集群
查看>>
安装k8s Master高可用集群
查看>>
忽略图片透明区域的事件(Flex)
查看>>
忽略图片透明区域的事件(Flex)
查看>>
AS3 Flex基础知识100条
查看>>
Flex动态获取flash资源库文件
查看>>
flex4 中创建自定义弹出窗口
查看>>
01Java基础语法-16. while循环结构
查看>>
01Java基础语法-18. 各种循环语句的区别和应用场景
查看>>
01Java基础语法-19. 循环跳转控制语句
查看>>
Django框架全面讲解 -- Form
查看>>
socket,accept函数解析
查看>>
今日互联网关注(写在清明节后):每天都有值得关注的大变化
查看>>
”舍得“大法:把自己的优点当缺点倒出去
查看>>
[今日关注]鼓吹“互联网泡沫,到底为了什么”
查看>>
[互联网学习]如何提高网站的GooglePR值
查看>>
[关注大学生]求职不可不知——怎样的大学生不受欢迎
查看>>
[关注大学生]读“贫困大学生的自白”
查看>>