zoukankan      html  css  js  c++  java
  • poj3757 Training little cats

    Time Limit: 2000MS   Memory Limit: 65536K
    Total Submissions: 11496   Accepted: 2815

    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 klines 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 6
    g 1
    g 2
    g 2
    s 1 2
    g 3
    e 2
    0 0 0

    Sample Output

    2 0 1

    题意:

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

                   m次后,每只猫咪有多少颗花生?

    可以构建一个1*(n+1)大小的辅助矩阵,即1 0 0 0,然后根据操作构造转置矩阵。

    转置矩阵的构造:

    转置矩阵一开始初始化为(n+1)*(n+1)大小的单位矩阵,然后每一次操作都要变化。

    1.g i 第0行第i列的元素加1

    2.e i 第i列的元素都变为0

    3.s i j 第i列和第j列的元素都换一下


    #include<iostream>
    #include<stdio.h>
    #include<stdlib.h>
    #include<string.h>
    #include<math.h>
    #include<vector>
    #include<map>
    #include<set>
    #include<queue>
    #include<stack>
    #include<string>
    #include<algorithm>
    using namespace std;
    typedef long long ll;
    #define inf 99999999
    #define pi acos(-1.0)
    
    struct matrix{
        ll n,m,i;
        ll data[105][105];
        void init_danwei(){
            for(i=0;i<n;i++){
                data[i][i]=1;
            }
        }
    };
    
    matrix multi(matrix &a,matrix &b){
        ll i,j,k;
        matrix temp;
        temp.n=a.n;
        temp.m=b.m;
        for(i=0;i<temp.n;i++){
            for(j=0;j<temp.m;j++){
                temp.data[i][j]=0;
            }
        }
        for(i=0;i<a.n;i++){
            for(k=0;k<a.m;k++){
                if(a.data[i][k]>0){
                    for(j=0;j<b.m;j++){
                        temp.data[i][j]=temp.data[i][j]+a.data[i][k]*b.data[k][j];
                    }
                }
            }
        }
        return temp;
    }
    
    matrix fast_mod(matrix &a,ll n){
        matrix ans;
        ans.n=a.n;
        ans.m=a.m;
        memset(ans.data,0,sizeof(ans.data));
        ans.init_danwei();
        while(n>0){
            if(n&1)ans=multi(ans,a);
            a=multi(a,a);
            n>>=1;
        }
        return ans;
    }
    
    int main()
    {
        ll n,k,m,i,j,e,c,d,h;
        while(scanf("%lld%lld%lld",&n,&m,&k)!=EOF)
        {
            if(n==0 && m==0 && k==0)break;
            matrix a;
            a.n=a.m=n+1;
            memset(a.data,0,sizeof(a.data));
            a.init_danwei();
    
    
            char s[10];
            ll temp;
            for(i=1;i<=k;i++){
                scanf("%s",s);
                if(s[0]=='g'){
                    scanf("%lld",&c);
                    a.data[0][c]++;
                }
                else if(s[0]=='s'){
                    scanf("%lld%lld",&c,&d);
                    for(j=0;j<n+1;j++){
                        temp=a.data[j][c];
                        a.data[j][c]=a.data[j][d];
                        a.data[j][d]=temp;
                    }
                }
                else if(s[0]=='e'){
                    scanf("%lld",&c);
                    for(j=0;j<n+1;j++){
                        a.data[j][c]=0;
                    }
    
                }
            }
    
    
            matrix cnt;
            cnt=fast_mod(a,m);
            matrix ant;
            ant.n=1;
            ant.m=n+1;
            memset(ant.data,0,sizeof(ant.data));
            ant.data[0][0]=1;
    
            matrix juzhen;
            juzhen=multi(ant,cnt);
    
    
            for(i=1;i<=n;i++){
                if(i==n)printf("%lld
    ",juzhen.data[0][n]);
                else printf("%lld ",juzhen.data[0][i]);
            }
    
        }
        return 0;
    }
    
    
    
    
    
    
    



  • 相关阅读:
    [导入]郎咸平在复旦大学的精彩演讲
    [导入]郎咸平介绍
    [导入]郎咸平网络文摘 在线交流论坛
    Delphi7,创建 RGB CMYK颜色分量圆 完整代码
    Excel分类汇总/平均/计数的办法
    Word插入 分节符实现多个页码
    Delphi 7中如何修改Form1: TForm1;
    在delphi 7里怎么给已有的单元文件改名
    Mac下安装mysql后无法启动,提示. ERROR! The server quit without updating PID file(就是账户权限不够)
    pip虚拟环境下倒出安装文件,或者批量导入安装。
  • 原文地址:https://www.cnblogs.com/herumw/p/9464602.html
Copyright © 2011-2022 走看看