poj 3233 Matrix Power Series

本文介绍了一种使用二分法优化矩阵快速幂求和的方法,针对给定的矩阵A和正整数k,高效计算S=A+A^2+A^3+…+A^k,并通过代码实现展示了具体的算法流程。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

http://poj.org/problem?id=3233
Description

Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.

Input

The input contains exactly one test case. The first line of input contains three positive integers n (n ≤ 30), k (k ≤ 109) and m (m < 104). Then follow n lines each containing n nonnegative integers below 32,768, giving A’s elements in row-major order.

Output

Output the elements of S modulo m in the same way as A is given.

Sample Input

2 2 4
0 1
1 1
Sample Output

1 2
2 3
题目大意:给你一个矩阵A,算出 S = A + A^2 + A^3 + … + A^k.并输出矩阵S;

解题思路:由于这题的k比较的而时间只有3秒,除了算快速幂还要求累加,这样很容易超时,所以我们要想办法优化。这里我们可以用二分的方法来优化。
如:S(10) = A + A^2 + A^3 + A^4 + A^5+ A^6 + A^7 + A^8 + A^9+ A^10;
= S(5) + A^5 *S(5);
S(5) = A + A^2 + A^3 + A^4 + A^5;
= S(2) + A^3 + A^3*S(2);
S(2) = S(1) + A * S(1);
S(1) = A;
看到这,我们不难发现规律,那就是当k为奇数时s(k) = s(k/2)+A^(k/2+1)+A^(k/2+1)* s(k/2);
当k为偶数时 s(k) = s(k/2)+A^(k/2)* s(k/2);
利用这个规律就可以吧线性的时间复杂度简化为log(n)了。
具体代码实现:

#include <cstdio>
using namespace std;
struct Matrix
{
    int m[35][35];
}a;
int m;
Matrix Mult(Matrix a,Matrix b,int n)
{
    Matrix temp;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)
            temp.m[i][j] = 0;
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)
        for(int k=1;k<=n;k++)
         temp.m[i][j] = (temp.m[i][j]+(a.m[i][k]*b.m[k][j])%m)%m;
    }
    return temp;
}
Matrix Add(Matrix a,Matrix b,int n)
{
    Matrix temp;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)
            temp.m[i][j] = 0;
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)
            temp.m[i][j] = (a.m[i][j]+b.m[i][j])%m;
    }
    return temp;
}
Matrix QuickPower(Matrix a,int n,int k)
{
    Matrix ans,res;
    for(int i=1;i<=n;i++)
       for(int j=1;j<=n;j++)
       if(i==j)ans.m[i][j] = 1;
       else ans.m[i][j] = 0;
    res = a;
    while(k){
        if(k&1)
            ans = Mult(ans,res,n);
        res = Mult(res,res,n);
        k>>=1;
    }
    return ans;
}
Matrix work(int n,int k)
///核心代码,因为这里我们只要二分分到1时可以直接得到答案,所以直到k=1时返回到上一层,
///然后判断这一层的k是计数还是偶数,按照各自的规则算出他们的值然后再返回到上一层,直到结束,这样就可以算出s(k)了。
{
    if(k == 1)
        return a;
    Matrix res = work(n,k/2);
    Matrix temp;
    if(k&1){
        temp = Add(res,QuickPower(a,n,k/2+1),n);
        res = Add(temp,Mult(QuickPower(a,n,k/2+1),res,n),n);
    }
    else{
        res = Add(res,Mult(QuickPower(a,n,k/2),res,n),n);
    }
    return res;
}
int main()
{
    int n,k;
     scanf("%d %d %d",&n,&k,&m);
     for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)
            scanf("%d",&a.m[i][j]);
     }
    Matrix S = work(n,k);
    for(int i=1;i<=n;i++)
    {
        printf("%d",S.m[i][1]);
        for(int j=2;j<=n;j++)
            printf(" %d",S.m[i][j]);
        printf("\n");
    }
    return 0;
}
资源下载链接为: https://pan.quark.cn/s/9648a1f24758 这个HTML文件是一个专门设计的网页,适合在告白或纪念日这样的特殊时刻送给女朋友,给她带来惊喜。它通过HTML技术,将普通文字转化为富有情感和创意的表达方式,让数字媒体也能传递深情。HTML(HyperText Markup Language)是构建网页的基础语言,通过标签描述网页结构和内容,让浏览器正确展示页面。在这个特效网页中,开发者可能使用了HTML5的新特性,比如音频、视频、Canvas画布或WebGL图形,来提升视觉效果和交互体验。 原本这个文件可能是基于ASP.NET技术构建的,其扩展名是“.aspx”。ASP.NET是微软开发的一个服务器端Web应用程序框架,支持多种编程语言(如C#或VB.NET)来编写动态网页。但为了在本地直接运行,不依赖服务器,开发者将其转换为纯静态的HTML格式,只需浏览器即可打开查看。 在使用这个HTML特效页时,建议使用Internet Explorer(IE)浏览器,因为一些老的或特定的网页特效可能只在IE上表现正常,尤其是那些依赖ActiveX控件或IE特有功能的页面。不过,由于IE逐渐被淘汰,现代网页可能不再对其进行优化,因此在其他现代浏览器上运行可能会出现问题。 压缩包内的文件“yangyisen0713-7561403-biaobai(html版本)_1598430618”是经过压缩的HTML文件,可能包含图片、CSS样式表和JavaScript脚本等资源。用户需要先解压,然后在浏览器中打开HTML文件,就能看到预设的告白或纪念日特效。 这个项目展示了HTML作为动态和互动内容载体的强大能力,也提醒我们,尽管技术在进步,但有时复古的方式(如使用IE浏览器)仍能唤起怀旧之情。在准备类似的个性化礼物时,掌握基本的HTML和网页制作技巧非常
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值