|
@@ -1772,7 +1772,6 @@ Vt*vxor(St*st,vt*self,Vt*x,Vt*y){if(xt==NUM&&yt==NUM)R Vnn(((int)xn)^((int)yn));
|
|
Vt*vband(St*st,vt*self,Vt*x,Vt*y){if(xt==NUM&&yt==NUM)R Vnn(((int)xn)&((int)yn));R _NAN;}
|
|
Vt*vband(St*st,vt*self,Vt*x,Vt*y){if(xt==NUM&&yt==NUM)R Vnn(((int)xn)&((int)yn));R _NAN;}
|
|
Lst*find_primes(U64 limit){B sieve[limit+1];
|
|
Lst*find_primes(U64 limit){B sieve[limit+1];
|
|
for(U64 i=0;i<=limit;i++)sieve[i]=F;
|
|
for(U64 i=0;i<=limit;i++)sieve[i]=F;
|
|
-
|
|
|
|
if(limit>2)sieve[2]=T;
|
|
if(limit>2)sieve[2]=T;
|
|
if(limit>3)sieve[3]=T;
|
|
if(limit>3)sieve[3]=T;
|
|
for(U64 x=1;x*x<=limit;x++)for(U64 y=1;y*y<=limit;y++){U64 n=(4*x*x)+(y*y);
|
|
for(U64 x=1;x*x<=limit;x++)for(U64 y=1;y*y<=limit;y++){U64 n=(4*x*x)+(y*y);
|
|
@@ -1784,7 +1783,6 @@ Lst*find_primes(U64 limit){B sieve[limit+1];
|
|
for(U64 r=5;r*r<=limit;r++)if(sieve[r])for(I i=r*r;i<=limit;i+=r*r)sieve[i]=F;
|
|
for(U64 r=5;r*r<=limit;r++)if(sieve[r])for(I i=r*r;i<=limit;i+=r*r)sieve[i]=F;
|
|
Lst*r=ln();
|
|
Lst*r=ln();
|
|
for(U64 a=1;a<=limit;a++)if(sieve[a])lp(r,Vnn(a));R r;}
|
|
for(U64 a=1;a<=limit;a++)if(sieve[a])lp(r,Vnn(a));R r;}
|
|
-
|
|
|
|
Vt*vprimes(St*st,vt*self,Vt*x){if(xt==NUM&&!spnp(xn))R Vna(find_primes(fabs(xn)+1));R st->udf;}
|
|
Vt*vprimes(St*st,vt*self,Vt*x){if(xt==NUM&&!spnp(xn))R Vna(find_primes(fabs(xn)+1));R st->udf;}
|
|
Vt*vparts(St*st,vt*self,Vt*x,Vt*y){if(xt!=NUM)R st->udf;
|
|
Vt*vparts(St*st,vt*self,Vt*x,Vt*y){if(xt!=NUM)R st->udf;
|
|
if(yt!=ARRAY)y=venlist(st,N,y);
|
|
if(yt!=ARRAY)y=venlist(st,N,y);
|