矩阵和线性方程组
HAI__嗨I起来
IT
展开
-
UVA 10870 递推关系 矩阵快速幂
题意:f(n) = a1 *f(n - 1) + a2 *f(n - 2) + a3 *f(n - 3) + … + ad* f(n - d), n > d.求f(n)见白书 155页。 由于n太大,不能直接递推,需要用矩阵快速幂来解决,时间复杂度为O(d^3logn) 举例,d=5的矩阵关系式为: |a1 a2 a3 a4 a5|原创 2016-08-15 19:48:58 · 620 阅读 · 0 评论 -
HDU 5833 高斯消元 异或方程组
题意:给出n个正整数,从中选出1个或者多个,使得选出来的整数乘积是完全平方数,一共有多少种选法。思路:用01向量表示一个数,再用n个01向量来表示我们的选择,因为完全平方数要求素因子的次数一定要是偶数的,所以我们可以统计的将奇数当作1,偶数当作0,那么就是一组可以变换成oxr的方程组,最后的结果有自由变量f个,答案是2^f-1,f求解就是求n-方程组的秩,(本题不允许一个都不选,所以减原创 2016-08-16 10:02:11 · 359 阅读 · 0 评论 -
HDU 5667 矩阵快速幂关于指数的递推
SequenceTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 1498 Accepted Submission(s): 505Problem Description Holion Augus原创 2016-08-16 15:58:08 · 261 阅读 · 0 评论 -
POJ 1166 枚举或者高斯消元
给出9个钟表的状态,给出九种操作,问最少要操作几次能把所有的钟表调回12点。思路:对于9个钟表分别列方程,然后高斯消元即可。然后用每个未知量表示是否进行此操作。所以取值是1或者0接下来就是枚举每个未知量得系数使得每个等式都成立。#include#include#include#include#define LL long long#define bug p原创 2016-08-17 21:37:52 · 262 阅读 · 0 评论 -
POJ 1222 开关问题高斯消元法
一个01矩阵,表示灯的亮灭状态,每次操作可以改变一个十字形状内的五个灯的状态。问能否将所有灯熄灭。最后输出开关的状态。思路:开灯关灯问题,5*6的灯阵,将每一个位置上开关状态看做一个变元,30个变元,而对于30个位置的灯的状态当成开关状态的解,列出30个异或方程,高斯消元解方程即可(增广矩阵)。对于每个状态,影响他的条件与五个位置的开关状态有关,让这五个状态的开关系数为1 ,其原创 2016-08-24 17:01:25 · 482 阅读 · 0 评论 -
POJ 2947 (高斯消元解同模方程)
公司被吞并,老员工几乎全部被炒鱿鱼。一共有n种不同的工具,编号1-N(代码中是0—N-1), 每种工具的加工时间为3—9天 ,但是现在老员工不在我们不知道每种工具的加工时间,庆幸的是还保留着一些对工人制造工具的记录,对于每个老员工,他的记录包括,他开始工作的时间(在某个星期的星期几),被炒鱿鱼的时间(某个星期的星期几),在第几个星期不知道.....在这段时间里,他正好加工了k件物品,给出了这k件物原创 2016-08-24 22:04:06 · 310 阅读 · 0 评论 -
POJ 2891 Strange Way to Express Integers解线性同余方程组(中国剩余定理不互质版)
题意:给出k个模方程组:x mod ai = ri。求x的最小正值。如果不存在这样的x,那么输出-1.题解:由于这道题目里面的ai、ri之间不满足两两互质的性质,所以不能用中国剩余定理直接求解。不过,我们可以模仿中国剩余定理的做法来解决这个问题。如果只有一个方程:x mod a0 = r0。那么,显然x的最小正值为a0+r0。根据模的性质,我们容易得知,x+a0*k原创 2016-08-25 19:35:11 · 624 阅读 · 0 评论 -
POJ 3101 Astronomy 轨道相遇问题,求n个分数的最小公倍数
题意:给你一个n,然后给你n个星球的周期。让你求出经过多长时间可以使所有的星球可以在条直线上。思路: 求轨道相遇问题,设经过了t个时间星球A转的角度为 (2π/T1)*t 星球B转的角度为(2π/T2)*t 他们在一条直线上所以可得:(2π/T1)*t-(2π/T2)*t=n*π,半圈就相遇,所以n取1。-->t=T1*T2/(2T1-2T2);然后需要求出所有的星球与第原创 2016-08-26 08:56:38 · 572 阅读 · 0 评论