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≤x<n. 例如, 1, 2, 和3是6的positive proper divisor, 但是6不是. Peter给你两个正整数n和d. 他想要知道有多少小于n的整数, 满足他们的最大positive proper divisor恰好是d.
输入描述
输入包含多组数据, 第一行包含一个整数T (1≤T≤106)表示测试数据组数. 对于每组数据: 第一行包含两个整数n和d (2≤n,d≤109).
输出描述
对于每组数据, 输出一个整数.
输入样例
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≤mp(d), 同时有y<n, 那么x≤⌊dn−1⌋. 于是就是计算有多少个素数x满足x≤min{mp(d),⌊dn−1⌋}.
当d比较大的时候, ⌊dn−1⌋比较小, 暴力枚举x即可. 当d比较小的时候, 可以直接预处理出答案. 阈值设置到106∼107都可以过.
当然判断最小质因子很简单,详情请看代码
【个人数论技巧,对于有着最大或者最小值时,一般是通过将公式推导为一个或者多个不等式,让他们同时满足条件,从而达到最终的目的,其中一定要有效的建立起所有能够运用关系,要懂得建立相关方程进行公式化简】
#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; }