任意給定一個正整數N,求一個最小的正整數M(M>1),使得N*M的十進位制表示形式裡只含有1和0。
阿新 • • 發佈:2019-01-24
解決這個問題首先考慮對於任意的N,是否這樣的M一定存在。可以證明,M是一定存在的,而且不唯一。
簡單證明:因為
這是一個無窮數列,但是數列中的每一項取值範圍都在[0, N-1]之間。所以這個無窮數列中間必定存在迴圈節。即假設有s,t均是正整數,且s<t,有 。於是迴圈節長度為t-s。於是10^s = 10^t。因此有:
,所以
例如,取N=3,因為10的任何非負次方模3都為1,所以迴圈節週期為1.有:
給定N,求M的方法:
方法一:給定N,令M從2開始,列舉M的值直到遇到一個M使得N*M的十進位制表示中只有1和0.
方法二:求出10的次方序列模N的餘數序列並找出迴圈節。然後搜尋這個餘數序列,搜尋的目的就是要在這個餘數序列中找到一些數出來讓它們的和是N的倍數。例如N=13,這個序列就是1,10,9,12,3,4然後不斷迴圈。很明顯有1+12=13,而1是10的0次方,12是10的3次方,所以這個數就是1000+1=1001,M就是1001/13=77。
方法三
上圖中括號內表示模3的餘數。括號外表示被搜尋的數。左子樹表示0,右子樹表示1.上圖中搜索到第二層(根是第0層)時遇到111,它模3餘數為0.所以N*M=111, M=111/3=37。
方法四:對方法三的改進。將方法三的搜尋空間按模N餘數分類,使得搜尋時間和空間都由原來的指數級降到了O(N)。改進的原理:假設當前正在搜尋由0,1組成的K位十進位制數,這樣的K位十進位制數共有2^k個。假設其中有兩個數X、Y,它們模N同餘,那麼在搜尋由0、1組成的K+1位十進位制數時,X和Y會被擴展出四個數:10X, 10X+1, 10Y, 10Y+1。因為X和Y同餘(同餘完全可以看作相等),所以10X與10Y同餘,10X+1與10Y+1同餘。也就是說由Y擴展出來的子樹和由X擴充套件產生出來的子樹產生完全相同的餘數,如果X比Y小,那麼Y肯定不是滿足要求的最小的數,所以Y這棵子樹可以被剪掉。這樣,2^K個數按照模N餘數分類,每類中只保留最小的那個數以供擴充套件。原來在這一層需要搜尋2^K個數,現在只需要搜尋O(N)個數。例如,當N=9時,第0層是1(1),
如上圖所示,第2層的110,第三層的1010、1110都因為同一層有和它同餘且更小的數而被剪掉。如果按照方法三搜尋,第三層本來應該有8個結點,但現在只有4個結點。
方法一的原始碼:
- #include <stdio.h>
- int HasOnlyOneAndZero(unsigned int n) {
- while(n) {
- if(n % 10 >= 2) return 0;
- n /= 10;
- }
- return 1;
- }
- int main() {
- int n, m;
-
while(scanf(
- for(m = 1;;m++) {
- if(HasOnlyOneAndZero(n*m)) {
- printf("n = %d, m = %d, n*m = %d/n", n, m, n*m);
- break;
- }
- }
- }
- return 0;
- }
方法三的原始碼:
- // 解法三(1):廣度優先搜尋
- #define _CRT_SECURE_NO_WARNINGS 1
- #include <cstdio>
- #include <queue>
- usingnamespace std;
- int main() {
- int N;
- while(scanf("%d", &N) != EOF) {
- queue<int> q;
- q.push(1);
- while(!q.empty()) {
- int t = q.front();
- q.pop();
- if(t % N == 0) {
- printf("n = %d, m = %d, n*m = %d/n", N, t/N, t);
- break;
- }
- q.push(t * 10);
- q.push(t * 10 + 1);
- }
- }
- return 0;
- }
方法四原始碼:
- // 解法四:將搜尋空間分過類的廣度搜索,這樣空間佔用是O(N)而不是
- // 指數級。分類的原則是按照模N的餘數分類
- #define _CRT_SECURE_NO_WARNINGS 1
- #include <cstdio>
- #include <bitset>
- #include <vector>
- #include <queue>
- usingnamespace std;
- struct QNode {
- int v, r; // v is value, r is remainder
- QNode(int vv, int rr): v(vv), r(rr) {}
- QNode(): v(0), r(0) {}
- };
- int main() {
- int N;
- while(scanf("%d", &N) != EOF) {
- //bitset<N> bn;
- queue<QNode> q;
- q.push(QNode(1, 1));
- while(!q.empty()) {
- //bn.reset();
- vector<bool> bn(N, false);
- int s = q.size();
- while(s--) {
- QNode t = q.front();
- if(t.r == 0) {
- printf("n = %d, m = %d, n*m = %d/n", N, t.v/N, t.v);
- goto ok;
- }
- q.pop();
- if(!bn[t.r * 10 % N]) {
- bn[t.r * 10 % N] = true;
- q.push(QNode(t.v * 10, t.r * 10 % N));
- }
- if(!bn[(t.r * 10 + 1) % N]) {
- bn[(t.r * 10 + 1) % N] = true;
- q.push(QNode(t.v * 10 + 1, (t.r * 10 + 1) % N));
- }
- }
- }
- ok:;
- }
- return 0;
- }