A 和 A + 1 肯定是互质的,因此A必然以 9 结尾。然后找规律:
末尾一个9, m 与 n 相差 8;
末尾两个9, m 与 n 相差 17;
可得 m - n + 1 为 9 的倍数,从而可以较快遍历找到 n, 找到 n 后,根据倍数可以知道末尾 9 的个数,对 A 头部非 9 的部分进行遍历,快速得到 A;
#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;
struct node{
int A, m, n;
};
bool cmp(node a, node b){
return a.n != b.n? a.n < b.n : a.A < b.A;
}
int gcd(int a, int b){
return b == 0? a : gcd(b, a % b);
}
bool isprime(int n){
if(n <= 1) return false;
for(int i = 2; i <= sqrt(n); i++){
if(n % i == 0) return false;
}
return true;
}
int getsum(int n){
int ans = 0;
while(n){
ans += n % 10;
n /= 10;
}
return ans;
}
int main() {
int N, k, m;
cin >> N;
for(int i = 1; i <= N; i++){
cin >> k >> m;
vector<node> ans;
printf("Case %d\n", i);
for(int j = 2; j <= m + 1; j++){
if(isprime(gcd(j, m)) && gcd(j, m) > 2 && (m - j + 1) % 9 == 0){
int a = (m - j + 1) / 9;
int A = 0;
for(int i = 0; i < a; i++) A = A * 10 + 9;
for(int u = pow(10, k - a - 1); u < pow(10, k - a); u++){
if(getsum(u) + a * 9 == m){
int temp = pow(10, a) * u + A;
ans.push_back({temp, m, j});
}
}
}
}
if(ans.size() == 0) printf("No Solution\n");
else{
sort(ans.begin(), ans.end(), cmp);
for(int j = 0; j < ans.size(); j++){
printf("%d %d\n", ans[j].n, ans[j].A);
}
}
}
return 0;
}