【P1203】买花

时间:2024-04-18 20:04:07

我先在已经弱到连高精乘单精都能写错的地步了QAQ

原题:

求一个小于等于N的数M,使得phi(M)/M最小,
其中phi(M)是与M互质且比M小的数的个数。
例如phi(4)=2,因为1,3和4互质。

N<=10^40

先推欧拉函数打表,发现答案一定是乘积<=m的前k个质数的乘积(比如50的答案就是2*3*5<=50)

然后就是高精水体

然而我写了2h QAQ,果然实力是会随着时间的推移降低的QAQ

主要错误就是关于乘出来结果的长度的问题,这个长度要在a[i]>=10往后推得时候变推边更新,不能先全部算完了再更新,因为推完后可能会因为%=10而出现中间有0的情况,这时候推长度推到这里就断了,但是边推边更新就不会出现这种情况,因为还没%=10,不会==0

我好弱啊QAQ

代码:

 #include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
void read_big(int *x,int &lx){lx=; char ch=getchar();
while(ch<''||ch>'') ch=getchar();
while(ch>=''&&ch<=''){ x[++lx]=ch-''; ch=getchar();}
}
int n[],ln;
int a[],la,b[],lb,c[];
bool kang[]; int zhi[],ztop=;
void shai(){
memset(kang,,sizeof(kang));
for(int i=;i<=;i++)if(!kang[i]){
zhi[++ztop]=i;
int temp=; while(i*(temp+)<=) kang[i*(++temp)]=true;
}
}
bool bi(int *x,int lx,int *y,int ly){
if(lx!=ly) return (lx>ly);
for(int i=lx;i>=;i--)if(x[i]!=y[i]) return (x[i]>y[i]);
return false;
}
void cheng(int *x,int &lx,int y){
for(int i=;i<=lx;i++) x[i]*=y;
for(int i=;i<=lx;i++) x[i+]+=x[i]/,x[i]%=;
for(;a[lx+];lx++) x[lx]+=x[lx]/,x[lx]%=;
}
int main(){//freopen("ddd.in","r",stdin);
shai();
read_big(c,ln);
for(int i=ln;i>=;i--) n[i]=c[ln-i+];
a[la=]=;
for(int i=;i<=;i++){
lb=la;
for(int j=;j<=lb;j++) b[j]=a[j];
cheng(a,la,zhi[i]);
if(bi(a,la,n,ln)){ for(int j=lb;j>=;j--)cout<<b[j]; cout<<endl; return ;}
}
return ;
}