poj6313(构造)

构造一个矩阵，任意子矩阵中没有四个点都为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;
}