A "doublely-recursive" sequence
Leroy Quet
qqquet at mindspring.com
Sun Aug 3 00:34:51 CEST 2003
Farideh Firoozbakht wrote:
>...
>for n=0,...,100
>
> a(n)= :
> 1,1,2,4,1,2,4,9,2,4,9,30,15,4,9,30,97,84,84,26,15,127,
> 308,30,15,127,898,24,913,97,24,913,308,69,2,4,9,30,2996,
> 4217,308,560,97,69,1040,11,69,868,9,30,2996,7327,14566,
> 13618,39,11544,26,1938,913,2875,64695,4292,97,4288,1,2,4,9,
> 30,2996,64721,11852,13618,11852,84,30,2996,211582,13698,3909,
> 146861,26,1938,2972,4217,39494,433,676711,481203,146861,223434,
> 2972,2875,676711,11,69,146861,8,26,1938,7327
>
>and {n,a(n)}=:
>
> {0, 1},{1, 1},{2, 2},{3, 4},{4, 1},{5, 2},{6, 4},{7, 9},{8, 2},{9,4},
> {10, 9},{11, 30},{12, 15},{13, 4},{14, 9},{15, 30},{16, 97},{17, 84},
> {18, 84},{19, 26},{20, 15},{21, 127},{22, 308},{23, 30},{24, 15},
> {25, 127},{26, 898},{27, 24},{28, 913}, {29, 97},{30, 24},{31, 913},
> {32, 308},{33, 69},{34, 2},{35, 4},{36, 9},{37, 30},{38, 2996},
> {39, 4217},{40, 308},{41, 560},{42, 97},{43, 69},{44, 1040},{45, 11},
> {46, 69},{47, 868},{48, 9},{49, 30},{50, 2996},{51, 7327},{52, 14566},
> {53, 13618},{54, 39}, {55, 11544},{56, 26},{57, 1938},{58, 913},
> {59, 2875}, {60,64695},{61, 4292},{62, 97},{63, 4288},{64, 1},{65, 2},
> {66, 4},{67, 9},{68, 30},{69, 2996},{70, 64721}, {71, 11852},{72,13618}
> {73, 11852},{74, 84},{75, 30},{76, 2996},{77, 211582},{78, 13698},
> {79, 3909},{80, 146861},{81, 26},{82, 1938},{83, 2972}, {84, 4217},
> {85, 39494},{86, 433},{87, 676711},{88, 481203},{89, 146861},
> {90, 223434},{91, 2972},{92, 2875},{93, 676711}, {94, 11},{95, 69},
> {96, 146861},{97, 8},{98, 26},{99, 1938},{100, 7327}
>
>
>If the sequence has closed form, finding of it,is difficult.
>
>
>Regards,
>
>Farideh Firoozbakht
>UniverSity of Isfahan,Iran
>f.firoozbakht at sci.ui.ac.ir
>
>
>Quoting Leroy Quet <qqquet at mindspring.com>:
>>>...
>> If I did not make a mistake calculating the first few terms by hand,
>> here is a recursively-defined sequence which is not in the EIS yet.
>>
>> a[0] = 1;
>>
>> and for m >= 1,
>>
>> a[m] = sum{j=0 to a[m-1](mod m)} a[j]
>>
>> Ascii-art:
>> a[m-1](mod m)
>> ---
>> \
>> a[m] = / a[k]
>> ---
>> k=0
>>
>>
>> And, 0 <= a[m-1](mod m) <= m-1.
>>
>>...
Interesting.
I am slightly curious about such things such as:
for which indexes k, does a[k] = 1,
or when is a[k] < a[k-1] ?
And what would be the sequence of differences between consecutive k's,
where
a[k] < a[k-1]?
(The k's where a[k] = 1 are 0, 1, 4, and 64, for k < 100 {as seen from
the list generated by Farideh Firoozbakht}. These k's, except the 0, all
happen to be powers of 2; but this is much too small a sample to deduce
anything, obviously.)
Hmmm...I also wonder what the plot of this sequence, perhaps scaled
logarithmically, looks like.
Thank you,
Leroy Quet
More information about the SeqFan
mailing list