【HDOJ】1153 Magic Bitstrings【组合数学】

传送门:
Magic Bitstrings

思路

菜逼终究是菜逼,首先表示题目没读懂。看了别人的翻译之后,总算读懂题了。
然后公式推不出来,菜哭了。

把矩阵列出来
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[3%n], a[6%n], a[9%n], ..., a[3(n-1)%n]       (3)
...

比较 (1), (2)
发现:
如果 a[1%n] != a[2%n],那么 a[2%n] != a[4%n],那么 a[1%n] == a[4%n];
如果 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[1%n] == a[9%n],
a[1%n] == a[16%n],
....
所有下标是 i 平方 mod n 都相等。
下标不是 i 平方 mod n 都相等。

为了字典顺序最小,并且避免整个数列都是同一个数(题目要求 non-constant),令 a[1%n] = a[4%n] = a[9%n] = ... = 0,其他数都是 1,这样符合题意。

My AC Code Relying Network

参考 https://blog.csdn.net/chengouxuan/article/details/6877054

/*
核心:推导+找规律
solution:令 a[1%n] = a[4%n] = a[9%n] = ... = 0,其他数都是 1
*/
#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=1e5+5;

int A[maxn];

int main()
{
	int n;
	while(scanf("%d",&n))
	{
		if(n==0)
			break;
		if(n==2)
		{
			printf("Impossible
");
			continue;
		}
		for(int i=1;i<=n-1;i++)
		{
			A[i]=1;
		}
		//这里要用long long 10^10 会溢出
		for(long long i=1;i<=n-1;i++)
		{
			A[(i*i)%n]=0;
		}
		for(int i=1;i<=n-1;i++)
			printf("%d",A[i]);
		printf("
");
	}		
    return 0;
}
原文地址:https://www.cnblogs.com/shengwang/p/9768534.html