poj 1222 EXTENDED LIGHTS OUT 高斯消元

本文介绍了一种算法,用于解决通过点击灯来改变其自身及其周围灯的状态的问题。算法涉及矩阵运算、高斯消元等技术,适用于电子设备或游戏场景中灯的控制逻辑。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

传送门

题目就是给你一些灯, 点击一个灯, 就会改变它自身和它四周的所有的灯的状态。 给你所有灯的初始状态, 要你输出所有灯的点击情况, 1代表点击, 0代表不点击。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define mem(a) memset(a, 0, sizeof(a))
 4 const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
 5 int a[35][35], ans[35];
 6 void gauss() {
 7     for(int k = 0; k<30; k++) {
 8         int i = k;
 9         while(i<30&&!a[i][k])
10             i++;
11         if(i>k) {
12             for(int r = 0; r<=30; r++)
13                 swap(a[i][r], a[k][r]);
14         }
15         for(int i = 0; i<30; i++) {
16             if(i!=k && a[i][k]) {
17                 for(int j = 0; j<=30; j++) {
18                     a[i][j] ^= a[k][j];
19                 }
20             }
21         }
22     }
23     for(int i = 0; i<30; i++) {
24         if(a[i][30]) {
25             int j = 0;
26             while(j<30&&!a[i][j])
27                 j++;
28             ans[j] = a[i][30];
29         }
30     }
31 }
32 int main()
33 {
34     int t;
35     cin>>t;
36     for(int cnt = 1; cnt<=t; cnt++) {
37         mem(a);
38         mem(ans);
39         for(int i = 0; i<30; i++) {
40             scanf("%d", &a[i][30]);
41         }
42         for(int i = 0; i<30; i++)
43             a[i][i] = 1;
44         for(int i = 0; i<5; i++) {
45             for(int j = 0; j<6; j++) {
46                 for(int k = 0; k<4; k++) {
47                     int x = i+dir[k][0];
48                     int y = j+dir[k][1];
49                     if(x>=0&&x<5&&y>=0&&y<6) {
50                         a[i*6+j][x*6+y] = 1;
51                     }
52                 }
53             }
54         }
55         gauss();
56         printf("PUZZLE #%d\n", cnt);
57         for(int i = 0; i<30; i++) {
58             int x = i/6;
59             int y = i%6;
60             a[x][y] = ans[i];
61         }
62         for(int i = 0; i<5; i++) {
63             printf("%d", a[i][0]);
64             for(int j = 1; j<6; j++) {
65                 printf(" %d", a[i][j]);
66             }
67             puts("");
68         }
69     }
70 }

 

转载于:https://www.cnblogs.com/yohaha/p/5009691.html

### 矩阵在扩展版 Lights Out 游戏中的实现 #### 背景介绍 关灯游戏(Lights Out)是一种基于逻辑的益智游戏,其核心机制涉及通过一系列操作使所有灯光熄灭。该游戏可以通过线性代数的方法建模并求解,其中矩阵扮演了重要角色。具体来说,游戏的状态可以用二进制向量表示,而每次点击的操作可以看作是对当前状态施加的一个变换。 对于扩展版本的游戏(如 POJ 1222),通常是一个 \( M \times N \) 的网格,初始状态下某些位置可能亮起或熄灭。目标是找到一种最小化翻转次数的方式使得整个网格变为全零状态(即所有灯均关闭)。这一过程可通过构建一个对应的布尔方程组来描述,并利用高斯消法或其他数值技术求解[^3]。 #### 数学模型建立 假设我们有一个大小为 \( m=5, n=6 \) 的棋盘,则总共有 30 个独立变量代表各个格子是否被按压过。设这些未知数构成列向量 \( x=[x_1,x_2,...,x_{30}]^T \),如果某位 i 对应的位置需要按下一次就令 xi=1 否则等于 0 。接着定义另一个长度相同的响应矢量 y ,它记录的是最终期望达到的目标配置 —— 这里就是全是 'off' 或者说都是 0 值的情况下的布局图样 [y₁,y₂,…,yn]=₀₆₀⁰[^4]. 为了表达上述关系,我们需要创建一个系数矩阵 A 来捕捉每一个按钮动作如何影响周围邻居节点的变化情况: \[ Ax = b \] 这里, - **A** 是一个稀疏矩阵,每一行对应于某个特定单格及其邻域的影响模式; - **b** 表示初始条件下的光源分布状况; - 解决这个系统意味着寻找合适的输入序列 x 以满足给定的要求 b. 由于涉及到异或运算特性(xor operation), 实际上是在有限字段 GF(2)=Z/2Z 上执行计算而不是普通的实数空间 R^n 中进行处理. 因此标准算法比如 Gaussian elimination 需要做适当调整以便适应这种特殊的算术环境. #### 使用高斯消去法解决问题 一旦建立了这样的数学框架之后,就可以采用经典的 Gauss-Jordan Elimination 方法逐步简化增广矩阵直到得到唯一可行解或者是证明无解存在为止。特别注意当面对大规模实例时效率问题变得至关重要因此也可能考虑其他更高效的替代策略例如 Bitmasking Techniques 结合快速傅立叶变换 FFT 加速乘法步骤等等优化手段提升性能表现水平. 以下是用 Python 编写的简单例子展示如何运用 NumPy 库完成基本功能演示: ```python import numpy as np def create_matrix(m,n): """Create coefficient matrix for lights out game.""" size=m*n mat=np.zeros((size,size)) # Fill diagonal elements and neighbors based on grid structure. for r in range(m): for c in range(n): idx=r*n+c # Current cell itself always toggles. mat[idx,idx]=1 # Toggle top neighbor if exists. if r>0: mat[(r-1)*n+c,idx]=mat[idx,(r-1)*n+c]=1 # Bottom neighbor similarly handled... ... return mat.astype(int) # Example usage creating small test case setup. if __name__=="__main__": M,N=(5,6) initial_state=[[int(c)for c in line.strip().split()]for _in range(M)] target_vector=sum(initial_state,[]) coeff_mat=create_matrix(M,N) augmented=np.hstack([coeff_mat,np.array(target_vector).reshape(-1,1)]) rank_before_elim=np.linalg.matrix_rank(coeff_mat) reduced_form,rref_rows,_=gauss_jordan(augmented) solution_exists=len(rref_rows)==len(set(map(tuple,reduced_form[:,-1]))) print("Solution Exists:",solution_exists) ``` 以上脚本片段仅作为概念验证用途并未完全覆盖边界情形检测等功能完善需求实际部署前还需进一步测试改进确保鲁棒性和兼容性良好。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值