#include <iostream>
using namespace std;
long long p;
int a[100010];
int main()
{
while(cin >> p)
{
if(p==0)break;
if(p==2)cout << "Impossible" << endl;
else
{
for(int i=1;i<p;i++)a[i]=1;
for(long long i=1;i<p;i++)a[i*i%p]=0;//用int,RE一次
for(int i=1;i<p;i++)cout << a[i];
cout << endl;
}
}
return 0;
}
思路:必须承认这道题目是看了网上的解题报告,才知道怎么做的。
a[1%n], a[2%n], a[3%n], ..., a[n-1] (1)
a[2%n], a[4%n], a[6%n], ..., a[2(n-1)%n] (2)
若a[1%n]=a[2%n],则原字符串或者补字符串(complement)第1,2位都是相等的,所以(2)式中a[2%n]=a[4%n](题目中说每一行不是原字符串就是补字符串:If each row of the matrix is either the original bitstring or its complement, the original bitstring is magic. 这句话很关键)
同理若a[1%n]!=a[2%n],则a[2%n]!=a[4%n]
由上可知a[1%n]=a[4%n],同理可推出a[1%n]=a[4%n]=a[9%n]=...=a[i*i%n] (*)
令(*)=0,其他位置为1,即可得到答案。
因为(*)包含a[1]所以要想字典序最小则必须为零
其他位置必须为1的理由如下:
若a[1%n]=a[2%n],则a[2%n]=a[4%n],则a[4%n]=a[8%n],a[3%n]=a[6%n],a[4%n]=a[6%n]...
最终会导致整个矩阵都为零,与题目要求non-constant magic bitstring不符。