HDU 6313 Hack It（构造+数论）

Description

给出下面一个代码用来判断一个$n\times n$的$01$矩阵中是否存在一个子矩形的四个角的值均为$1$

let count be a 2d array filled with 0s
iterate through all 1s in the matrix:
  suppose this 1 lies in grid(x,y)
  iterate every row r:
    if grid(r,y)=1:
      ++count[min(r,x)][max(r,x)]
      if count[min(r,x)][max(r,x)]>1:
        claim there is a rectangle satisfying the condition
claim there isn't any rectangle satisfying the condition

要求构造一个$n\times n$的$01$矩形使得$1\le n\le 2000$且至少有$85000$个位置是$1$，使得其不存在四个角均为$1$的子矩形但是上述代码的运行结果会判断为存在

Input

无

Output

输出$n$以及这个$n\times n$的$01$矩阵

Sample Input

无

Sample Output

无

Solution

显然只要一行有两个$1$且不存在四个角都为$1$的矩形即为反例，问题转化为构造$n$个序列使得任意两个序列不存在两个都为$1$的位置

给定$0\le a,b,c,d<p$，其中$p$为素数，假设对于两个不同的$k$满足$a+kb=c+kd\ mod\ p$，则有$a=c,b=d$，故取不同的二元组$(a,b)$，通过$a+kb\ mod\ p$可以构造出来$p^2$个长度为$p$的序列，且这些序列中任意两个序列不存在两个位置的值相同，那么我们取素数$47$满足$47^2>2000$，把$2000$分成若干长度为$47$的段，第$k$段的第$a+kb\ mod\ p$个值为$1$，其余值为$0$，则可以构造出满足条件的矩阵

Code

#include<cstdio>
using namespace std;
int main()
{
    int n=2000,k=47;
    printf("%d\n",n);
    int num=0;
    for(int i=0;i<k;i++)
        for(int j=0;j<k;j++)
        {
            int m=i*k+j+1;
            if(m<=n)
            {
                int res=0;
                for(int x=0;x<43;x++)
                {
                    int t=(i+x*j)%k+1;
                    for(int l=1;l<=k;l++)
                    {
                        res++;
                        if(l!=t)printf("0");
                        else printf("1");
                        if(res==n)
                        {
                            printf("\n");
                            break;
                        }
                    }
                }
            }
            else break;
        } 
    return 0;
}