HDU 5750 Dertouzos (数论)

Dertouzos

 

 Accepts: 76

 

 Submissions: 1357

 Time Limit: 7000/3500 MS (Java/Others)

 

 Memory Limit: 131072/131072 K (Java/Others)

问题描述

正整数$x$称为$n$的positive proper divisor, 当且仅当$x|n$并且$1\le x<n$. 例如, 1, 2, 和3是6的positive proper divisor, 但是6不是.

Peter给你两个正整数$n$和$d$. 他想要知道有多少小于$n$的整数, 满足他们的最大positive proper divisor恰好是$d$.

输入描述

输入包含多组数据, 第一行包含一个整数$T$ (1 \le T \le 10^6)(1≤T≤10​6​​)表示测试数据组数. 对于每组数据:

第一行包含两个整数$n$和$d$ (2 \le n, d \le 10^9)(2≤n,d≤10​9​​).

输出描述

对于每组数据, 输出一个整数.

输入样例

9
10 2
10 3
10 4
10 5
10 6
10 7
10 8
10 9
100 13

输出样例

1
2
1
0
0
0
0
0
4

随便推导下, 令$y=xd$, 如果$d$是$y$的maximum positive proper divisor, 显然要求$x$是$y$的最小质因子. 令$mp(n)$表示$n$的最小质因子, 那么就有$x\le mp(d)$, 同时有$y<n$, 那么x \le \lfloor \frac{n-1}{d} \rfloorx≤⌊​d​​n−1​​⌋. 于是就是计算有多少个素数$x$满足x \le \min{mp(d), \lfloor \frac{n-1}{d} \rfloor}x≤min{mp(d),⌊​d​​n−1​​⌋}.
当$d$比较大的时候, \lfloor \frac{n-1}{d} \rfloor⌊​d​​n−1​​⌋比较小, 暴力枚举$x$即可. 当$d$比较小的时候, 可以直接预处理出答案. 阈值设置到10^6 \sim 10^710​6​​∼10​7​​都可以过.
当然判断最小质因子很简单，详情请看代码
【个人数论技巧，对于有着最大或者最小值时，一般是通过将公式推导为一个或者多个不等式，让他们同时满足条件，从而达到最终的目的，其中一定要有效的建立起所有能够运用关系，要懂得建立相关方程进行公式化简】
#include <map>
#include <cmath>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>

#define uint unsigned int
using namespace std;
const int MAXN = 5e6 + 5;
int primes[MAXN];
bool vis[MAXN];
int cnt,T, n,d;
void init() {
	cnt = 0;
	memset(vis, false, sizeof(vis));
	for(int i = 2; i < MAXN; i ++) {
		if(vis[i]) continue;
		vis[i] = true;
		primes[cnt ++] = i;
		for(int j = i * 2; j < MAXN; j += i) {
			vis[j] = true;
		}
	}
}
int main() {
	init();
	while(~scanf("%d", &T)) {
		while(T --) {
			scanf("%d%d", &n, &d);
			int l = 0;
			while(primes[l] <= (n - 1) / d) {
				l ++;
				if(d % primes[l - 1] == 0) break;
			}
			printf("%d\n", l);
		}
	}
	return 0;
}