A046074 -id:A046074 - cp's OEIS frontend

A226064 Largest fixed point in base n for the sum of the fourth power of its digits.

1, 1, 243, 419, 1, 1, 273, 1824, 9474, 10657, 1, 8194, 1, 53314, 47314, 36354, 1, 246049, 53808, 378690, 170768, 185027, 1, 247507, 1, 1002324, 722739, 278179, 301299, 334194, 1004643, 959859, 1, 1538803, 1798450, 1, 4168450, 2841074, 1, 1877793, 5556355

Offset: 2

Kevin L. Schwartz and Christian N. K. Anderson, May 24 2013

All fixed points in base n have at most 5 digits. Proof: In order to be a fixed point, a number with d digits in base n must meet the condition n^d <= d*(n-1)^4, which is only possible for d < 5.

For 5-digit numbers vwxyz in base n, only numbers where v*n^4 + n^3 - 1 <= v^4 + 3*(n-1)^4 or v*n^4 + n^4 - 1 <= v^4 + 4*(n-1)^4 are possible fixed points. v <= 2 for n <= 250.

			The fixed points in base 8 are {1,16,17,256,257,272,273}, because in base 8, these are written as {1,20,21,400,401,420,421} and 1^4 = 1, 2^4 + 0^4 = 16, 2^4 + 1^4 = 17, 4^4 + 0^4 + 0^4 = 256, etc. The largest of these is 273 = a(8).

Christian N. K. Anderson, Table of n, a(n) for n = 2..250
Christian N. K. Anderson, Table of base, largest fixed point, number of fixed points, and a list of all fixed points in base 10 and base n for n = 1..250

Cf. A226063 (number of fixed points).
Cf. A052455 (fixed points in base 10).
Cf. A023052, A046074, A046197.

for(b in 2:50) {
    fp=c()
    for(w in 1:b-1) for(x in 1:b-1) if((v1=w^4+x^4)<=(v2=w*b^3+x*b^2))
        for(y in 1:b-1) if((u1=v1+y^4)<=(u2=v2+y*b) & u1+b^4>u2+b-1) {
            z=which(u1+(1:b-1)^4==u2+(1:b-1))-1
            if(length(z)) fp=c(fp,u2+z)
        }
    cat("Base",b,":",fp[-1],"\n")
}

A192636 Powerful sums of two powerful numbers.

Original entry on oeis.org

8, 9, 16, 25, 32, 36, 64, 72, 81, 100, 108, 121, 125, 128, 144, 169, 196, 200, 216, 225, 243, 256, 288, 289, 324, 343, 361, 392, 400, 432, 441, 484, 500, 512, 576, 625, 648, 675, 676, 729, 784, 800, 841, 864, 900, 961, 968, 972, 1000, 1024, 1089, 1125, 1152, 1156, 1225

Offset: 1

Charles R Greathouse IV, Jul 06 2011

nonn

Browning & Valckenborgh conjecture that a(n) ~ kn^2 with k approximately 0.139485255. See their Conjecture 1 and equation (14). Their Theorems 1 and 2 establish upper and lower asymptotic bounds.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5000 from Charles R Greathouse IV)
Tim D. Browning and K. Van Valckenborgh, Sums of three squareful numbers, Experimental Mathematics, Vol. 21, No. 2 (2012), pp. 204-211; arXiv preprint, arXiv:1106.4472 [math.NT], 2011.

Subsequence of A001694 and of A076871.

Cf. A001694, A007532, A005934, A005188, A003321, A014576, A023052, A046074, A013929, A076871, A143813. - Jonathan Vos Post, Jul 10 2011

With[{m = 1225}, pow = Select[Range[m], # == 1 || Min[FactorInteger[#][[;; , 2]]] > 1 &]; Intersection[pow, Plus @@@ Tuples[pow, {2}]]] (* Amiram Eldar, Feb 12 2023 *)

isPowerful(n)=if(n>3,vecmin(factor(n)[,2])>1,n==1)
sumset(a,b)={
  my(c=vectorsmall(#a*#b));
  for(i=1,#a,
    for(j=1,#b,
      c[(i-1)*#b+j]=a[i]+b[j]
    )
  );
  vecsort(c,,8)
}; selfsum(a)={
  my(c=vectorsmall(binomial(#a+1,2)),k);
  for(i=1,#a,
    for(j=i,#a,
      c[k++]=a[i]+a[j]
    )
  );
  vecsort(c,,8)
};
list(lim)={
  my(v=select(isPowerful, vector(floor(lim),i,i)));
  select(n->n<=lim && isPowerful(n), Vec(selfsum(v)))
};

Numbers k such that there exists some a, b, c with A001694(a) + A001694(b) = k = A001694(c).

Corrected (on the advice of Donovan Johnson) by Charles R Greathouse IV, Sep 25 2012

cp's OEIS Frontend

A226064 Largest fixed point in base n for the sum of the fourth power of its digits.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

R

A192636 Powerful sums of two powerful numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions