Displaying 1-4 of 4 results found.
page
1
1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 1, 1, 7, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 7, 1, 1, 3, 1, 1, 9, 7, 1, 1, 1, 3, 3, 1, 1, 3, 1, 1, 3, 1, 1, 3
PROG
(PARI)
A253553(n) = if(n<=2, 1, my(f=factor(n), k=#f~); if(f[k, 2]>1, f[k, 2]--, f[k, 1] = precprime(f[k, 1]-1)); factorback(f));
A356300sq(x, y) = if(1==x||1==y, 1, my(lista=List([]), i, k=x, stemvec, stemlen, h=y); while(k>1, listput(lista, k); k = A253553(k)); stemvec = Vecrev(Vec(lista)); stemlen = #stemvec; while(1, if((i=vecsearch(stemvec, h))>0, return(stemvec[i])); h = A253553(h)));
The nearest common ancestor of n and gcd(n, sigma(n)) in the Doudna tree ( A005940).
+10
4
1, 1, 1, 1, 1, 6, 1, 1, 1, 2, 1, 2, 1, 2, 3, 1, 1, 2, 1, 2, 1, 2, 1, 12, 1, 2, 1, 28, 1, 6, 1, 1, 3, 2, 1, 1, 1, 2, 1, 10, 1, 3, 1, 2, 3, 2, 1, 2, 1, 1, 3, 2, 1, 2, 1, 2, 1, 2, 1, 6, 1, 2, 1, 1, 1, 3, 1, 2, 3, 2, 1, 2, 1, 2, 1, 2, 1, 3, 1, 2, 1, 2, 1, 5, 1, 2, 3, 2, 1, 2, 5, 2, 1, 2, 5, 12, 1, 1, 3, 1, 1, 3, 1, 2, 3
PROG
(PARI)
Abincompreflen(n, m) = { my(x=binary(n), y=binary(m), u=min(#x, #y)); for(i=1, u, if(x[i]!=y[i], return(i-1))); (u); };
Abinprefix(n, k) = { my(digs=binary(n)); fromdigits(vector(k, i, digs[i]), 2); };
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); (t); };
A156552(n) = {my(f = factor(n), p, p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ From A156552
A356156(n) = A348041sq(n, gcd(n, sigma(n)));
1, 1, 1, 1, 3, 3, 1, 1, 3, 3, 3, 3, 7, 3, 3, 1, 3, 9, 5, 3, 1, 3, 3, 3, 5, 3, 3, 7, 3, 9, 1, 1, 3, 3, 3, 3, 19, 3, 3, 3, 3, 3, 11, 3, 9, 3, 3, 3, 3, 3, 9, 7, 3, 9, 3, 3, 3, 3, 3, 15, 31, 3, 3, 1, 3, 9, 17, 3, 3, 3, 3, 9, 37, 3, 3, 5, 3, 21, 5, 3, 3, 3, 3, 3, 3, 3, 15, 3, 3, 45, 7, 3, 1, 3, 3, 3, 7, 3, 9, 5, 3, 9, 13, 3, 3
PROG
(PARI)
A253553(n) = if(n<=2, 1, my(f=factor(n), k=#f~); if(f[k, 2]>1, f[k, 2]--, f[k, 1] = precprime(f[k, 1]-1)); factorback(f));
A356300sq(x, y) = if(1==x||1==y, 1, my(lista=List([]), i, k=x, stemvec, stemlen, h=y); while(k>1, listput(lista, k); k = A253553(k)); stemvec = Vecrev(Vec(lista)); stemlen = #stemvec; while(1, if((i=vecsearch(stemvec, h))>0, return(stemvec[i])); h = A253553(h)));
CROSSREFS
Positions of 1's in this sequence is given by the union of A000079 and A046528.
1, 1, 1, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 9, 1, 1, 1, 1, 1, 3, 5, 1, 3, 7, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 9, 1, 1, 3, 1, 1, 3, 1, 1, 9, 1, 1, 3, 1, 1, 15, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 9, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 1, 45, 7, 1, 1, 1, 1, 3, 1, 1, 9, 5, 1, 3, 1, 1, 3
PROG
(PARI)
A253553(n) = if(n<=2, 1, my(f=factor(n), k=#f~); if(f[k, 2]>1, f[k, 2]--, f[k, 1] = precprime(f[k, 1]-1)); factorback(f));
A356300sq(x, y) = if(1==x||1==y, 1, my(lista=List([]), i, k=x, stemvec, stemlen, h=y); while(k>1, listput(lista, k); k = A253553(k)); stemvec = Vecrev(Vec(lista)); stemlen = #stemvec; while(1, if((i=vecsearch(stemvec, h))>0, return(stemvec[i])); h = A253553(h)));
Search completed in 0.053 seconds
|