整除
同余
重要的性质在于:
只要运算过程只有乘法和加分,那么对每一步求余和对最后结果求余的余数都是一样的
GCD(最大公因子)
LCM(最小公倍数)
例题分析
状态分析
首先枚举所有b1的因子,然后判断是否满足条件
代码复现
def gcd(a, b):
if b == 0:
return a
return gcd(b,a%b)
def lcm(a, b):
return a*b//gcd(a,b)
T = int(input())
for _ in range(T):
ans = 0
a0,a1,b0,b1 = map(int,input().split())
#枚举所有b1的因子
for x in range(1, int(b1**0.5)+1):
if b1 % x == 0:
if gcd(x,a0) ==a1 and lcm(x, b0) == b1:
ans += 1
y = b1//x
if x != y:
if gcd(y, a0) == a1 and lcm(y, b0) == b1:
ans += 1
print(ans)
高斯消元
如果不知道高斯消元法是什么可以去看b站的【五分钟直觉理解高斯消元法(3)】 https://www.bilibili.com/video/BV1r5411J7YR/?share_source=copy_web&vd_source=fe59d667f9424d841d598d939d8c404a
通过高斯消元得到的阶梯形矩阵有无数个,但是最简的阶梯形矩阵只有一个,通过判断系数与常数的关系可以分为以下三种情况:
模拟高斯消元代码
from bitarray import bitarray
N = 2800
a = [bitarray(N+1) for _ in range(N+1)] # 初始化一个2800*2800的位数组,用于存储方程组的系数矩阵
n, m, x, ans = 0, 0, 0, 0 # 初始化变量,n是方程数,m是当前处理的方程的非零元素数,x用于读入数据,ans用于存放结果
def gauss():
cnt = 0 # 记录当前处理到第几个方程
for i in range(1, n+1):
maxx = cnt + 1
for j in range(i+1, n+1): # 找到第i列中系数最大的行
if a[j][i] and not a[maxx][i]:
maxx = j
a[cnt+1], a[maxx] = a[maxx], a[cnt+1] # 交换当前行和最大行
if not a[i][i]: # 如果当前行的主元是0,则跳过
continue
cnt += 1
for j in range(1, n+1): # 对每一行进行消元处理
if a[j][i] and i != j:
a[j] = a[j] ^ a[i] # 使用异或操作进行行变换,实现消元
if cnt < n: # 如果有方程组中有方程无解,则输出无解
for i in range(1, n+1):
if not a[i][i] and a[i][n+1]:
print("no solution")
# 输入处理部分
n = int(input()) # 读入方程数
for i in range(1, n+1):
m = int(input()) # 读入每个方程的非零元素数
if m % 2 != 0: # 如果非零元素个数为奇数,设置增广矩阵的对应项为1
a[i][n+1] = 1
a[i][i] = 1
while m > 0: # 读入每个非零元素的位置
x = int(input())
a[i][x] = 1 # 设置系数矩阵中对应的项为1
m -= 1
gauss() # 调用高斯消元法函数
感悟
坚持每天学一点点吧,积少成多,厚积薄发!
蓝桥杯云课学习笔记分享,欢迎大佬们批评指正!
一直在进步就好咯!
by 闻不多