构造一个矩阵,任意子矩阵中没有四个点都为1
先看构造的公式:
for (int i = 0; i < p; i++)
for (int j = 0; j < p; j++)
for (int k = 0; k < p; k++)
mat[i*p+j][k*p+(j*k+i)%p] = 1;
再贴上图:
循环的i相当于大行行数,j相当于小行行数,k相当于大列列数,那么k*p+(j*k+i)%p就是计算出来的小列列数
找到某一行中的两个值为1的点,设坐标为(k1,k2)(k3,k4)其中k1,k3是大列坐标,k2,k4小列坐标,
带入方程k2 = (j*k1+i) % p 相减可得到(k1-k3)*j = (k2-k4) mod p
可以唯一解出i和j,也就是说唯一存在这一行,使得这两列上的值是1,如果有存在两行不就构成矩阵了对吧
现在证明为什么当p为质数的时候有唯一解:
因为按我们这样构造,如果要出现子矩阵四个点唯一的情况,一定只能是在同一大行里才有可能,也就是i相同
i相等了,就要证明j唯一:(k1-k3)*j = (k2-k4) mod p这个式子里,j的范围是(0-p),因为p是质数所以下一个解应该是(j + p),超过了(0-p)的范围,所以p一定要为质数才能使解唯一
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int mx = 2505;
int mat[mx][mx];
int main() {
memset(mat, 0, sizeof(mat));
int p = 47;
for (int i = 0; i < p; i++)
for (int j = 0; j < p; j++)
for (int k = 0; k < p; k++)
mat[i*p+j][k*p+(j*k+i)%p] = 1;
puts("2000");
for (int i = 0; i < 2000; i++) {
for (int j = 0; j < 2000; j++)
printf("%d",mat[i][j]);
puts("");
}
return 0;
}