1104 天长地久 (20 分)

“天长地久数”是指一个 K 位正整数 A，其满足条件为：A 的各位数字之和为 m，A+1 的各位数字之和为 n，且 m 与 n 的最大公约数是一个大于 2 的素数。本题就请你找出这些天长地久数。

输入格式：

输入在第一行给出正整数 N（≤5），随后 N 行，每行给出一对 K（3<K<10）和 m（1<m<90），其含义如题面所述。

输出格式：

对每一对输入的 K 和 m，首先在一行中输出 Case X，其中 X 是输出的编号（从 1 开始）；然后一行输出对应的 n 和 A，数字间以空格分隔。如果解不唯一，则每组解占一行，按 n 的递增序输出；若仍不唯一，则按 A 的递增序输出。若解不存在，则在一行中输出 No Solution。

输入样例：

2
6 45
7 80

输出样例：

Case 1
10 189999
10 279999
10 369999
10 459999
10 549999
10 639999
10 729999
10 819999
10 909999
Case 2
No Solution

分析：

本题主要是数学分析部分比较 ♂，但陈越姥姥说暴搜就能解决，我看使用剪枝的小伙伴们也很多
1、若A为xxxx1~xxxx8，则A+1为 xxxx2~xxxx9，此时gcd(A,A+1)=1，不符合条件
2、若A为xxxx9，eg.gcd(89,90)={1,2,4,8} (因为8+9-9=8)，不符合条件
3、若A为xxx99，eg.gcd(99,100)={1,17} (因为9+9-1=17)，符合条件
4、注意set<pair<int,int> >的使用
（set会自动升序排列，然而这里每一个元素都是pair，它排列是先排pair里的first，再排second。)

代码：

#include<bits/stdc++.h>
using namespace std;
bool isPrime(int n) {	// 判断素数
	for(int i=2; i*i<=n; i++)
		if(n%i==0) return false;	// 非素数
	return true;	// 素数
}
int gcd(int a,int b) {	// 求最大公约数
	return !b?a:gcd(b,a%b);
}
int digit_sum(int n) { // 计算这个数的各位数字之和
	int sum=0;
	while(n) {
		sum+=n%10;
		n/=10;
	}
	return sum;
}
set<pair<int,int> > res;	// 结果集
int main() {
	int n,k,m,cnt;	// k位，各位数字之和为 m
	cin>>n;
	for(int i=0; i<n; i++) {
		cin>>k>>m;
		printf("Case %d\n",i+1);
		cnt=0;
		for(int j=99+pow(10,k-1); j<=pow(10,k)-1; j+=100) {
			if((digit_sum(j)==m)
					&&(gcd(m,digit_sum(j+1))>2)
					&&isPrime(gcd(m,digit_sum(j+1)))) { // 符合条件 
				res.insert(make_pair(digit_sum(j+1),j));
				cnt++;
			}
		}
		if(cnt==0) cout<<"No Solution"<<endl;
		else{
			for(auto it=res.begin();it!=res.end();it++){
				cout<<it->first<<" "<<it->second<<endl;
			}
		}
		res.clear();
	}
	return 0;
}

提交记录：

数学题，思路非常的精巧 ~