login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A227644 Perfect powers equal to the sum of 2 factorial numbers. 10
4, 8, 25, 121, 144, 5041 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
a(7), if it exists, is greater than 10^100.
a(7), if it exists, is greater than 10000!. - Filip Zaludek, Jul 18 2017
a(7), if it exists, is greater than 11750!. - Filip Zaludek, Sep 07 2018
a(7), if it exists, is greater than 20000!. - Filip Zaludek, Nov 04 2020
LINKS
EXAMPLE
5041 = 71^2 = 1! + 7!.
PROG
(C)
/* To compile: gcc -Wall -O2 A227644.c -o A227644 -lgmp */
#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>
int main()
{
int bsz=256, a=0;
mpz_t *f, t;
f = malloc(sizeof(mpz_t) * bsz);
mpz_init(t); mpz_init(f[0]); mpz_set_ui(f[0], 1);
while (1)
{
a += 1;
if (a == bsz)
{
bsz *= 2;
f = (mpz_t *) realloc(f, sizeof(mpz_t) * bsz);
}
mpz_init(f[a]);
mpz_mul_ui(f[a], f[a-1], a);
for (int i=1; i<=a; i++)
{
mpz_add(t, f[a], f[i]);
if (mpz_perfect_power_p(t))
{
gmp_printf("%Zd, ", t);
fflush(stdout);
}
}
}
return 0;
}
CROSSREFS
Sequence in context: A046736 A174171 A262042 * A074188 A126733 A297695
KEYWORD
nonn,more
AUTHOR
Giovanni Resta, Jul 19 2013
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 09:04 EDT 2024. Contains 371240 sequences. (Running on oeis4.)