login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = n^3 mod 33.
2

%I #24 Dec 18 2023 14:59:35

%S 0,1,8,27,31,26,18,13,17,3,10,11,12,19,5,9,4,29,24,28,14,21,22,23,30,

%T 16,20,15,7,2,6,25,32,0,1,8,27,31,26,18,13,17,3,10,11,12,19,5,9,4,29,

%U 24,28,14,21,22,23,30,16,20,15,7,2,6,25,32,0,1,8,27,31,26,18,13,17,3

%N a(n) = n^3 mod 33.

%C Equivalently: n^(10*m + 3) mod 33. - _G. C. Greubel_, Mar 30 2016

%H G. C. Greubel, <a href="/A070495/b070495.txt">Table of n, a(n) for n = 0..1000</a>

%H <a href="/index/Rec#order_33">Index entries for linear recurrences with constant coefficients</a>, signature (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1).

%F a(n) = a(n-33). - _G. C. Greubel_, Mar 30 2016

%t PowerMod[Range[0, 100], 3, 33] (* _G. C. Greubel_, Mar 30 2016 *)

%o (Sage) [power_mod(n,3,33) for n in range(0, 76)] # _Zerinvary Lajos_, Oct 30 2009

%o (Magma) [Modexp(n, 3, 33): n in [0..100]]; // _Vincenzo Librandi_, Mar 31 2016

%o (PARI) a(n)=n^3%33 \\ _Charles R Greathouse IV_, Apr 06 2016

%K nonn,easy

%O 0,3

%A _N. J. A. Sloane_, May 12 2002