zoj 2671 Cryptography(线段树+矩阵乘法)

原创 2012年03月25日 21:04:22

【题目大意】:给出数字范围r,n个矩阵,和m个查询,每个查询第i个矩阵到第j个矩阵的乘积,且乘后不超过r。


【解题思路】:利用线段树进行优化。

  原来case与case直接的矩阵也要有空行...-_-!!!!!我无限PE啊....


【代码】:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <cmath>
#include <string>
#include <cctype>
#include <map>
#include <iomanip>
                   
using namespace std;
                   
#define eps 1e-8
#define pi acos(-1.0)
#define inf 1<<30
#define linf 1LL<<60
#define pb push_back
#define lc(x) (x << 1)
#define rc(x) (x << 1 | 1)
#define lowbit(x) (x & (-x))
#define ll long long

struct Mart{
    int mat[4][4];
}Mar[2*32768+1];

int r,n,m;

Mart mart_multiply(Mart a,Mart b){  
    Mart c;  
    for(int i=0; i<2; i++)  
        for(int j=0; j<2; j++){  
            c.mat[i][j]=0;  
            for(int k=0; k<2; k++)  
                c.mat[i][j]=(c.mat[i][j]+a.mat[i][k]*b.mat[k][j]%r)%r;  
        }  
    return c;  
} 

Mart query(int p,int le,int ri,int l,int r){
    if (le==l && ri==r) return Mar[p];
    int mid=(le+ri)>>1;
    if (r<=mid) return query(lc(p),le,mid,l,r);
    else if (l>mid) return query(rc(p),mid+1,ri,l,r);
    Mart ans;
    ans=mart_multiply(query(lc(p),le,mid,l,mid),query(rc(p),mid+1,ri,mid+1,r));
    return ans;
}

int main() {
    bool flag=false;
    while (~scanf("%d%d%d",&r,&n,&m)){
        int tmp=1;
        while (tmp<n) tmp<<=1;
        for (int i=0; i<n; i++){
            for (int j=0; j<2; j++){
                for (int l=0; l<2; l++)
                    scanf("%d",&Mar[tmp+i].mat[j][l]);
            }
        }
        for (int i=tmp-1; i>0; i--){
            Mar[i]=mart_multiply(Mar[lc(i)],Mar[rc(i)]);
        }
        for (int i=0; i<m; i++){
            int l,r;
            if (flag) printf("\n");
            else flag=true;
            scanf("%d%d",&l,&r);
            Mart ans;
            ans=query(1,1,tmp,l,r);
            for (int j=0; j<2; j++){
                printf("%d %d\n",ans.mat[j][0],ans.mat[j][1]);
            }
        }
    }    
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

zoj 2671 Cryptography(矩阵相乘+线段树)

Cryptography Time Limit: 5 Seconds      Memory Limit: 32768 KB Young cryptoanalyst Georgie ...

【ZOJ】2671 Cryptography 线段树

传送门:【ZOJ】2671 Cryptography 题目分析:线段树水题。

ZOJ 3772 Calculate the Function [线段树+矩阵乘法]【思维?】

题目链接:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3772 ————————————————————————————...

zoj3772--【水题】线段树区间查询+矩阵乘法

思路来源:http://blog.csdn.net/u013654696/article/details/23037407#comments 【思路不是自己的,但代码是自己敲的,由于伦家不懂】

zoj 3772 Calculate the Function(线段树+矩阵乘法)

Calculate the Function Time Limit: 2 Seconds      Memory Limit: 65536 KB You are given a ...

ZOJ 2671 线段树 + 矩阵区间乘积

传送门 : zoj 2671summary 下次, 线段树不用cin, cout了有毒啊 AC code:#include #include #include using namespace st...
  • ADjky
  • ADjky
  • 2016-12-10 23:11
  • 221

HDU上的专题训练(背包问题-线段树+树状数组+DP优化+网络流+字符匹配+最短路+矩阵乘法+搜索专题++)

就是DIY上的专题训练 【背包问题】 2602 Bone Collector      1114 Piggy-Bank      1203 I NEED A OFFER! 1171 Big ...

【弱校胡策】2016.4.14 (bzoj2164)最短路+状压DP+矩阵乘法+高斯消元+树链剖分+线段树+背包DP

cyyz&qhyz&lwyz&gryz弱校胡策 命题人:cyyz ws_fqkT3暴力写挫了 50+10+0滚粗辣!奇妙的约会(appointment.cpp/c/pas)【问题描述】DQS和sxb在...

[bzoj 4540] [Hnoi2016]序列:离线,线段树,矩阵乘法

题意:给一个长为N的序列,每项绝对值不超过10^9。Q个询问,求某区间的所有子区间的最小值之和。(1 ≤N,Q ≤ 100000) UPDATE 2017.5.27: 这类离线线段树问题更好的理解方式...

hdu 5068 Harry And Math Teacher ( 线段树+矩阵乘法 )

Harry And Math Teacher Time Limit: 5000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/O...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)