|
%I #35 Jan 06 2023 19:50:15
%S 4,1,5,8,9,11,12,13,16,17,20,22,23,25,28,30,31,34,36,39,41,43,44,46,
%T 47,48,50,51,52,53,55,56,57,58,61,62,64,65,66,67,68,71,72,74,75,76,77,
%U 79,80,81,83,84,85,86,89,90,91,93,94,96,97,99,100,103,105
%N Sequence describing the positions of the consonants when the sequence is written in French (version 1).
%C The sequence depends on how the initial terms are selected - see the Angelini article for the precise definition.
%H Eric Angelini, <a href="https://www.pourlascience.fr/sr/article/jeux-de-suites-1637.php">Jeux de suites</a>, in Dossier Pour La Science, pp. 32-35, Volume 59 (Jeux math'), April/June 2008, Paris.
%o (PARI) lista(n) = frconspos([4])(n)
%o frconspos(startseq) = n -> {
%o my(consonants = Vec("bcdfghjklmnpqrstvwxyz"), charoffset = 0, seqindex = #startseq + 1, seq = Vec(startseq, n));
%o if(n <= #startseq, seq[1..n], for(i = 1, n, my(thisfrench = Vec(strchr([c|c<-Vecsmall(French(seq[i])), c>96])), consonantindexes = select(x->setsearch(consonants, x), thisfrench, 1));
%o for(j = 1, #consonantindexes, if(#select(I -> I == consonantindexes[j] + charoffset, startseq) == 0, if(seqindex + j - 1 <= n, seq[seqindex] = consonantindexes[j] + charoffset; seqindex++)));
%o charoffset += #thisfrench);
%o seq)}
%o /* see A167507 for French() */
%o \\ _Tyler Busby_, Dec 31 2022
%Y Cf. A139212, A139357, A139358, A167507.
%K nonn,word,easy
%O 1,1
%A _N. J. A. Sloane_ (based on Angelini's article), Jun 08 2008
%E a(12) corrected and more terms from _Tyler Busby_, Dec 31 2022
|
