HDU 5317(打表)
2015-10-05 23:29
176 查看
RGCDQ
Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 2261 Accepted Submission(s): 923
Problem Description
Mr. Hdu is interested in Greatest Common Divisor (GCD). He wants to find more and more interesting things about GCD. Today He comes up with Range Greatest Common Divisor Query (RGCDQ). What’s RGCDQ? Please let me explain it to you gradually. For a positive
integer x, F(x) indicates the number of kind of prime factor of x. For example F(2)=1. F(10)=2, because 10=2*5. F(12)=2, because 12=2*2*3, there are two kinds of prime factor. For each query, we will get an interval [L, R], Hdu wants to know maxGCD(F(i),F(j)) (L≤i<j≤R)
Input
There are multiple queries. In the first line of the input file there is an integer T indicates the number of queries.
In the next T lines, each line contains L, R which is mentioned above.
All input items are integers.
1<= T <= 1000000
2<=L < R<=1000000
Output
For each query,output the answer in a single line.
See the sample for more details.
Sample Input
2 2 3 3 5
Sample Output
1 1
//当时多校没做出来。。今天回过头来想想其实并不是很难
//题意:定义一个函数f(x)=x的素数因子的个数 给你L R 求 res=maxGCD(f(i),f(j)) L<=i<j<=R
//数据量大 只能打表
//首先一个素数表是肯定的 接着要求素数因子的个数 我们可以在打表的过程中遇到是素数 他的倍数都加1
//求GCD的时候我们自己算下或者测下数据可以发现最大的素数因子个数是7 可以开一个数组记录前缀和 接着特判下
#include <stdio.h> #include <string.h> #include <algorithm> using namespace std; int vis[1000010]; int prime[1000010]; int res[1000010]; int ans[1000010][9]; void Prime() { memset(vis,0,sizeof(vis)); memset(res,0,sizeof(res)); memset(prime,0,sizeof(prime)); for(int i=2;i<=1000000;i++) //素数打表 { if(!vis[i]) { prime[i]=1; for(int k=i;k<=1000000;k+=i) { vis[k]=1; if(prime[i]) 当i是素数时他的倍数都有一个素数的因子 res[k]++; } } } } void function() { memset(ans,0,sizeof(ans)); for(int i=2;i<=1000000;i++) //求前缀和 { for(int k=1;k<8;k++) ans[i][k]=ans[i-1][k]; ans[i][res[i]]++; } } int main() { Prime(); function(); int t,n,m; scanf("%d",&t); while(t--) { scanf("%d%d",&n,&m); int flag=0; for(int i=7;i>=1;i--) //特判 7 7 6 6 5 5 4 4 3 3 6 3 6 4 6 2 4 2 2 2 1 { if(i>=3) { if(ans[m][i]-ans[n-1][i]>=2) { flag=1; printf("%d\n",i); break; } } else { if(ans[m][6]-ans[n-1][6]==1&&ans[m][3]-ans[n-1][3]==1) { flag=1; printf("3\n"); break; } if((ans[m][6]-ans[n-1][6]==1&&ans[m][4]-ans[n-1][4]==1)||(ans[m][2]-ans[n-1][2]>=2)||(ans[m][4]-ans[n-1][4]==1&&ans[m][2]-ans[n-1][2]==1)||(ans[m][6]-ans[n-1][6]==1&&ans[m][2]-ans[n-1][2]==1)) //刚刚才发现我AC的代码少了个6 2 hdu 数据有点水啊 { flag=1; printf("2\n"); break; } } } if(!flag) printf("1\n"); } return 0; }
相关文章推荐
- Java 集合系列02之 Collection架构
- KVC 和KVO浅谈
- 马哥运维笔记----5
- hdoj 5457 Hold Your Hand 【字典树 + 最小割】
- git一般开发使用流程
- sleepSort和硬币概率问题
- merge vs rebase
- Java 集合系列01之 总体框架
- hdu 3577 Fast Arrangement(区间更新)
- Android面试准备:自定义控件(二)
- Java 集合系列03之 ArrayList详细介绍(源码解析)和使用示例
- 卓越管理的秘密(Behind Closed Doors)
- 树的子结构
- C/C++编程必备网址
- Incompatible pointer types returning 'NSArray *' from a function with result type 'xxx *'
- 好的文章网址
- 爬虫入门——02
- 如何在面试中面试官的提示下,积极地思考?
- ASP精华[转]
- 【动态域名解析工具】tunnel,国内版的ngrok,花生壳可以睡觉了