A298945 -id:A298945 - cp's OEIS frontend

A298946 a(n) = binomial(2*c-1, c-1) (mod c^4), where c is the n-th composite number.

35, 462, 2339, 4627, 2378, 4238, 5148, 1260, 57635, 85026, 64410, 100509, 163716, 171918, 93876, 309780, 148969, 444220, 370712, 532771, 652200, 938386, 816466, 907874, 569300, 1107298, 2470810, 2953692, 887812, 1341810, 2956584, 1941390, 589961, 6248628

Offset: 1

Felix Fröhlich, Jan 30 2018

nonn

Composites c where a(n) = 1 could be called "Wolstenholme pseudoprimes". Do any such composites exist?

A necessary condition for c to be a "Wolstenholme pseudoprime" would be that it is a term of A228562 or A267824.

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A088164, A228562, A244214, A267824, A281302, A298944, A298945.

R:= NULL:
count:= 0: F:= 10;
for n from 4 while count < 100 do
  F:= F * (4*n-2)/n;
  if not isprime(n) then
     count:= count+1;
     R:= R, F mod (n^4);
  fi
od:
R; # Robert Israel, Feb 02 2018

Table[Mod[Binomial[2 c - 1, c - 1], c^4], {c, Select[Range@ 50, CompositeQ]}] (* Michael De Vlieger, Feb 01 2018 *)

forcomposite(c=1, 200, print1(lift(Mod(binomial(2*c-1, c-1), c^4)), ", "))

from sympy import binomial, composite
def A298946(n):
    c = composite(n)
    return binomial(2*c-1,c-1) % c**4 # Chai Wah Wu, Feb 02 2018

A298944 a(n) = 2^(c-1) mod c^2, where c is the n-th composite number.

Original entry on oeis.org

8, 32, 0, 13, 12, 32, 156, 184, 0, 176, 288, 319, 464, 320, 341, 496, 40, 64, 212, 0, 301, 308, 9, 1040, 952, 472, 1088, 1544, 800, 391, 508, 2048, 1191, 1312, 922, 2608, 284, 2359, 1920, 688, 1800, 3488, 2668, 2524, 0, 2291, 428, 144, 3109, 2612, 1472, 2888

Offset: 1

Felix Fröhlich, Jan 30 2018

nonn

a(n) = 0 iff c is a term of A000079 > 4.
Composites c where a(n) = 1 could be called "Wieferich pseudoprimes". Do any such composites exist?
A necessary condition for c to be a "Wieferich pseudoprime" would be that it is a term of both A001567 and A270833 (see comments in A240719).

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A000079, A001220, A001567, A240719, A270833, A298945, A298946.

map(c -> 2&^(c-1) mod c^2, remove(isprime, [$4..1000])); # Robert Israel, Feb 27 2018

composite[n_Integer] := FixedPoint[n + PrimePi@ # + 1 &, n + PrimePi@ n + 1]; Array[With[{c = composite@ #}, Mod[2^(c - 1), c^2]] &, 52] (* Michael De Vlieger, Jan 31 2018, composite function by Robert G. Wilson v at A066277 *)

forcomposite(c=1, 200, print1(lift(Mod(2, c^2)^(c-1)), ", "))

cp's OEIS Frontend

A298946 a(n) = binomial(2*c-1, c-1) (mod c^4), where c is the n-th composite number.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python