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!)
A274518 Numbers k such that k^2 divides A000073(k). 1
1, 103, 112, 2621, 30576, 77168, 694512, 9919728, 24770928, 55638128, 57268848, 80995824, 1300820976 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
LINKS
EXAMPLE
tribonacci(103) = 331800673921785084815380861 = 103^2 * 31275395788649739354829.
PROG
(Ruby)
require 'matrix'
def power(a, n, mod)
return Matrix.I(a.row_size) if n == 0
m = power(a, n >> 1, mod)
m = (m * m).map{|i| i % mod}
return m if n & 1 == 0
(m * a).map{|i| i % mod}
end
def f(m, n, mod)
ary0 = Array.new(m, 0)
ary0[0] = 1
v = Vector.elements(ary0)
ary1 = [Array.new(m, 1)]
(0..m - 2).each{|i|
ary2 = Array.new(m, 0)
ary2[i] = 1
ary1 << ary2
}
a = Matrix[*ary1]
(power(a, n, mod) * v)[m - 1]
end
p (1..10 ** 6).select{|i| f(3, i, i * i) == 0}
CROSSREFS
Sequence in context: A178527 A144714 A140817 * A066131 A328805 A095639
KEYWORD
nonn,more
AUTHOR
Seiichi Manyama, Jun 26 2016
EXTENSIONS
a(9)-a(13) from Chai Wah Wu, Jan 29 2018
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 29 08:13 EDT 2024. Contains 371265 sequences. (Running on oeis4.)