设为首页 加入收藏

TOP

HDU 4965 Fast Matrix Calculation 矩阵快速幂
2015-07-20 17:51:53 来源: 作者: 【 】 浏览:3
Tags:HDU 4965 Fast Matrix Calculation 矩阵 快速

?

题意:一个矩阵N*K的矩阵A,一个K*N的矩阵B,(4 <= N <= 1000,2 <=K <= 6)。先进行矩阵乘法C=A*B,然后进行D=C^(N*N)的幂运算,然后对D的每个数对6取模,最后求出D的所有位置数字之和。

思路:像之前那道矩阵乘法一样,特别大的矩阵直接进行乘法在没有小规律的帮助时是不可能直接过的(目前看即使是Strassen矩阵算法也不会加速到要求以内)题目中给的C矩阵是1000*1000的矩阵进行快速幂是一定超时的,所以我注意到了A矩阵的列数和B矩阵的行数K是不超过六的,而6*6的矩阵进行快速幂是没问题的。D=A*B*A*B*...*A*B*A*B,可以看成D=A*(B*A*...*B*A)*B。这样既不影响矩阵左乘和右乘的关系,又能防止行列数太大爆掉。(题里的小细节是挺厉害的,首先爆栈的问题需要解决,然后还有矩阵的赋值运算也是O(N)的复杂度)

代码:

?

#pragma comment(linker, /STACK:1024000000,1024000000)
#include 
  
   
#include 
   
     #include 
    
      #include 
     
       #include 
      
        #include 
       
         #include 
        
          #include 
         
           #include
           #include 
           
             #include 
            
              #include 
             
               #include 
              
                #include 
               
                 #define eps 1e-8 #define INF 0x7fffffff #define PI acos(-1.0) #define seed 31//131,1313 typedef long long LL; typedef unsigned long long ULL; #define MOD 6 #define maxn 1005 #define maxm 1005 using namespace std; struct Matrix { int n,m; int a[maxn][maxm]; void init() { n=m=0; memset(a,0,sizeof(a)); } Matrix operator +(const Matrix &b) const { Matrix tmp; tmp.n=n; tmp.m=m; for(int i=0; i
                
                 >=1; m.Copy(m*m); } return tmp; } Matrix a,b; int main() { int n,k; while(scanf(%d%d,&n,&k)) { if(n==0&&k==0) break; a.n=n; a.m=k; for(int i=0; i
                 
                  

?

】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
分享到: 
上一篇hdu 4966 最小树形图 下一篇HDU4968-Improving the GPA

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容:

·Announcing October (2025-12-24 15:18:16)
·MySQL有什么推荐的学 (2025-12-24 15:18:13)
·到底应该用MySQL还是 (2025-12-24 15:18:11)
·进入Linux世界大门的 (2025-12-24 14:51:47)
·Download Linux | Li (2025-12-24 14:51:44)