bzoj3034: Heaven Cow与God Bull

Description

__int64 ago,there's a heaven cow called sjy...
A god bull named wzc fell in love with her...
As an OI & MOer,wzc gave sjy a quesiton...

给定一个整数n，求一个整数m，满足m<=n，并且m/phi(m)的值最大。
注：phi(m)代表m的欧拉函数，即不大于m且与m互质的数的个数。

Input

第一行是一个整数T，表示该测试点有T组数据。
接下来T行，每行一个整数n，意义如上所述。

Output

输出一共T行，每行一个整数m。
若对于某个n，有不止一个满足条件的m，则输出最小的m。

贪心考虑，答案为质数的前缀积，筛出60000以内的质数，离线处理询问
高精度需要压位优化

#include<cstdio>
#include<cstring>
int ps[60000],pp=0;
bool isnp[60004];
int T;
char str[30000];
int p10[]={1,10,100,1000};
struct integer{
int x[7000],id,ed;
void read(){
scanf("%s",str);
int l=strlen(str)-1;
for(int i=0;i<=l;i++){
x[i>>2]+=p10[i&3]*(str[l-i]-'0');
}
}
void print(){
int p=6999;
while(p&&!x[p])--p;
printf("%d",x[p]);
for(int i=p-1;~i;--i)printf("%04d",x[i]);
putchar(10);
}
void operator*=(int a){
for(int i=0;i<7000;i++)x[i]*=a;
for(int i=0;i<6999;i++)x[i+1]+=x[i]/10000,x[i]%=10000;
}
}q[100],v,v1;
bool operator>(integer&a,integer&b){
for(int i=6999;~i;--i){
if(a.x[i]!=b.x[i])return a.x[i]>b.x[i];
}
return 0;
}
int main(){
for(int i=2;i<=60000;i++){
if(!isnp[i])ps[pp++]=i;
for(int j=0;j<pp&&i*ps[j]<=60000;j++){
isnp[i*ps[j]]=1;
if(i%ps[j]==0)break;
}
}
scanf("%d",&T);
for(int i=0;i<T;i++){
q[i].id=i;
q[i].read();
}
v.x[0]=v1.x[0]=1;
for(int i=0;i<pp;i++){
v*=ps[i];
for(int j=0;j<T;j++)if(!q[j].ed&&v>q[j]){
q[j]=v1;
q[j].ed=1;
}
v1=v;
}
for(int i=0;i<T;i++){
for(int j=0;j<T;j++)if(q[j].id==i){
q[j].print();
}
}
return 0;
}