2019杭电多校 第六场 6641（原1008） TDL（规律题）

2019杭电多校 第六场 6641（原1008） TDL【规律+gcd】

题目：http://acm.hdu.edu.cn/showproblem.php?pid=6641

大致题意：

给你两个正整数k和m，要求你找出符合 (f(n,m)−n)⊕n=k 这个条件式的 n，（⊕表示异或xor，f(n,m)表示求出比n大的第m小的互质数）。

思路：

首先，分析 f(n,m)-n 可以知道，我们求出n的第m小的互质数，然后又减去了n，题目中可以看到m的范围是1~100，所以说m最大不超过100，通过暴力代码可以发现这样一个规律：


用暴力代码随便跑几个数的最近的100个互质数，可以发现几乎都是在1000范围以内（甚至更小的范围）就可以找完这个数的100个互质数，那么我们可以知道 f(n,m)-n 的取值范围大概不超过 1000。

接下来，从二进制的角度分析 (f(n,m)−n)⊕n=k 。

由题可知，1≤k≤1e18 ，假设k很大，取1e15，为了方便，下面我们令

 p = (f(n,m)−n)    即 p ⊕ n = k

因为p的值非常小，k值非常大，然而 p ⊕ n = k，从二进制异或的角度来看，存在 n ≈ k 。所以得出结论 n就在k的附近 。

所以我们用暴力代码，直接搜索 k±1000 范围内的数字是否符合条件，如果存在则输出 n，不存在则输出 -1 。

代码部分：

暴力代码：

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;

//这是一个很神奇的 gcd() 代码
int gcd(int a, int b) {
	while(b^=a^=b^=a%=b);
	return a;
}

int main() {
	int n;
	const int m = 100;
	while(~scanf("%d", &n)) {
		int count = 0;
		int num = 1;
		//count 是用来控制输出100个互质数的
		while(count != 100) {
		//补充一下概念：两个数的最大公约数为1，则为互质数
			if(gcd(n+num, n) == 1) {
				printf("%d ", n+num);
				count++;
			}
			num++;
		}
		cout << endl << endl;
	}
	return 0;
}

ac代码：

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;
#define ll long long
ll k, m;

ll gcd(ll a, ll b) {
	while(b^=a^=b^=a%=b);
	return a;
}

void solve() {
	ll n = k-1000;
	if(n <= 0)
		n = 1;
	
	ll result = -1;
	//p用来暴力寻找n的互质数
	//count表示第几小互质数
	while(n <= k+1000) {
		ll count = 0, p = 0;
		while(count != m) {
			p++;
			if(gcd(n + p, n) == 1)
				count++;
		}
		ll temp = p ^ n;
		if(temp == k) {
			result = n;
			break;
		}
		n++;
	}
	
	if(result != -1)
		printf("%lld\n", result);
	else
		printf("%lld\n", result);
	return;
}

int main() {
	int T;
	while(~scanf("%d",&T)) {
		while(T--) {
			scanf("%lld %lld",&k, &m);
			solve();
		}
	}
	return 0;
}

总结

当时做，几乎要放弃的一道题目，因为完全没有头绪，但就是总感觉(f(n,m)−n)是一个非常神奇非常关键的东西，看到那么多人过自己过不了好不甘心啊，就一直盯着看，幸亏最后做出来了哈哈哈，ac的那一刻的喜悦无法言说，爽！rua！！！冲冲冲。