『专题』十个利用矩阵乘法解决的经典题目(转)

最新推荐文章于 2021-12-17 23:29:33 发布
最新推荐文章于 2021-12-17 23:29:33 发布
阅读量416 收藏 1
点赞数
分类专栏： 『数据结构与算法』

转载声明:http://www.cnblogs.com/kuangbin/archive/2012/08/17/2643347.html
（没找到原址,转了kuangbin巨巨的保留而用）

个人感想:最近在敲10道经典矩阵题目,我搞了前9道,最后一道看似看懂了,但也没搞出来,贴下一战绩吧,感觉很有用处.

poj3233

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e5+10;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }


ll N,m,k;

struct mat
{
     ll g[35][35];
     void Clear()
     {
         memset(g,0,sizeof(g));
     }

};

mat mul(mat A,mat B)
{
    mat C;
    C.Clear();
    for(int i=0;i<N;i++)
    {
        for(int k=0;k<N;k++)
        {
            for(int j=0;j<N;j++)
            {
                C.g[i][j]=(C.g[i][j]+A.g[i][k]*B.g[k][j])%m;
            }
        }
    }
    return C;
}
mat cal(mat A,mat B)
{
    mat C;
    C.Clear();

    for(int i=0;i<N;i++)
    {
        for(int j=0;j<N;j++)
        {
            C.g[i][j]=(A.g[i][j]+B.g[i][j])%m;
        }
    }

    return C;
}
mat pow(mat A,ll n)
{
    mat B;
    B.Clear();
    for(int i=0;i<N;i++)
    {
        B.g[i][i]=1;
    }

    while(n>0)
    {
        if(n&1) B=mul(B,A);
        A=mul(A,A);
        n>>=1;
    }
    return B;
}
mat deal(mat x,ll k)
{

    if(k==1) return x;

    mat B=pow(x,(k+1)/2);


    mat T=deal(x,k/2);



    if(k%2==0) return cal(T,mul(B,T));//exp: s(6) = (1 + A^3) * s(3)
    return cal(x,mul(cal(x,B),T));//exp: s(7) = A + (A + A^4) * s(3)
}
void solve()
{
    mat A;
    for(int i=0;i<N;i++)
    {
        for(int j=0;j<N;j++)
        {
            scanf("%I64d",&A.g[i][j]);
        }
    }




    mat P=deal(A,k);

    for(int i=0;i<N;i++)
    {
        for(int j=0;j<N;j++)
        {
            if(j==0)printf("%I64d",P.g[i][j]);
            else printf(" %I64d",P.g[i][j]);
        }
        printf("\n");
    }
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("coco.txt","r",stdin);
    freopen("lala.txt","w",stdout);
#endif
    while(scanf("%I64d%I64d%I64d",&N,&k,&m)!=EOF)
    {
        solve();
    }

    return 0;
}

hdu 2254 奥运

#include <stdio.h>
#include <string.h>
#include <cmath>

#define maxn 30
const int MOD = 2008;
int len;
__int64 city[maxn];
struct Mat
{
    int mat[maxn][maxn];
}e,init;
int FindCity(__int64 x,bool flag)
{
    int i;
    for (i=0;i<len;i++)
        if(city[i]==x)
            return i;
    if (flag)
    {
        city[len++] = x;
        return len-1;
    }
    return 50;/*没找到城市*/
}
Mat operator+(Mat a,Mat b)
{
    Mat c;
    int i,j;
    for (i=0;i<len;i++)
    {
        for(j=0;j<len;j++)
            c.mat[i][j] = (a.mat[i][j]+b.mat[i][j])%MOD;
    }
    return c;
}
Mat operator*(Mat a,Mat b)
{
    int i,j,k;
    Mat c;
    for (i=0;i<len;i++)
    {
        for (j=0;j<len;j++)
        {
            c.mat[i][j] = 0;
            for(k=0;k<len;k++)
                c.mat[i][j]+=(a.mat[i][k]*b.mat[k][j])%MOD;
        }
    }
    return c;
}
Mat operator^(Mat a,int x)
{
    Mat p = e,q = a;
    while (x)
    {
        if(x&1)
            p = p*q;
        x>>=1;
        q = q*q;
    }
    return p;
}

Mat solve(Mat x,int k)
{

    if(k==1) return x;

    Mat B=x^((k+1)/2);


    Mat T=solve(x,k/2);



    //if(k%2==0) return cal(T,mul(B,T));//exp: s(6) = (1 + A^3) * s(3)
    if(k%2==0) return T+(B*T);
    //return cal(x,mul(cal(x,B),T));//exp: s(7) = A + (A + A^4) * s(3)
    return x+((x+B)*T);
}
int main()
{
    __int64 v1,v2;
    int t1,t2,i,j,n,k,c,locate1,locate2,sum;
    Mat tem1,tem2;
    for(i=0;i<maxn;i++)
        for(j=0;j<maxn;j++)
            e.mat[i][j] = (i==j);
    while (scanf("%d",&n)!=EOF)
    {
        len = 0;
        memset(&init,0,sizeof(init));
        while (n--)
        {
            scanf("%I64d%I64d",&v1,&v2);
            locate1 = FindCity(v1,1);
            locate2 = FindCity(v2,1);
            init.mat[locate1][locate2]++;
        }
        scanf("%d",&k);
        while (k--)
        {
            scanf("%I64d%I64d%d%d",&v1,&v2,&t1,&t2);
            if(t1>t2)
            {
                int temp=t1;
                t1=t2;
                t2=temp;
            }

            locate1 = FindCity(v1,0);
            locate2 = FindCity(v2,0);
            if(locate1==50||locate2==50||(t1==0&&t2==0))
                printf("0\n");
            else
            {
                if (t1>1)
                {
                    tem1 = solve(init,t2);
                    tem2 = solve(init,t1-1);
                    sum = (tem1.mat[locate1][locate2]%MOD-tem2.mat[locate1][locate2]%MOD)%MOD;
                }
                else
                {
                    tem1 = solve(init,t2);
                    sum = tem1.mat[locate1][locate2]%MOD;
                }
                if(sum<0)
                    sum+=MOD;
                printf("%d\n",sum);
            }
        }
    }
    return 0;
}

vijosP1067 守望者的烦恼:

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e5+10;
const long long mod=7777777;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }

 typedef vector<ll>vec;
 typedef vector<vec>mat;
 ll K,N;
 mat mul(mat &A,mat & B)
 {
     mat C(A.size(),vec(B[0].size()));

     for(int i=0;i<A.size();i++)
     {
         for(int k=0;k<B.size();k++)
         {
             for(int j=0;j<B[0].size();j++)
             {
                 C[i][j]=(C[i][j]+A[i][k]*B[k][j])%mod;
             }
         }
     }
     return C;
 }
 mat pow(mat A,ll n)
 {
     mat B(A.size(),vec(A.size()));

     for(int i=0;i<A.size();i++)
     {
         B[i][i]=1;
     }
     while(n>0)
     {
         if(n&1) B=mul(B,A);
         A=mul(A,A);
         n>>=1;
     }
     return B;
 }
 void solve()
 {
     mat A(K,vec(K));
     for(int i=0;i<K;i++)
     {
         A[0][i]=1;
     }
     for(int i=1;i<K;i++)
     {
         A[i][i-1]=1;
     }

     A=pow(A,N);

     mat B(K,vec(K));
     B[0][0]=1;

     A=mul(A,B);

     cout<<A[0][0]<<endl;

 }
int main()
{

    while(cin>>K>>N)
    {
        solve();
    }
    return 0;
}

vijosP1049 送给圣诞夜的礼物（m次置换）

/* Author:GavinjouElephant
 * Title:
 * Number:
 * main meanning：
 *
 *
 *
 */


#include <iostream>
using namespace std;
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <cctype>
#include <vector>
#include <set>
#include <cstdlib>
#include <map>
#include <queue>
//#include<initializer_list>
//#include <windows.h>
//#include <fstream>
//#include <conio.h>
#define MaxN 0x7fffffff
#define MinN -0x7fffffff
#define lson 2*k
#define rson 2*k+1
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e5+10;
int Scan()//读入整数外挂.
{
    int res = 0, ch, flag = 0;

    if((ch = getchar()) == '-')             //判断正负
        flag = 1;

    else if(ch >= '0' && ch <= '9')           //得到完整的数
        res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9' )
        res = res * 10 + ch - '0';
    return flag ? -res : res;
}
 void Out(int a)    //输出外挂
 {
     if(a>9)
         Out(a/10);
     putchar(a%10+'0');
 }

 int N,m,k;

 int G[15][105];
 class mat
 {
public:
     int g[105][105];
     void Clear()
     {
         memset(g,0,sizeof(g));
     }
 };
 mat mul(mat A,mat B)
 {
     mat C;
     C.Clear();


     for(int i=0;i<N;i++)
     {
         for(int k=0;k<N;k++)
         {
             for(int j=0;j<N;j++)
             {
                 C.g[i][j]=(C.g[i][j]+A.g[i][k]*B.g[k][j]);
             }
         }
     }
     return C;
 }
 mat pow(mat A,ll n)
 {
     mat B;
     B.Clear();

     for(int i=0;i<N;i++)
     {
         B.g[i][i]=1;
     }
     while(n>0)
     {
         if(n&1) B=mul(B,A);
         A=mul(A,A);
         n>>=1;
     }
     return B;
 }

 void solve()
 {
     int t=k/m;

     mat A;
     A.Clear();

     for(int j=0;j<N;j++)
     {
         A.g[j][G[m][j+1]-1]=1;

        // cout<<G[m][j+1]<<endl;
     }

     A=pow(A,t);






     int l=k%m;//最后的转换

     if(l)
     {
         mat B;
         B.Clear();

         for(int j=0;j<N;j++)
         {
            B.g[j][ G[l][j+1]-1 ]=1;
         }

         A=mul(B,A);
     }


     //ans----->
     mat C;
     C.Clear();

     for(int j=0;j<N;j++)
     {
         C.g[j][0]=j+1;
     }




     A=mul(A,C);

//     cout<<"------->"<<endl;
//     for(int i=0;i<N;i++)
//     {
//         for(int j=0;j<N;j++)
//         {
//             printf("%d ",A.g[i][j]);
//         }
//         printf("\n");
//     }

     for(int j=0;j<N;j++)
     {
         if(j==0)printf("%d",A.g[j][0]);
         else printf(" %d",A.g[j][0]);
     }
     printf("\n");
 }

int main()
{


    while(scanf("%d%d%d",&N,&m,&k)!=EOF)
    {
          memset(G,0,sizeof(G));
          for(int j=1;j<=N;j++)
          {
             G[0][j]=j;
          }

          int x;
          for(int i=1;i<=m;i++)
          {
              for(int j=1;j<=N;j++)
              {
                  scanf("%d",&x);
                  G[i][j]=G[i-1][x];
              }
          }
          solve();
    }

    return 0;
}

ps:多米诺骨牌的判断疑问
1.((~i)&j) == ((~i)&(SIZE-1))
对于第n-1列,那么我们就 例如 （从左往右看 转成从上往下看）{001} -反码代表第n层状态{110} 因为只能打横铺不能打竖,否则就冲突,就判断j这个状态存不存在(110）这个状态.
2. for (int k=0; k<8; k++)isValid=isValid||(i&j)==validSet[k];
我们可以发现kuangbin巨巨写的const int validSet[]={0,3,6,12,15,24,27,30};两个两个1在移动为什么需要那么判呢,主要是判断是否可以放一个打竖的多米诺骨牌,这样的状态是允许的,因为是为了防止状态重复计算,而把打竖的抽离出来再这样判. 然后建图即可.

Gavinjou大笨象
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
十个利用矩阵乘法解决的经典题目
matrix67的专栏
08-04 2478
    好像目前还没有这方面题目的总结。这几天连续看到四个问这类题目的人，今天在这里简单写一下。这里我们不介绍其它有关矩阵的知识，只介绍矩阵乘法和相关性质。    不要以为数学中的矩阵也是黑色屏幕上不断变化的绿色字符。在数学中，一个矩阵说穿了就是一个二维数组。一个n行m列的矩阵可以乘以一个m行p列的矩阵，得到的结果是一个n行p列的矩阵，其中的第i行第j列位置上的数等于前一个矩阵第i行
Matrix67_ 十个利用矩阵乘法解决的经典题目1
08-04
本篇文章主要介绍了十个利用矩阵乘法解决的经典编程题目，通过这些题目，我们可以深入理解矩阵乘法的特性和应用。 首先，矩阵乘法的基本规则是：一个n行m列的矩阵可以与一个m行p列的矩阵相乘，得到一个n行p列的新...
矩阵特征值的求解例子
算法channel
12-07 2万+
请点击上面公众号，免费订阅。　 《实例》阐述算法，通俗易懂，助您对算法的理解达到一个新高度。包含但不限于：经典算法，机器学习，深度学习，LeetCode 题解，Kaggle 实战。期待您的到来！ 01 — 求矩阵特征值的例子 矩阵的特征值为：2，0.4，分别对应的特征向量如上所述。 请记住：每天一小步，日积月累一大步！
KMP算法个人理解
weixin_42250399的博客
12-17 265
参考的原文：详解KMP算法 KMP：解决字符串匹配问题，字符串s是否是主串的子集？是，返回第一个元素坐标；否，返回-1。 核心思想，尽量不移动主串的i指针，回溯子串的j指针。 整个KMP的重点就在于当某一个字符与主串不匹配时，我们应该知道j指针要移动到哪？ 接下来我们自己来发现j的移动规律： 如图：C和D不匹配了，我们要把j移动到哪？显然是第1位。为什么？因为前面有一个A相同啊： 如下图也是一样的情况： 可以把j指针移动到第2位，因为前面有两个字母是一样的： 至此我们可以大概看出一点端倪，当匹配失败
【矩阵乘法专题】解题报告
weixin_30385925的博客
08-09 110
HDU 1575Tr A 　　水题。直接矩阵快速幂后求对角线元素和。 HDU 1757A Simple Math Problem 　　水题。公式非常容易推。 HDU 2604Queuing 　　以mm结尾的只能由fm结尾的或者mm结尾的推来。以mf结尾的只能由mm结尾的推来，以fm结尾的只能由mf或者ff推来，以ff结尾的只能由mf推来。于是定义状态dp[i][j...
超硬核十万字！全网最全 数据结构 代码，随便秒杀老师/面试官，我说的
热门推荐
少说，多做。
04-20 11万+
本文代码实现基本按照《数据结构》课本目录顺序，外加大量的复杂算法实现，一篇文章足够。能换你一个收藏了吧？
力扣经典题目解法记录（更新中）
skkw111的博客
03-25 2325
力扣算法题库引用自：https://cloud.tencent.com/developer/article/1671018 基础 排序算法 快速排序 归并排序 冒泡排序 二分查找 数组 4. 寻找两个正序数组的中位数 33. 搜索旋转排序数组 直接使用二分法 287. 寻找重复数 dict 34. 在排序数组中查找元素的区间 取开始下标(mid = (l + r) // 2); 取结束下标(mid = (a + b + 1) // 2) 矩阵 240. 搜索二维...
矩阵快速幂专题
weixin_43184669的博客
08-26 116
poj3233 经典题 因为k的值太大，如果直接算，肯定会超时。 因此采用二分的方法。 #include <cstdio> #include <cstring> using namespace std; int n,k,m; struct matrix { int mat[35][35]; }; matrix A; void init(matrix &a) ...
专题：强大的运算工具——矩阵
倘若天下着雨。。。你还在
07-10 998
什么是矩阵 对于 矩阵的的初步认识这里有一些讲解。 一开始我对矩阵奇葩的运算方式感到奇怪，它有什么用。。。 DPDP问题的加速可以利用矩阵的结合律来实现。 同时一些图论问题，利用矩阵有关图的联通和边权问题我们又有了一个新招了。（讲道理，要用到矩阵的题目似乎都是为它专门设计的。） 同时还有其他的作用可以参考大牛的博客十个利用矩阵乘法解决的经典题目，codecode 矩阵乘法在图联通中的
十大利用矩阵乘法解决的经典题目
09-08
在计算机科学和编程领域，矩阵乘法是一种基本的数学运算，尤其在图形处理、线性代数、机器...通过阅读提供的《十个利用矩阵乘法解决的经典题目.pdf》文档，你将能更深入地探索这些主题，并在实践中提升你的编程技巧。
ASP基于BS结构的学生在线选课系统的实现(源代码+论文).rar
08-10
资源包主要包含以下内容： ASP项目源码：每个资源包中都包含完整的ASP项目源码，这些源码采用了经典的ASP技术开发，结构清晰、注释详细，帮助用户轻松理解整个项目的逻辑和实现方式。通过这些源码，用户可以学习到ASP的基本语法、服务器端脚本编写方法、数据库操作、用户权限管理等关键技术。 数据库设计文件：为了方便用户更好地理解系统的后台逻辑，每个项目中都附带了完整的数据库设计文件。这些文件通常包括数据库结构图、数据表设计文档，以及示例数据SQL脚本。用户可以通过这些文件快速搭建项目所需的数据库环境，并了解各个数据表之间的关系和作用。 详细的开发文档：每个资源包都附有详细的开发文档，文档内容包括项目背景介绍、功能模块说明、系统流程图、用户界面设计以及关键代码解析等。这些文档为用户提供了深入的学习材料，使得即便是从零开始的开发者也能逐步掌握项目开发的全过程。 项目演示与使用指南：为帮助用户更好地理解和使用这些ASP项目，每个资源包中都包含项目的演示文件和使用指南。演示文件通常以视频或图文形式展示项目的主要功能和操作流程，使用指南则详细说明了如何配置开发环境、部署项目以及常见问题的解决方法。 毕业设计参考：对于正在准备毕业设计的学生来说，这些资源包是绝佳的参考材料。每个项目不仅功能完善、结构清晰，还符合常见的毕业设计要求和标准。通过这些项目，学生可以学习到如何从零开始构建一个完整的Web系统，并积累丰富的项目经验。
基于SSM的高校实验室预约管理系统的设计与实现源码
08-10
本系统结合当前实验室预约、使用和管理中存在的问题和现状进行分析，设计一套可以解决实验室使用冲突、预约流程简单的实验室预约管理系统，满足实验室的在线预约申请与审核的基本需求，可以方便对实验室和实验室设备进行信息化管理。本系统中的教师可以查看当前实验室以及设备的基本信息，在线进行实验室的预约。管理员可以对实验室和实验室的设备进行管理，以及实验室的预约申请进行审核。
毕业设计javajsp超市会员管理系统ssm源码含文档工具包
08-10
毕业设计javajsp超市会员管理系统ssm源码含文档工具包 后台是ssm框架，页面是jsp，数据库mysql，jdk1.8，开发工具用ecplise、myecplise、sts、idea都可以 会员管理 公告育理 评论管理 前台 邮箱管理 积分管理 帮动 修改密码 退出系统 会员管理 公告管理 评论管理 工具管理 包含：源码、数据库脚本、论文、答辩ppt、开题报告、环境工具包、相同框架项目的安装教程（在说明文档中）
STM32F103C8T6 CAN通讯HAL库代码，采用中断接收，含Stm32CubeMX工程
最新发布
08-10
STM32F103C8T6 CAN通讯HAL库例程代码，采用中断接收，含Stm32CubeMX工程
c++的概要介绍与分析
08-10
C++是一种广泛使用的编程语言，它以其强大的功能和灵活性在计算机科学领域占据重要地位。以下是对C++资源的详细描述： 一、C++概述 C++是在C语言的基础上发展而来的，它保留了C语言的简洁、高效和接近硬件等特点，并增加了面向对象编程的支持。C++支持面向过程和面向对象的程序设计方法，使得开发者能够根据需要灵活选择编程范式。 二、C++特点 高效性：C++生成的代码质量高，程序执行效率高，仅比汇编语言慢10%～20%。 面向对象：支持封装、继承和多态等面向对象编程特性，便于构建复杂系统。 兼容性：与C语言高度兼容，绝大多数C语言程序可以不经修改在C++环境中运行。 可扩展性：支持泛型编程、模板等高级特性，提高了代码的可重用性和可维护性。 三、学习资源 官方文档：C++标准委员会发布的官方文档是了解C++语言特性的权威资源。 书籍：《C++ Primer》、《Effective C++》等经典书籍是学习C++的必备资料，它们详细介绍了C++的语法、特性和最佳实践。 在线课程：Coursera、网易云课堂等在线教育平台提供了丰富的C++课程，适合不同水平的学习者。 社区和论坛：Stac
AE文字动画预设，AE文字动画预设
08-10
AE文字动画预设
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值