OFFSET
1,3
COMMENTS
The German name of n ends in a vowel iff n==2 or n==3 (mod 100), so the German analog of this sequence would be {2, 3, 102, 103, 202, 203, ...}. - M. F. Hasler, Sep 20 2014
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (1,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,0,0,0,0,0,0,0,0,0,0,0,0,1,-1).
FORMULA
From Chai Wah Wu, Apr 18 2024: (Start)
a(n) = a(n-1) + a(n-46) - a(n-47) for n > 48.
G.f.: x^2*(x^46 + 4*x^45 + 2*x^44 + x^43 + x^42 + 2*x^41 + 4*x^40 + 2*x^39 + x^38 + x^37 + 2*x^36 + 4*x^35 + 2*x^34 + x^33 + x^32 + 2*x^31 + 4*x^30 + 2*x^29 + x^28 + x^27 + 2*x^26 + 4*x^25 + 2*x^24 + x^23 + x^22 + 2*x^21 + 4*x^20 + 2*x^19 + x^18 + x^17 + 2*x^16 + 4*x^15 + 2*x^14 + x^13 + x^12 + 2*x^11 + 4*x^10 + 2*x^9 + x^8 + x^7 + 9*x^6 + 3*x^5 + 4*x^4 + 2*x^3 + x^2 + x + 1)/(x^47 - x^46 - x + 1). (End)
PROG
(PARI) print1(0", "); for(n=1, 3, for(k=10^(n-1), 10^n-1, m100=k%100; m10=k%10; if(m100==12 || ((m100>20 || m100<10) && ((m10==1)||(m10==2)||(m10==3)||(m10==5)||(m10==9))), print1(k", ")))) \\ Herman Jamke (hermanjamke(AT)fastmail.fm), Nov 20 2006
(PARI) is(n, S=[1, 2, 3, 5, 9, 12])=setsearch(S, if(n%100>19, n%10, n%100))||!n \\ Add 0 in S (and optionally remove trailing ||...) for the A059437 variant. Add "Set()" around [...] for PARI versions < 2.6. - M. F. Hasler, Sep 14 2014
CROSSREFS
KEYWORD
word,nonn,easy
AUTHOR
EXTENSIONS
Corrected by Logan K. Young (lkyfella(AT)yahoo.com), Mar 19 2005
More terms from Herman Jamke (hermanjamke(AT)fastmail.fm), Nov 20 2006
STATUS
approved