矩阵快速幂小结&反思 不断补充ing

最新推荐文章于 2020-05-26 17:22:17 发布
最新推荐文章于 2020-05-26 17:22:17 发布
阅读量573 收藏
点赞数
分类专栏： ACM
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/YZBYZZ/article/details/12101631
版权
ACM 同时被 3 个专栏收录
35 篇文章 0 订阅
订阅专栏
hdu
14 篇文章 0 订阅
订阅专栏
poj
11 篇文章 0 订阅
订阅专栏

等整理完网络流的题目，再整理矩阵的，此坑待填

 

====================

 

poj3233 Matrix Power Series

 

#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

#define FD "%d"
#define FLL "%lld"

int matsz,mod;
const int N=70;

struct matrix
{
    int mat[N][N];
    matrix(int v)
    {
        memset(mat,0,sizeof(mat));
        for(int i=0;i<matsz;i++)
        mat[i][i]=v;
    }
    matrix operator *(matrix b)
    {
        matrix tmp(0);
        for(int i=0;i<matsz;i++)
        for(int k=0;k<matsz;k++)
        if(mat[i][k])
        {
          for(int j=0;j<matsz;j++)
          {
              tmp.mat[i][j]+=mat[i][k]*b.mat[k][j];
              tmp.mat[i][j]%=mod;
          }
        }
        return tmp;
    }
    void print()
    {
        for(int i=0;i<matsz;i++)
            for(int j=0;j<matsz;j++)
            printf(j==matsz-1?"%d\n":"%d ",mat[i][j]);
    }
};

matrix qpow(matrix t,int k)
{
    matrix ret(1);
    while(k)
    {
        if(k&1)
        ret=ret*t;
        t=t*t;
        k>>=1;
    }
    return ret;
}

int main()
{
    int n,k,m;
    while(~scanf("%d%d%d",&n,&k,&m))
    {
        matsz=n<<1;
        mod=m;
        matrix t(0);
        for(int i=0;i<n;i++)
        t.mat[i][i]=t.mat[i][i+n]=1;

        for(int i=n;i<matsz;i++)
        for(int j=n;j<matsz;j++)
        scanf("%d",&t.mat[i][j]);

        matrix ret=qpow(t,k+1);
        int ans;
        for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        t.mat[i][j]=t.mat[i+n][j+n];
        for(int i=n;i<matsz;i++)
        for(int j=0;j<n;j++)
        t.mat[i][j]=t.mat[i][j+n]*t.mat[i][j+n]%mod;

        for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        {
            if(i==j)ret.mat[i][j+n] = (ret.mat[i][j+n]-1+mod)%mod;
            printf(j==n-1?"%d\n":"%d ",ret.mat[i][j+n]);
        }
    }
    return 0;
}

 

 

 

 

 

 

loj1006  Hex-a-bonacci

 

题意：

          给出递推公式，求第n项。

做法：

         构造矩阵，然后快速幂。

 

#include <iostream>
#include <cstdio>
#include <cstring>
 
using namespace std;
const int Mod=10000007;
int a, b, c, d, e, f;
int fn( int n ) {
    if( n == 0 ) return a;
    if( n == 1 ) return b;
    if( n == 2 ) return c;
    if( n == 3 ) return d;
    if( n == 4 ) return e;
    if( n == 5 ) return f;
    return( fn(n-1) + fn(n-2) + fn(n-3) + fn(n-4) + fn(n-5) + fn(n-6) );
}
const int N=6;
struct matrix
{
   long long mat[N][N];
   matrix(){memset(mat,0,sizeof(mat));for(int i=0;i<N;i++)mat[0][i]=1;for(int i=1;i<N;i++)mat[i][i-1]=1;}
   matrix(int n){memset(mat,0,sizeof(mat));for(int i=0;i<N;i++)mat[i][i]=1;}
   void print(){for(int i=0;i<N;i++){for(int j=0;j<N;j++)cout<<mat[i][j]<<" ";printf("\n"); } }
   matrix operator *(matrix a)
   {
       matrix ret;
       for(int i=0;i<N;i++)
       for(int j=0;j<N;j++)
       {
           ret.mat[i][j]=0;
           for(int k=0;k<N;k++)
           ret.mat[i][j]+=mat[i][k]*a.mat[k][j];
           ret.mat[i][j]%=Mod;
       }
       return ret;
   }
};
 
int solve(int n)
{
    int ret=0;
    matrix a,t(n);
 
    n-=5;
    while(n)
    {
        if(n&1)
        t=t*a;
        a=a*a;
        n>>=1;
    }
    for(int i=0;i<N;i++)
    ret=(ret+t.mat[0][i]*fn(5-i))%Mod;
    return ret;
}
 
int main() {
//    freopen("in","r",stdin);
    int n, caseno = 0, cases,ret;
    scanf("%d", &cases);
    while( cases-- ) {
        scanf("%d %d %d %d %d %d %d", &a, &b, &c, &d, &e, &f, &n);
        if(n<=5)
        ret=fn(n);
        else
        ret=solve(n);
        printf("Case %d: %d\n", ++caseno,ret % 10000007);
    }
    return 0;
}

 

 

 

 

 

 

hdu4686 Arc of Dream（2013多校题目）

题意：

          给出递推公式，求第n项。

做法：

         构造矩阵，然后快速幂。

反思：

          注意及时取模，赋值不要写错了。

 

#include <cstdio>
#include <cstring>

using namespace std;

typedef long long LL;
#define FLL "%I64d"
const LL Mod=1e9+7;
const int Maxn=7;
struct martix
{
    LL mat[Maxn][Maxn];
    martix()
    {
        memset(mat,0,sizeof(mat));
    }
    martix(int a)
    {
        memset(mat,0,sizeof(mat));
        for(int i=0; i<Maxn; i++)mat[i][i]=1;
    }
    void operator *=(const martix &a)
    {
        martix tmp;
        for(int i=0; i<Maxn; i++)
          for(int k=0; k<Maxn; k++)
            if(mat[i][k])
            for(int j=0; j<Maxn; j++)
                {
                    tmp.mat[i][j]+=(mat[i][k]*a.mat[k][j])%Mod;
                    tmp.mat[i][j]%=Mod;
                }


        for(int i=0; i<Maxn; i++)
            for(int j=0; j<Maxn; j++)
                mat[i][j]=tmp.mat[i][j];
    }
};

martix qpow(martix a,LL n)
{
    martix ret(1);
    while(n)
    {
        if(n&1)
        ret*=a;
        a*=a;
        n>>=1;
    }
    return ret;
}

int main()
{
    LL n,a0,ax,ay,b0,bx,by,s1;

    while(~scanf(FLL,&n))
    {
        scanf(FLL FLL FLL FLL FLL FLL,&a0,&ax,&ay,&b0,&bx,&by);
        if(!n)
        {
            puts("0");
            continue;
        }
        s1=a0*b0%Mod;
        if(n==1)
        {
            printf(FLL,s1);
            puts("");
            continue;
        }

        martix ret;
        ret.mat[0][0]=ret.mat[0][1]=1;
        ret.mat[1][1]=ax*bx%Mod;
        ret.mat[1][2]=ax*by%Mod;
        ret.mat[1][3]=ay*bx%Mod;
        ret.mat[1][4]=1;

        ret.mat[2][2]=ax;
        ret.mat[2][5]=1;//!!!

        ret.mat[3][3]=bx;
        ret.mat[3][6]=ret.mat[4][4]=ret.mat[5][5]=ret.mat[6][6]=1;

        ret=qpow(ret,n-1);
        LL ans=0,a1,b1;
        a1=(a0%Mod*ax%Mod+ay)%Mod;
        b1=(b0%Mod*bx%Mod+by)%Mod;
        LL tmp[Maxn]={s1,a1*b1%Mod,a1,b1,ay*by%Mod,ay,by};
        for(int i=0;i<Maxn;i++)
        {
            ans+=(tmp[i]*ret.mat[0][i])%Mod;
            ans%=Mod;
        }
        printf(FLL,ans);
        puts("");
    }
    return 0;
}

 

 

 

 

 

 

 

poj3735 Training little cats

题意：

          给出n个数，一开始都是0。

          有三种操作：g i：第i个数加一；e i：第i个数变成0；s i j：第i和第j个数交换。

          每次给你k个操作，然后重复这k个操作m次，求最后n个数分别是多少。

做法：

           把三种操作用矩阵乘法来表示

反思：

           没用long long wa了一次。

 

 

#include <cstdio>
#include <cstring>

using namespace std;

#define SET(a,x) memset(a,x,sizeof(a))
#define RE(i,n) for(long long i = 1; i <= n; i++)
#define _RE(i,n) for(long long i = 0; i <= n; i++)

long long n,m,k;
const long long MAX = 110;

struct martix {
    long long mat[MAX][MAX];
    long long a,b;
    martix(long long f=0) {
        SET(mat,0);
        _RE(i,n)mat[i][i]=f;
    }
    void operator *= (martix t) {
        martix ret;
        _RE(i,n)
        _RE(k,n)if(mat[i][k])
            _RE(j,n) {
            ret.mat[i][j] += mat[i][k] * t.mat[k][j];
//              ret.mat[i][j] %= MOD;
        }
        _RE(i,n)
        _RE(j,n)
        mat[i][j] = ret.mat[i][j];
    }
    void dis(){
    _RE(i,n){
        _RE(j,n)
        printf("%I64d ",mat[i][j]);
        puts("");
    }

    }
};

martix qpow(martix ret,long long m){
    martix ans(1);
   while(m){
    if(m & 1)
    ans *= ret;
    ret *= ret;
    m >>= 1;
   }
    return ans;
}

int main() {
    while(scanf("%I64d%I64d%I64d",&n,&m,&k),n||m||k) {
        char op[5];
        long long x,y;
        martix ret(1);
        while(k--) {

            scanf("%s%I64d",op,&x);
            martix cur(1);
            if(op[0] == 'g') {
                cur.mat[x][0] = 1;
            } else if(op[0] == 'e') {
                cur.mat[x][x] = 0;
            } else if(op[0] == 's') {
                scanf("%I64d",&y);
                cur.mat[x][x] = cur.mat[y][y] = 0;
                cur.mat[x][y] = cur.mat[y][x] = 1;
            }
            cur *=  ret;
            ret = cur;
        }
       if(m == 0){
        RE(i,n)
        printf(i == n?"0\n":"0 ");
        continue;
       }
       ret = qpow(ret,m);
        RE(i,n) {
            printf(i == n?"%I64d\n":"%I64d ",ret.mat[i][0]);
        }
    }
}

 

 

 

 

 

YZBYZZ
关注
专栏目录
【矩阵分析】矩阵幂级数 发散 条件 | 幂级数 与 解析函数 的关系 | 幂级数 收敛半径r 的求法 | 人工智能 面试题：请介绍一下自编码器（Autoencoder）及其在无监督学习中的应用。
追光者♂：记录、分享、总结、提升，现象级专栏《Python从入门到人工智能》作者，无惧黑暗，坚信曙光
11-18 1863
【矩阵分析】矩阵幂级数 发散 条件 | 幂级数 与 解析函数 的关系 | 幂级数 收敛半径r 的求法 | 人工智能 面试题：请介绍一下自编码器（Autoencoder）及其在无监督学习中的应用。
ACM算法整理(不断补充ing)
Griffin_0的专栏
04-12 1521
动态规划 1.背包问题  (1)01背包 voidZeroOnepack(int F[],int C[],int W[]){ FOR(i,1,n) DFR(v,V,C[i]) F[v]=max(F[v],F[v-C[i]]+W[i]);}//初始化时//若背包不一定装满F全初始化为0//若装满 F[0]=0 其它为-inf   (2)完全背包 voi
矩阵快速幂在ACM中的应用
weixin_33736649的博客
12-14 110
矩阵快速幂在ACM中的应用 16计算机2黄睿博 首发于个人博客http://www.cnblogs.com/BobHuang/ 作为一个acmer，矩阵在这个算法竞赛中还是蛮多的，一个优秀的算法可以影响到一个程序的运行速度的快慢，在算法竞赛中常常采用快速幂算法，因为有些递推式及有些问题都可以化为矩阵，所以矩阵快速幂也就有了意义。 首先引入快速幂： 我们要求a^b，那么其实b是可以拆成二进...
ACM-矩阵之快速幂
我了个去
10-10 1396
ACM-矩阵之快速幂
2016ACM-ICPC亚洲区沈阳站 C-Recursive sequence 矩阵快速幂
Katapeltes
10-31 679
递推式为F[n]=2*F[n-2]+F[n-1]+n^4.用矩阵快速幂求解:构造 { F[n-2],F[n-1],n^4,n^3,n^2,n^1,1 } * X = { F[n-1],F[n],(n+1)^4,(n+1)^3,(n+1)^2,(n+1)^1,1 }.可求得矩阵 X = { 0 2 0 0 0 0 0 } { 1 1 0 0 0 0 0 } { 0 1 1 0
ACM-数论-矩阵快速幂 解斐波那契数列Yet another Number Sequence（UVA）
CS_W98的博客
12-21 445
题意：就是求某斐波那契数列的第n个数值，同时要进行取余运算，避免爆数据。value of n ranges in [0, 1000000000] value of m ranges in [1, 4]input4 0 1 11 3 0 1 42 4 0 1 22 4 0 1 21 4output89 4296 7711 946矩阵快速幂的原理 详情点链接这里介绍比较全，我主要讲一下这题的我踩
ACM-数论-矩阵快速幂 HDU5015&HDU2276
CS_W98的博客
05-05 227
HDU2276题目链接 HDU5015题目链接 2276 矩阵：（只有左边和本身有关联（1），其他都无关（0）） 第一个样例1 0 0 0 0 0 1 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 1代码：#include<iostream> #include<string.
一年级语文上册 ang eng ing ong 1教学反思 北京版 教案.doc
10-21
这篇文档主要讨论了一年级语文教学中关于汉语拼音后鼻韵母"ang eng ing ong"的教学反思和教案设计。在教学实践中，教师注重学生主体地位的体现，采取多种阅读方式如指名读、开火车读、同桌读，旨在激发学生参与度和...
小白积累本：stata矩阵计算初级操作的汇总（更新ing）
热门推荐
Hannahzyh的博客
05-26 1万+
作者是金融工程学生，正在学习stata，主要课程任务是通过矩阵运算实现一些stata内置算法包，所以积累了一些矩阵运算的常规操作。 刚接触时最大的困难是不懂这个软件的语言习惯，所以用备忘录记下了自己曾经困惑的地方，一些操作很小白，本学期会持续更新~ 写在前面的注意事项： 1、引用格式 `i’ $b 2、矩阵不显示为变量 3、直接输入Reg遇到缺省值时，会直接删除这一行。 设置系统能容纳的矩阵大小（有时stata不能容纳已有的矩阵，无法运算） set matsize 11000 诸如此类还有 set
新人教统编版一年级上册语文 汉语拼音13ɑng eng ing ong 课堂教学反思.docx
10-12
本文件是针对"新人教统编版一年级上册语文"的一堂汉语拼音教学课程，具体涉及的是第13课的ɑng、eng、ing、ong四个韵母的教学反思。教师在课堂中尝试通过创造有趣的情境和多样化的教学方式来提高一年级学生的学习...
ACM-数论-矩阵快速幂 HDU6030
CS_W98的博客
04-29 325
这里是题面 r:red b:blue 【题意】 有一串珍珠，长度为n（1e18） 每个珍珠要不染色成红色，要不染色成蓝色。 要求任何连续素数长度的珍珠，都必须是红色个数&amp;gt;=蓝色个数 让你求出有多少种对这串珍珠的染色方案。 一开始以为要打质数表，但是后来发现不需要，套路特别像斐波那契数列。 要求任何连续素数长度的珍珠，都必须是红色个数&amp;gt;=蓝色个数 1.用f[...
【ACM】矩阵快速幂+快速幂+费马小定理
小小小草儿的博客哟
11-08 1146
一：矩阵快速幂算法 矩阵快速幂的思想和数的快速幂的思想是一样的，但是需要自己实现矩阵的乘法，然后套用数的快速幂模板即可。 核心： 难点在于构造矩阵，一般用于可以推出递推公式的题目，发现时间复杂度为O（n），因此可以构造一个矩阵，利用矩阵快速幂算法把时间复杂度降低到O（logn）。 矩阵的快速幂算法是用于高效的计算矩阵的高次方的。比如矩阵A*A*A*A*A*A可以变为（A*A）（A*A）（A*...
hdoj 1575 Tr A 51nod oj 1113 矩阵快速幂【矩阵快速幂】
leibniz_zhang的博客
08-12 430
传送门：1575 矩阵快速幂--与一般快速幂一样-- 代码： #include #include #include using namespace std; struct node{ int n; long long shu[12][12]; node operator * (node &b){ node c; c.n=n; for (int i=1;i<=
hdu 5411 CRB and Puzzle（矩阵快速幂）
小G的ACM之路
08-21 604
题意： 给一个有向图，从任意点开始，最多走m步，求形成的图案总数。 思路： 令dp[i][j]dp[i][j]表示：走jj步最后到达i的方法数， 则dp[i][j]=∑dp[k][j−1]dp[i][j]=\sum{dp[k][j-1]}，其中k表示：走jj步可以直接到达ii的点。 ans=∑dp[i][j]ans=\sum{dp[i][j]}。 此题的关键在于，如何减少状
矩阵快速幂基础讲解
weixin_30408739的博客
06-03 402
1.基础知识储备篇 矩阵的相关运算会再线性代数中学到。 1.1矩阵的定义：  N阶方阵（N阶矩阵）：行数m与列数n相同的矩阵，如下图所示就是一个44的方阵：* 行矩阵（行向量）：只有一行的矩阵，下图就是一个行矩阵： ​ 列矩阵（列向量）：只有一列的矩阵，下图就是一个列矩阵： ​ 同型矩阵：设先有矩阵A和矩阵B，矩阵A的行数与列数和矩阵B的相同，则矩阵A、B是同型矩阵。 ​ 单位矩阵：在矩...
hdoj 1575 Tr A (矩阵快速幂)
R
11-13 1874
Tr A Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 4549    Accepted Submission(s): 3428 Problem Description A为一个方阵，则Tr A表示A的迹（就是主
矩阵 快速幂
weixin_30407099的博客
11-29 820
矩阵的快速幂是用来高效地计算矩阵的高次方的。将朴素的o（n）的时间复杂度，降到log（n）。 这里先对原理（主要运用了矩阵乘法的结合律）做下简单形象的介绍： 一般一个矩阵的n次方，我们会通过连乘n-1次来得到它的n次幂。 但做下简单的改进就能减少连乘的次数，方法如下： 把n个矩阵进行两两分组，比如：A*A*A*A*A*A => (A*A)*(A*A)*(A*A) 这...
poj 1002 487-3279 （（含易错点、测试数据） 简单模拟 ）
程序人生
04-08 3870
题目大意： 把大写字母按照对应关系转换成数字，统计每个电话号码出现的次数，按号码升序输出那些有重复的。 思路： 用字符串数组接收输入，转换成七个数字的同时表示成一个七位数，方便比较。 全部转换完后，进行排序。 然后遍历一次，边统计次数，边输出。 思路应该是比较清晰的，提交后却是wa。 然后重新看了一下，发现自己一开始看题的时候有注意到前导零的控制以及“No duplicates“的情况，但写的时候...
scau 1138 代码等式 （ 并查集 ）
程序人生
09-28 2524
  1138 代码等式   时间限制:500MS  内存限制:65536K 提交次数:59 通过次数:21 题型: 编程题   语言: 无限制   Description 一个代码等式就是形如x1x2...xi=y1y2...yj，这里xi和yj是二进制的数字（0或1）或者是一个变量（如英语中的小写字母）。每一个变量都是一个有固定长度的二进制代码。例如： a,b,c,d,e是变且它...
ING家族基因在结直肠癌中的低表达与临床关系
然而，关于ING1至ING5这五个家族成员在结直肠癌中的表达情况及其与临床特征的关联，尚不明确。 研究表明，ING1-5在结直肠癌组织中的mRNA表达水平显著降低。这种低表达状态与肿瘤的一些关键临床参数紧密相关。具体来...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值