POJ 1426 bfs / dp

大致题意

给定一个正整数 n，找到一个非零的 n 的倍数 m，使其十进制表示每一位都有 0, 1 组成，且不超过 100 位。1 <= n <= 200

dfs

bfs 枚举每一位会超时，n 值较小，可以打表。发现在 long long 范围内即可搜索出答案，也可以使用限制搜索深度的 dfs。

#include <cstdio>
#include <stdlib.h>
#include <algorithm>
#include <iostream>
#define min(a,b)    (((a) < (b)) ? (a) : (b))
#define max(a,b)    (((a) > (b)) ? (a) : (b))
#define abs(x)    ((x) < 0 ? -(x) : (x))
#define INF 0x3f3f3f3f
#define delta 0.85
#define eps 1e-3
#define M_PI 3.14159265358979323846
#define MAX_N 10
using namespace std;
typedef long long ll;

bool dfs(ll m, int dep, int n){
	if(m % n == 0){
		cout << n << endl;
		return 1;
	}
	if(dep > 18) return 0;
	return (dfs(m * 10, dep + 1, n) || dfs(m * 10 + 1, dep + 1, n)) ? 1 : 0;
}

int main(){
	int n;
	while(cin >> n && n){
		dfs(1, 0, n);
	}
	return 0;
}

dp

dfs 耗时仍然不理想，于是观摩了 dp 写法。

手算除法中（解 res 除以 n），被除数的高位对 n 的余数乘以 10 加上低一位的数，即为低一位的被除数。dp 下标的二进制表示代表对应的十进制，比如下标 (110)₂ 代表可能解 (110)₁₀ ，从高位到低位树状动态规划。

这方法局限与 dfs 一致，在不知道可能解的最高位范围时，仅根据题中 m 不超过 100 位的条件，估复杂度会爆。

#include <cstdio>
#include <stdlib.h>
#include <algorithm>
#include <iostream>
#include <string>
#define min(a,b)    (((a) < (b)) ? (a) : (b))
#define max(a,b)    (((a) > (b)) ? (a) : (b))
#define abs(x)    ((x) < 0 ? -(x) : (x))
#define INF 0x3f3f3f3f
#define delta 0.85
#define eps 1e-3
#define M_PI 3.14159265358979323846
#define MAX_N 10
using namespace std;
int dp[1 << 19]; //最高 18 位

int solve(int n){
	dp[1] = 1 % n;
	for(int i = 1;; i++){
		// 十进制该位为 0 或 1 的情况，余数为 0 表示可被 n 整除
		if(!(dp[i << 1] = dp[i] * 10 % n)) return i << 1;
		if(!(dp[(i << 1) + 1] = (dp[i] * 10 + 1) % n)) return (i << 1) + 1;
	}
}

int main(){
	int n;
	while(cin >> n && n){
		int loc = solve(n);
		string res = "";
		while(loc > 0){
			if(loc % 2 == 1) res = "1" + res;
			else res = "0" + res;
			loc /= 2;
		}
		cout << res << endl;
	}
	return 0;
}

为了保证复杂度可以求解，可以如下考虑：由 0, 1 组成的非零十进制数字，看成由 10ⁱ (0 <= i <= 100) 组合而成。考虑到求余的性质以及 m % n == 0 ，则

$\sum 10^i=0mod(n)$

转化为多重部分和问题，按位 dp 求解即可。

#include <cstdio>
#include <stdlib.h>
#include <iostream>
#include <string>
#define min(a,b)    (((a) < (b)) ? (a) : (b))
#define max(a,b)    (((a) > (b)) ? (a) : (b))
#define abs(x)    ((x) < 0 ? -(x) : (x))
#define INF 0x3f3f3f3f
#define delta 0.85
#define eps 1e-3
#define M_PI 3.14159265358979323846
#define MAX_N 10
using namespace std;

void solve(int n){
	if(n == 1){
		cout << 1 << endl;
		return;
	}
	string dp[101][200];
	int mod = 1 % n;
	dp[1][0] = "0", dp[1][mod] = "1";
	for(int i = 1; i < 100; i++){
		mod = (mod * 10) % n;
		for(int j = 0; j < n; j++){
			if(dp[i][j] != ""){
				//不选择 10 ^ i
				dp[i + 1][j] = "0" + dp[i][j];
				//选择 10 ^ i
				int nxt = (mod + j) % n;
				dp[i + 1][nxt] = "1" + dp[i][j];
				//出现模为 0 且第一位为 1 的数即满足条件
				if(nxt == 0){
					cout << dp[i + 1][nxt] <<endl;
					return;
				}
			}
		}
	}
}

int main(){
	int n;
	while(cin >> n && n){
		solve(n);
	}
	return 0;
}