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: 输入只有一行, ...
随机推荐
- 搭建自己本地yum源
1.挂载系统光盘(注:medi下的cdrom是我自己创建的,可以挂载在任意目录) [root@liutao ~]# mount /dev/cdrom /media/cdrom/ 2.修改yum配置文件 ...
- (6) 深入理解Java Class文件格式(五)
前情回顾 本专栏的前几篇博文, 对class文件中的常量池进行了详细的解释. 前文讲解了常量池中的7种数据项, 它们分别是: CONSTANT_Utf8_info CONSTANT_NameAndTy ...
- SQL MD5加密
) 加密结果:
- HDU 5795 A Simple Nim（简单Nim）
p.MsoNormal { margin: 0pt; margin-bottom: .0001pt; text-align: justify; font-family: Calibri; font-s ...
- 忙了好一阵，今天随便写篇关于canvas的小东西
前几天在朋友圈发了几条3D demo视频,其中就有3D空间组成各种图形.如上! 那么这些图形的每个mesh的坐标可不是手动去写,如果你愿意我当然不拦着!所以今天这篇就来介绍如何获得这些图形的坐标数据. ...
- 《BI项目笔记》创建计算成员
也可以利用脚本来生成计算成员
- HTML基础（1）
1.前端开发语言介绍 HTML(Hypertext Markup Language)—— 结构—— 超文本标记语言 CSS(Cascading Style Sheets)—— 样式—— 层叠样式表 J ...
- 6个html5页面适配iphone6的技巧
iphone6及iphone6plus已经出来一段时间了.很多移动端网站,以前写死body为320px的,现在估计也忙着做适配了. 大屏幕手机其实一直有,只是以前大家没怎么重视,移动端的H5页面大部分 ...
- linux进程编程：子进程创建及执行函数简介
linux进程编程:子进程创建及执行函数简介 子进程创建及执行函数有三个: (1)fork();(2)exec();(3)system(); 下面分别做详细介绍.(1)fork() 函数定 ...
- #pragma once与 #ifndef的区别
为了避免同一个文件被include多次 1 #ifndef方式2 #pragma once方式 在能够支持这两种方式的编译器上,二者并没有太大的区别,但是两者仍然还是有一些细微的区别. ...