题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5317
题意:给定一个区间[l,r],求出区间内每个数的质因数的种类,并求种类数gcd的最大值
思路:刚开始做的时候考虑的是打表求出前1000000的每个数的质因数种类数记录在一个数组里,打表的时候发现前1000000个数中,种类最多是7种,那么我们可以把种类数相同的数字,记录在同一个数组中,或者数组的同一层,所以求gcd的最大值的时候,就可以直接根据种类数的多少直接判断,比如如果种类数是7的数字大于1个,那么gcd最大值就是7,否则如果种类数是6的大于1,那么gcd最大值就是6,以此类推,然而求质因数种类的算法不够优越,直接t了,借鉴一下别人求质因数的方法,
代码:
#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
#define maxn 1000005
using namespace std;
int f[maxn];
bool flag[maxn];
int s[10][maxn];
void init()
{
memset(f,0,sizeof(f));
memset(flag,0,sizeof(flag));
memset(s,0,sizeof(s));
for(int i=2;i<=1000000;i++)
if(!flag[i])
{
f[i]++;
for(int j=i+i;j<=1000000;j+=i)
{
flag[j]=1;
f[j]++;
}
}//找到每个数字对应的质因数种类数
for(int i=2;i<=1000000;i++)
for(int j=1;j<=7;j++)
{
if(f[i]==j)
s[j][i]=s[j][i-1]+1;
else
s[j][i]=s[j][i-1];
}
}
//记录前i个种类中,j有几个、、
int main()
{
init();
int t,l,r;
int cnt[8];
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&l,&r);
memset(cnt,0,sizeof(cnt));
for(int i=1;i<=7;i++)
{
cnt[i]=s[i][r]-s[i][l-1];
}
if(cnt[7]>=2)
printf("7\n");
else if(cnt[6]>=2)
printf("6\n");
else if(cnt[5]>=2)
printf("5\n");
else if(cnt[4]>=2)
printf("4\n");
else if(cnt[3]+cnt[6]>=2)
printf("3\n");
else if(cnt[2]+cnt[4]+cnt[6]>=2)
printf("2\n");
else
printf("1\n");
}
return 0;
}