User:Anatoly E. Voevudko/MyOEISpariScripts

My sequences and my PARI Scripts

A245713

  A245713(lim)={my(L=List(1),lim2=logint(lim,2));for(p=3,lim2,
  for(b=2,sqrtnint(lim,p),listput(L, b^p);));listsort(L); print(L)}
  \\ _Anatoly E. Voevudko_, Sep 21 2015

  b245713(lim)={my(L=List(1),lim2=logint(lim,2),ln,i);for(p=3,lim2,
  for(b=2,sqrtnint(lim,p),listput(L, b^p))); listsort(L);ln=#L;
  for(i=1,ln,print(i ," ",L[i]))}  \\ _Anatoly E. Voevudko_, Sep 29 2015
  [deleted by me]

  A245713(lim,bflag=0)={my(L=List(1),lim2=logint(lim,2),ln,i);for(p=3,lim2,
  for(b=2,sqrtnint(lim,p),listput(L, b^p))); listsort(L);ln=#L;
  if(bflag, for(i=1,ln,print(i ," ",L[i]))); if(!bflag, return(Vec(L))); }
  \\ _Anatoly E. Voevudko_, Nov 23 2015 [not published]
  

A261782

  list(lim)=my(v=List(), u=v, t);
  for(z=3, logint(lim\=1, 2), for(C=2, sqrtnint(lim, z), listput(v, C^z)));
  v=Set(v); for(i=1, #v, for(j=i, #v, t=v[i]+v[j]; if(t>lim, break);
  if(setsearch(v, t), listput(u, t)))); Set(u) \\ _Charles R Greathouse IV_, Sep 03 2015
 
  b261782(lim)={my(Lcz=List(1),Lb=List(),czn,lczn,lbn,lim2=logint(lim,2),lim3);
  for(z=3,lim2,lim3=sqrtnint(lim,z); for(C=2,lim3,listput(Lcz,C^z)));
  Lcz=Set(Lcz); lczn=#Lcz; if(lczn==0,return(-1));
  for(i=1,lczn, for(j=i,lczn,czn=Lcz[i]+Lcz[j]; if(czn>lim, break);
  if(setsearch(Lcz,czn), listput(Lb,czn)))); listsort(Lb,1); lbn=#Lb;
  for(i=1,lbn, print(i," ",Lb[i]));}
  \\ _Anatoly E. Voevudko_, Sep 23 2015 [not published, based on CRG4 script]
  
  A261782(lim,bflag=0)={my(Lcz=List(1),Lb=List(),czn,lczn,lbn,lim2=logint(lim, 2),lim3);
  for(z=3,lim2, lim3=sqrtnint(lim, z); for(C=2,lim3, listput(Lcz, C^z)));
  Lcz=Set(Lcz); lczn=#Lcz; if(lczn==0, return(-1));
  for(i=1,lczn, for(j=i,lczn, czn=Lcz[i]+Lcz[j]; if(czn>lim, break);
  if(setsearch(Lcz,czn), listput(Lb,czn)))); listsort(Lb,1); lbn=#Lb;
  if(bflag, for(i=1,lbn, print(i," ",Lb[i]))); if(!bflag, return(Vec(Lb))); }
  \\ _Anatoly E. Voevudko_, Nov 23 2015 [not published, based on CRG4 script]

A264901

  b264901(lim)=
  {my(Lc=List(1),Lb=List(),La=Lb,czn,lan,lbn,lcn,lim2=logint(lim,2),lim3);
  for(z=3,lim2, lim3=sqrtnint(lim, z); for(C=2,lim3, listput(Lc, C^z)));
  lcn=#Lc; if(lcn==0, return(-1));
  for(i=1,lcn, for(j=i,lcn, czn=Lc[i]+Lc[j]; if(czn>lim, next);
  La=findinlista(Lc, czn); lan=#La; if(!lan, next);
  for(k=1,lan, listput(Lb,czn)); )); lbn=#Lb; listsort(Lb);
  for(i=1,lbn, print(i," ",Lb[i]))}
  findinlista(list,item,sind=1)=
  {my(ln=#list,Li=List()); if(ln==0 || sind<1 || sind>ln, return(Li));
  for(i=sind,ln, if(list[i]==item, listput(Li,i))); return(Li); }
  \\ _Anatoly E. Voevudko_, Nov 21 2015
  
  A264901(lim,bflag=0)=
  {my(Lc=List(1),Lb=List(),La=Lb,czn,lan,lbn,lcn,lim2=logint(lim,2),lim3);
  for(z=3,lim2, lim3=sqrtnint(lim, z); for(C=2,lim3, listput(Lc,C^z)));
  lcn=#Lc; if(lcn==0, return(-1));
  for(i=1,lcn, for(j=i,lcn, czn=Lc[i]+Lc[j]; if(czn>lim, next);
  La=findinlista(Lc, czn); lan=#La; if(!lan, next);
  for(k=1,lan, listput(Lb,czn)); )); lbn=#Lb; listsort(Lb);
  if(bflag, for(i=1,lbn, print(i ," ",Lb[i]))); if(!bflag, return(Vec(Lb))); }
  findinlista(list, item, sind=1)=
  {my(ln=#list,Li=List()); if(ln==0 || sind<1 || sind>ln, return(Li));
  for(i=sind,ln, if(list[i]==item, listput(Li,i))); return(Li); }
  \\ _Anatoly E. Voevudko_, Nov 23 2015 [not published]
    

A265731

  A265731(lim,bflag=0)=
  {my(Lcz=List(1),Lb=List(),czn,lczn,lbn,lim2=logint(lim,2),lim3);
  for(z=2,lim2, lim3=sqrtnint(lim,z); for(C=2,lim3, listput(Lcz, C^z)));
  Lcz=Set(Lcz); lczn=#Lcz; if(lczn==0, return(-1));
  for(i=1,lczn, for(j=i,lczn, czn=Lcz[i]+Lcz[j]; if(czn>lim, break);
  if(setsearch(Lcz,czn), listput(Lb,czn)))); listsort(Lb,1); lbn=#Lb;
  if(bflag, for(i=1,lbn, print(i," ",Lb[i]))); if(!bflag, return(Vec(Lb))); }
  \\ _Anatoly E. Voevudko_, Nov 23 2015
  

A265732

  A265732(lim,bflag=0)=
  {my(Lc=List(1),Lb=List(),La=Lb,czn,lcn,lan,lim2=logint(lim,2), lim3,k);
  for(z=2,lim2, lim3=sqrtnint(lim, z); for(C=2,lim3, listput(Lc, C^z)) );
  lcn=#Lc; if(lcn==0, return(-1));
  for(i=1,lcn, for(j=i,lcn, czn=Lc[i]+Lc[j]; if(czn>lim, next);
  La=findinlista(Lc,czn); lan=#La; if(!lan, next);
  for(k=1,lan, listput(Lb,czn)))); lcn=#Lb; listsort(Lb,0);
  if(bflag, for(i=1,lcn,print(i," ",Lb[i]))); if(!bflag, return(Vec(Lb)));
  }
  findinlista(list,item,sind=1)={my(ln=#list,Li=List());
  if(ln==0||sind<1||sind>ln, return(Li));
  for(i=sind,ln, if(list[i]==item, listput(Li,i))); return(Li);
  } \\ _Anatoly E. Voevudko_, Nov 23 2015
  

PARI Scripts and/or links, etc., are mine

A076467

  A076467(lim)={my(L=List(1),lim2=logint(lim,2)); for(k=3,lim2,
  for(m=2,sqrtnint(lim, k),listput(L,m^k))); listsort(L,1); print(L)}
  \\ _Anatoly E. Voevudko_, Sep 29 2015

  b076467(lim)={my(L=List(1), lim2=logint(lim, 2), ln);for(k=3,lim2,
  for(m=2, sqrtnint(lim, k), listput(L, m^k)));
  listsort(L, 1); ln=#L; for(i=1, ln, print(i , " ", L[i]))}
  \\ _Anatoly E. Voevudko_, Sep 29 2015

A000290

  b000290(maxn)={for(n=0, maxn, print(n, " ", n^2); )} \\ _Anatoly E. Voevudko_, Nov 11 2015

  a000290(maxn)={my(i=0);for(n=0,maxn,i++;print(i," ",n^2," = ",n,"^2"))}
  \\ _Anatoly E. Voevudko_, Nov 11 2015 [not published]
  A000290(maxn)={my(V=Vec(0,maxn+1));for(n=0, maxn, V[n+1]=n^2); print(V)}
  \\ _Anatoly E. Voevudko_, Nov 11 2015 [not published]

A003309

  A003309(maxn,bflag=0)={my(Vw=vector(maxn, x, x+1),Vl=Vec([1]),vwn=#Vw,i);
  while(vwn>0, i=Vw[1]; Vl=concat(Vl,[i]);
        Vw=vector((vwn*(i-1))\i,x,Vw[(x*i+i-2)\(i-1)]); vwn=#Vw);
  if(bflag, vwn=#Vl; for(i=1,vwn, print(i," ",Vl[i]))); if(!bflag, return(Vl)); }
  \\ _Anatoly E. Voevudko_, Feb 28 2016 [Not published. Upgraded Phil Carmody script from A003309]

  A003309(maxn,bflag=0)={my(Vw=vector(maxn, x, x+1),Vl=Vec([1]),vwn=#Vw,i,fn="b003309.txt");
  while(vwn>0, i=Vw[1]; Vl=concat(Vl,[i]);
        Vw=vector((vwn*(i-1))\i,x,Vw[(x*i+i-2)\(i-1)]); vwn=#Vw);
  if(bflag, vwn=#Vl; for(i=1,vwn, write(fn,concat([i," ",Vl[i]])))); if(!bflag, return(Vl)); }
  \\ _Anatoly E. Voevudko_, Feb 28 2016 [Not published. Same as above, but writing ready to use b-file (instead of printing).]

A192503

  A192503(maxn,bflag=0)={my(Vw=vector(maxn, x, x+1),Vl=Vec([1]),vwn=#Vw,i,vj,L=List());
  while(vwn>0, i=Vw[1]; Vl=concat(Vl,[i]);
        Vw=vector((vwn*(i-1))\i,x,Vw[(x*i+i-2)\(i-1)]); vwn=#Vw);
  kill(Vw); vwn=#Vl;
  for(j=1,vwn, vj=Vl[j]; if(isprime(vj),listput(L,vj))); kill(Vw); vwn=#L;
  if(bflag, for(i=1,vwn, print(i," ",L[i]))); if(!bflag, return(Vec(L)));
  } \\ _Anatoly E. Voevudko_, Feb 28 2016

A192504

  A192504(maxn,bflag=0)={my(Vw=vector(maxn, x, x+1),Vl=Vec([1]),vwn=#Vw,i,vj,L=List());
  while(vwn>0, i=Vw[1]; Vl=concat(Vl,[i]);
        Vw=vector((vwn*(i-1))\i,x,Vw[(x*i+i-2)\(i-1)]); vwn=#Vw);
  kill(Vw); vwn=#Vl;
  for(j=1,vwn, vj=Vl[j]; if(!isprime(vj),listput(L,vj))); kill(Vw); vwn=#L;
  if(bflag, for(i=1,vwn, print(i," ",L[i]))); if(!bflag, return(Vec(L)));
  } \\ _Anatoly E. Voevudko_, Feb 28 2016

A008882, A008883, A008884

  Collatz(n,lim=0)={
  my(c=n,e=0,L=List(n)); if(lim==0, e=1; lim=n*10^6); 
  for(i=1,lim, if(c%2==0, c=c/2, c=3*c+1); listput(L,c); if(e&&c==1, break));
  return(Vec(L)); } \\_Anatoly E. Voevudko_, Mar 26 2016
  
  print(Collatz(99,60)) \\A008882 (from 99 first 60)
  print(Collatz(51,60)) \\A008883 (from 51 first 60)
  print(Collatz(27)) \\A008884 (from 27 to the first 1)

A070165

  Collatz(n,lim=0)={
  my(c=n,e=0,L=List(n)); if(lim==0, e=1; lim=n*10^6); 
  for(i=1,lim, if(c%2==0, c=c/2, c=3*c+1); listput(L,c); if(e&&c==1, break));
  return(Vec(L)); } \\_Anatoly E. Voevudko_, Mar 26 2016
  Collatzaf(ns,nf,fn="")={
  my(V,vn,fng="aCollatz.txt"); if(fn=="", fn=fng);
  for(i=ns,nf, V=Collatz(i); vn=#V; write(fn,Str(i,"/",vn,": ",V)); kill(V); );
  } \\_Anatoly E. Voevudko_, Mar 30 2016
  Collatzaf(1,10000,"a070165.txt"); \\ a070165.txt file with 10000 rows