poj 3735 （矩阵快速幂）

题意：有n只猫咪，开始时每只猫咪有花生0颗，现有一组操作，由下面三个中的k个操作组成：
               1. g i 给i只猫咪一颗花生米
               2. e i 让第i只猫咪吃掉它拥有的所有花生米
               3. s i j 将猫咪i与猫咪j的拥有的花生米交换

               现将上述一组操作做m次后，问每只猫咪有多少颗花生？

解题报告：http://www.cnblogs.com/acSzz/archive/2012/08/20/2648087.html

// File Name: 3735.cpp
// Author: Missa
// Created Time: 2013/2/7 星期四 12:35:55

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<cstdlib>
#include<map>
using namespace std;
#define ll long long
ll n,m,k;
//**********矩阵快速幂********************
const int maxn=105;//矩阵最大范围
//int n,m;//矩阵n*m
struct Mat
{
    ll mat[maxn][maxn];
    void clear()
    {
        memset(mat,0,sizeof(mat));
    }
    void unit()//单元矩阵n*n矩阵的时候需要
    {
        clear();
        for(int i=0;i<=n;i++)
            mat[i][i]=1;
    }
};
Mat a;
void init()
{
    a.unit();
}
Mat mul(Mat a,Mat b)
{
    Mat c;
    c.clear();
    for(int i=0;i<=n;i++)//n+1 * n+1 的矩阵
    {
        for(int k=0;k<=n;k++)
        {
            if(a.mat[i][k])
            {
                for(int j=0;j<=n;j++)
                    c.mat[i][j]+=a.mat[i][k]*b.mat[k][j];
            }
        }
    }
    return c;
}
Mat mpow(Mat a,ll k)//非递归实现
{
    if(k==1) return a;
    Mat e;
    e.unit();
    while(k)
    {
        if(k&1)     e=mul(e,a);
        k>>=1;
        a=mul(a,a);
    }
    return e;
}
//***********矩阵快速幂*******************

int main()
{
    while(~scanf("%lld%lld%lld",&n,&m,&k))
    {
        if(!n && !m && !k) break;
        init();
        while(k--)
        {
            ll x,y;
            char o[10];
            scanf("%s",&o);
            if(o[0]=='g')
            {
                scanf("%lld",&x);
                a.mat[0][x]++;
            }
            else if(o[0]=='e')
            {
                scanf("%lld",&x);
                for(int i=0;i<=n;i++)
                    a.mat[i][x]=0;
            }
            else
            {
                scanf("%lld%lld",&x,&y);
                for(int i=0;i<=n;i++)
                    swap(a.mat[i][x],a.mat[i][y]);
            }
        }
        if(m==0)
        {
            printf("0");
            for(int i=2;i<=n;i++)
                printf(" 0");
            printf("\n");
            continue;
        }
        a=mpow(a,m);
        printf("%lld",a.mat[0][1]);
        for(int i=2;i<=n;i++)
            printf(" %lld",a.mat[0][i]);
        printf("\n");
    }
    return 0;
}