hdu 5750 Dertouzos(数论:质数枚举)

2024-03-28 12:48

本文主要是介绍hdu 5750 Dertouzos(数论:质数枚举),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!


大意:求1到n内有多少个数(以下称为x),x的最大因数是d,也可以理解为x能整除d且能被x整除的最大数为d

思路:分两种情况:

1、当d为质数,直接枚举一个k,使得k * d = x <n,即可,考虑k为质数,如果k为合数,那么k的因子乘d会比d大。比如最后一组样例 k = 4,d = 13,那么4 * 13 = 52,52的最大因子是2*13 = 26而非13,所以枚举素数k到n/d。

    2、当d为合数,那么枚举的k会小于等于d的最小因子。例如k = 5,d = 12,k * d = 60。其中5*3比12大,且为60的因子(60还有更大因子),15是由5 * 3得来的,即12分解为3*2*2,所以这里同样枚举k <=2。


综上所述  k < min(n/d,d的最小因子) 即可 


#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int MAXN = 100001;/*
bool vis[N];
int prime[N>>1],tot;
void getprime(){for(int i=2;i*i<=N-5;++i){if(vis[i])continue;for(int j=i*i;j<=N-5;j+=i)vis[j]=true;}for(int i=2;i<=N-5;++i)if(!vis[i])prime[++tot]=i;
}
*/int prime[MAXN+1];void getprime(){memset(prime,0,sizeof prime);for(int i = 2;i<=MAXN;i++){if(!prime[i]) prime[++prime[0]] = i;for(int j =1;j<= prime[0] && prime[j]<=MAXN/i;j++){prime[prime[j]*i] = 1;if(i % prime[j] == 0) break;}}
}
int main()
{getprime();//cout << prime[2] <<endl;int t;scanf("%d",&t);while(t --){ll n,d;scanf("%I64d%I64d",&n,&d);ll ans = 0;for(int i =1;i<MAXN&&d*prime[i]<n;i++){ans++;if(d % prime[i] == 0) break;}printf("%I64d\n",ans);}return 0;
}


这篇关于hdu 5750 Dertouzos(数论:质数枚举)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/855665

相关文章

C#实现获得某个枚举的所有名称

《C#实现获得某个枚举的所有名称》这篇文章主要为大家详细介绍了C#如何实现获得某个枚举的所有名称,文中的示例代码讲解详细,具有一定的借鉴价值,有需要的小伙伴可以参考一下... C#中获得某个枚举的所有名称using System;using System.Collections.Generic;usi

Java 枚举的常用技巧汇总

《Java枚举的常用技巧汇总》在Java中,枚举类型是一种特殊的数据类型,允许定义一组固定的常量,默认情况下,toString方法返回枚举常量的名称,本文提供了一个完整的代码示例,展示了如何在Jav... 目录一、枚举的基本概念1. 什么是枚举?2. 基本枚举示例3. 枚举的优势二、枚举的高级用法1. 枚举

Rust中的Option枚举快速入门教程

《Rust中的Option枚举快速入门教程》Rust中的Option枚举用于表示可能不存在的值,提供了多种方法来处理这些值,避免了空指针异常,文章介绍了Option的定义、常见方法、使用场景以及注意事... 目录引言Option介绍Option的常见方法Option使用场景场景一:函数返回可能不存在的值场景

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2093 考试排名(sscanf)

模拟题。 直接从教程里拉解析。 因为表格里的数据格式不统一。有时候有"()",有时候又没有。而它也不会给我们提示。 这种情况下,就只能它它们统一看作字符串来处理了。现在就请出我们的主角sscanf()! sscanf 语法: #include int sscanf( const char *buffer, const char *format, ... ); 函数sscanf()和

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

hdu 3790 (单源最短路dijkstra)

题意: 每条边都有长度d 和花费p,给你起点s 终点t,要求输出起点到终点的最短距离及其花费,如果最短距离有多条路线,则输出花费最少的。 解析: 考察对dijkstra的理解。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstrin