ACM 矩阵题目整理

先从最基础的矩阵快速幂加速递推开始。

HDU 1005 Number Sequence

|f[n-2],f[n-1]|* |0 B| =|f[n-1], B*f[n-2]+A*f[n-1]|=|f[n-1],f[n]|

　　　　　　　 |1 A|

建立矩阵如上然后利用快速幂求解即可。答案是mat[0][1]。

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
struct Matrix
{
    int n,m;
    int a[MAXN][MAXN];
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<b.m; ++j)
                ; k<m; ++k)
                {
                    tmp.a[i][j]+=a[i][k]*b.a[k][j];
                    tmp.a[i][j]%=;
                }
        return tmp;
    }
};
Matrix quickPow(Matrix mat,int n)
{
    Matrix res;
    res.clear(,);
    res.a[][]=res.a[][]=;
    while(n)
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    int A,B,n;
    while(scanf("%d%d%d",&A,&B,&n)!=EOF)
    {
        if(!A&&!B&&!n) break;
        ||n==)
        {
            puts(");
            continue;
        }
        Matrix mat;
        mat.clear(,);
        mat.a[][]=;
        mat.a[][]=B;
        mat.a[][]=;
        mat.a[][]=A;
        Matrix f;
        f.clear(,);
        f.a[][]=;
        f.a[][]=;
        n-=;
        Matrix s=f*quickPow(mat,n);
        printf(][])%);
    }
    ;
}

HDU 4920 Matrix multiplication

矩阵乘法，直接算会超时。由于是模3，所以0很多，也就是稀疏矩阵。可以优化。

参见下文：http://www.cnblogs.com/jackiesteed/articles/2021604.html

#include<iostream>
#include<vector>
#include<map>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<string>
#define LL long long
using namespace std;
][],B[][];
][];
int n;
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        ; i<=n; ++i)
            ; j<=n; ++j)
            {
                scanf("%d",&A[i][j]);
                A[i][j]%=;
            }
        ; i<=n; ++i)
            ; j<=n; ++j)
            {
                scanf("%d",&B[i][j]);
                B[i][j]%=;
            }
        memset(C,,sizeof(C));
        ; i<=n; ++i)
            ; j<=n; ++j)
            {
                ) continue;
                ; k<=n; ++k)
                    C[i][k]+=A[i][j]*B[j][k];
            }
        ; i<=n; ++i)
        {
            ; j<=n; ++j)
                ) printf();
                );
            printf("\n");
        }
    }
    ;
}

HDU 1575 Tr A

求矩阵的迹。直接快速幂即可。

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
struct Matrix
{
    int n,m;
    int a[MAXN][MAXN];
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<b.m; ++j)
            {
                ; k<b.m; ++k)
                {
                    ==a[i][j]) continue;//优化！
                    tmp.a[i][k]+=a[i][j]*b.a[j][k];
                    tmp.a[i][k]%=;
                }
            }
        return tmp;
    }
    void setOne(int x)
    {
        clear(x,x);
        ; i<x; ++i)
            a[i][i]=;
    }
};
Matrix one;
Matrix quickPow(Matrix mat,int n)
{
    Matrix res=one;
    while(n)
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        int n,k;
        scanf("%d%d",&n,&k);
        one.setOne(n);
        Matrix mat;
        mat.clear(n,n);
        ;i<n;++i)
            ;j<n;++j)
            scanf("%d",&mat.a[i][j]);
        mat=quickPow(mat,k);
        ;
        ;i<n;++i)
        {
            ans+=mat.a[i][i];
            ans%=;
        }
        printf("%d\n",ans);
    }
    ;
}

HDU 1757 A Simple Math Problem

一个简单的递推，建立矩阵后，快速幂。

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
int M;
struct Matrix
{
    int n,m;
    int a[MAXN][MAXN];
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<b.m; ++j)
            {
                ; k<b.m; ++k)
                {
                    ==a[i][j]) continue;//优化！
                    tmp.a[i][k]+=a[i][j]*b.a[j][k];
                    tmp.a[i][k]%=M;
                }
            }
        return tmp;
    }
    void setOne(int x)
    {
        clear(x,x);
        ; i<x; ++i)
            a[i][i]=;
    }
};
Matrix one;
Matrix quickPow(Matrix mat,int n)
{
    Matrix res=one;
    while(n)
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    int k;
    one.setOne();
    while(scanf("%d%d",&k,&M)!=EOF)
    {
        Matrix mat;
        mat.clear(,);
        ; i<; ++i)
            mat.a[i+][i]=;
        ; i<; ++i)
            scanf(-i][]);
        ) printf("%d\n",k%M);
        else
        {
            Matrix f;
            f.clear(,);
            ; i<; ++i)
                f.a[][i]=i;
            k-=;
            f=f*quickPow(mat,k);
            printf(][]);
        }
    }
    ;
}

UVa 10870 Recurrences

与上面类型相似的题，建立矩阵，再快速幂加速。注意会超int。

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
int M;
struct Matrix
{
    int n,m;
    LL a[MAXN][MAXN];
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<b.m; ++j)
            {
                ; k<b.m; ++k)
                {
                    ==a[i][j]) continue;//优化！
                    tmp.a[i][k]+=a[i][j]*b.a[j][k];
                    tmp.a[i][k]%=M;
                }
            }
        return tmp;
    }
    void setOne(int x)
    {
        clear(x,x);
        ; i<x; ++i)
            a[i][i]=;
    }
};
Matrix one;
Matrix quickPow(Matrix mat,LL n)
{
    Matrix res=one;
    while(n)
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    int d,n;
    while(scanf("%d%d%d",&d,&n,&M)!=EOF)
    {
        if(!d&&!n&&!M) break;
        one.setOne(d);
        Matrix mat;
        mat.clear(d,d);
        ; i<d; ++i)
            mat.a[i+][i]=;
        ; i<d; ++i)
            scanf(-i][d-]);
        Matrix f;
        f.clear(,d);
        ; i<d; ++i)
            scanf(][i]);
        if(n<=d)
            printf(][n-]%M);
        else
        {
            n-=d;
            f=f*quickPow(mat,n);
            printf(][d-]%M);
        }
    }
    ;
}

POJ 3734 Blocks

递推+矩阵快速幂。

设计状态的时候要考虑所有可能的情况，保证无重复无遗漏，再考虑转移。最后用矩阵加速。

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
;
struct Matrix
{
    int n,m;
    int a[MAXN][MAXN];
    Matrix(,):n(x),m(y)
    {
        memset(a,,sizeof(a));
    }
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<b.m; ++j)
            {
                ; k<b.m; ++k)
                {
                    ==a[i][j]) continue;//ÓÅ»¯£¡
                    tmp.a[i][k]+=a[i][j]*b.a[j][k];
                    tmp.a[i][k]%=M;
                }
            }
        return tmp;
    }
    void setOne(int x)
    {
        clear(x,x);
        ; i<x; ++i)
            a[i][i]=;
    }
};
Matrix one;
Matrix quickPow(Matrix mat,int n)
{
    Matrix res=one;
    while(n)
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    one.setOne();
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        Matrix mat(,);
        mat.a[][]=,mat.a[][]=,mat.a[][]=;
        mat.a[][]=,mat.a[][]=,mat.a[][]=;
        mat.a[][]=,mat.a[][]=,mat.a[][]=;
        Matrix f(,);
        f.a[][]=;
        Matrix ans=f*quickPow(mat,n);
        printf(][]);
    }
    ;
}

POJ 3420 Quad Tiling

这个题的关键是在求递推式。

设dp[i]表示填满前i个方块的种数。那么从dp[i-1]转移有1种情况，从dp[i-2]转移有4种情况，在之后从dp[i-a]，如果a是奇数有2种情况，a是偶数有3种情况。

这样最终化简得dp[i]=dp[i-1]+5dp[i-2]+dp[i-3]-dp[i-4]

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
int M;
struct Matrix
{
    int n,m;
    int a[MAXN][MAXN];
    Matrix(,):n(x),m(y)
    {
        memset(a,,sizeof(a));
    }
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<b.m; ++j)
            {
                ; k<b.m; ++k)
                {
                    ==a[i][j]) continue;//优化
                    tmp.a[i][k]+=a[i][j]*b.a[j][k];
                    tmp.a[i][k]%=M;
                }
            }
        return tmp;
    }
    void setOne(int x)
    {
        clear(x,x);
        ; i<x; ++i)
            a[i][i]=;
    }
};
Matrix one;
Matrix quickPow(Matrix mat,int n)
{
    Matrix res=one;
    )
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    int n;
    one.setOne();
    while(scanf("%d%d",&n,&M)!=EOF)
    {
        if(!n&&!M) break;
        Matrix mat(,);
        ; i<; ++i)
            mat.a[i+][i]=;
        mat.a[][]=-,mat.a[][]=,mat.a[][]=,mat.a[][]=;
        Matrix f(,);
        f.a[][]=,f.a[][]=,f.a[][]=,f.a[][]=;
        ) printf(][n]%M);
        else
        {
            n-=;
            Matrix s=f*quickPow(mat,n);
            printf(][]);
        }
    }
    ;
}

ZOJ 3690 Choosing number

a[i]表示第i个人选小于等于k的数字的合法种数，a[i]=(k-1)*a[i]+k*b[i]

b[i]表示第i个人选大于等于k的数字的合法种数，b[i]=(m-k)*a[i]+(m-k)*b[i]

a[0]=0,b[0]=1,答案是a[n]+b[n]

建立矩阵快速幂计算即可。

#include<iostream>
#include<vector>
#include<cmath>
#include<map>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#include<cmath>
#include<set>
#include<string>
#define INF 1000000005
#define LL long long
using namespace std;
;
;
struct Matrix
{
    int n,m;
    LL a[MAXN][MAXN];
    Matrix(,):n(x),m(y)
    {
        memset(a,,sizeof(a));
    }
    ,)
    {
        n=x;
        m=y;
        memset(a,,sizeof(a));
    }
    Matrix operator * (const Matrix &b)const
    {
        Matrix tmp;
        tmp.clear(n,b.m);
        ; i<n; ++i)
            ; j<m; ++j)
            {
                ; k<b.m; ++k)
                {
                    ==a[i][j]) continue;//优化！
                    tmp.a[i][k]+=a[i][j]*b.a[j][k];
                    tmp.a[i][k]%=M;
                }
            }
        return tmp;
    }
    void setOne(int x)
    {
        clear(x,x);
        ; i<x; ++i)
            a[i][i]=;
    }
};
Matrix one;
Matrix quickPow(Matrix mat,int n)
{
    Matrix res=one;
    while(n)
    {
        ) res=res*mat;
        mat=mat*mat;
        n=n>>;
    }
    return res;
}
int main()
{
    int n,m,k;
    one.setOne();
    while(scanf("%d%d%d",&n,&m,&k)!=EOF)
    {
        Matrix mat(,);
        mat.a[][]=k-,mat.a[][]=k;
        mat.a[][]=m-k,mat.a[][]=m-k;
        Matrix f(,);
        f.a[][]=,f.a[][]=;
        Matrix s=quickPow(mat,n)*f;
        printf(][]+s.a[][])%M);
    }
    ;
}

ACM 矩阵题目整理的更多相关文章

ACM 字符串题目整理
AC自动机 UVa 11468 Substring AC自动机+概率DP. 注意要补全不存在的边. 为什么要补全不存在的边呢?补全以后可以直接找到状态的转移,即从所有子节点就可以实现所有状态转移. ...
Noip往年题目整理
Noip往年题目整理张炳琪一．历年题目按时间倒序排序年份 T1知识点 T2知识点 T3知识点得分总体 2016day1 模拟 Lca,树上差分期望dp 144 挺难的一套题目,偏思维难度 ...
NOIp初赛题目整理
NOIp初赛题目整理这个 blog 用来整理扶苏准备第一轮 csp 时所做的与 csp 没有关系的历年 noip-J/S 初赛题目,记录了一些我从不知道的细碎知识点,还有一些憨憨题目,不定期 ...
ACM 暴力搜索题题目整理
UVa 129 Krypton Factor 注意输出格式,比较坑爹. 每次要进行处理去掉容易的串,统计困难串的个数. #include<iostream> #include<vec ...
ACM - 动态规划专题题目整理
CodeForces 429B Working out 预处理出从四个顶点到某个位置的最大权值,再枚举相遇点,相遇的时候只有两种情况,取最优解即可. #include<iostream> ...
BZOJ 题目整理
bzoj 500题纪念总结一发题目吧,挑几道题整理一下,(方便拖板子) 1039:每条线段与前一条线段之间的长度的比例和夹角不会因平移.旋转.放缩而改变,所以将每条轨迹改为比例和夹角的序列,复制一份 ...
第八届蓝桥杯c/c++省赛题目整理
第一题标题: 购物单小明刚刚找到工作,老板人很好,只是老板夫人很爱购物.老板忙的时候经常让小明帮忙到商场代为购物.小明很厌烦,但又不好推辞. 这不,XX大促销又来了!老板夫人开出了长长的购物单,都 ...
ACM金牌选手整理的【LeetCode刷题顺序】
算法和数据结构知识点图首先,了解算法和数据结构有哪些知识点,在后面的学习中有大局观,对学习和刷题十分有帮助. 下面是我花了一天时间花的算法和数据结构的知识结构,大家可以看看. 后面是为大家精心挑 ...
内工大acm校赛--整理代码
题目:小明搜到一行无缩进无换行代码,请帮小明整理代码.无for语句和case语句,而且只有一个主函数.你只要控制注意“:”“{”“}”这三个符号带来的缩进和换行效果就行. Input: 输入只有一行, ...

随机推荐

【原创】js中利用cookie实现记住密码功能
在登录界面添加记住密码功能,我首先想到的是在java后台中调用cookie存放账号密码,大致如下: HttpServletRequest request HttpServletResponse res ...
如何使用百度音乐搜索接口API
百度有开放音乐搜索的api 比如: http://box.zhangmen.baidu.com/x?op=12&count=1&title=大约在冬季$$齐秦$$$$ http://b ...
再谈Bellman-Ford
这几天学校女生节,挺累的,感觉还是挺好玩的,前几天看了一下最短路,Bellman-fort算法果然比较厉害,今天又参考了刘汝佳的两本书,有了一点新的认识. 废话不说,先上代码: #include &l ...
QT笔记之VS开发添加类
1. 2. 3.
C语言面试题（三）
这篇主要聚焦在排序算法,包括常见的选择排序,插入排序,冒泡排序,快速排序.会对这四种排序的时间复杂度和空间复杂度进行探究. a.选择排序 int main(int argc,char **argv){ ...
【树莓派】基于TinyProxy搭建HTTP代理服务器
一.前言关于为什么要玩玩HTTP代理就不用我多说了. 二.搭建环境 * Linux laptop 2.6.32-45-generic #100-Ubuntu SMP Wed Nov 14 10:4 ...
implicit operator
class Digit { public Digit(double d) { val = d; } public double val; // ...other members // User-def ...
mysql命令导出导入数据和结构
在命令行下mysql的数据导出有个很好用命令mysqldump,它的参数有一大把,可以这样查看: mysqldump 最常用的: mysqldump -uroot -pmysql databasefo ...
zoj3228Searching the String（ac自动机)
链接这个题把病毒分为了两种,一种包含可以覆盖,另一种不可以,需要分别求出包含他们的个数,可以把两种都建在一颗tire树上,在最后求得时候判断一下当前节点是属于哪种字符串,如果是不包含的需要判断一下p ...
STM32学习笔记(四) RCC外设的学习和理解
RCC时钟模块并不好理解,初次接触我也是一头雾水,而且我真正掌握它的时候也比较晚,是我在学习uC/os-II,需要分析时钟时才有了深刻认识.但在学习中我却一定要把放在了前列,因为这是整个嵌入式最重要的 ...

ACM 矩阵题目整理

ACM 矩阵题目整理的更多相关文章

随机推荐

热门专题