login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A289834 Number of perfect matchings on n edges which represent RNA secondary folding structures characterized by the Lyngso and Pedersen (L&P) family and the Cao and Chen (C&C) family. 0
1, 1, 3, 11, 39, 134, 456, 1557, 5364, 18674, 65680, 233182, 834796, 3010712, 10929245, 39904623, 146451871, 539972534, 1999185777, 7429623640, 27705320423, 103636336176, 388775988319, 1462261313876, 5513152229901, 20832701135628, 78884459229627 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Table of n, a(n) for n=0..26.

Aziza Jefferson, The Substitution Decomposition of Matchings and RNA Secondary Structures, PhD Thesis, University of Florida, 2015.

FORMULA

a(n) = Sum_{i=0..n-2} C_i*(Sum_{j=1..n-i} C_j - (n-i)) + C_n where C is A000108.

From Vaclav Kotesovec, Jul 13 2017: (Start)

Recurrence (of order 3): (n+2)*(41*n^3 - 228*n^2 + 391*n - 180)*a(n) = 6*(41*n^4 - 187*n^3 + 192*n^2 + 120*n - 160)*a(n-1) - 3*(3*n - 4)*(41*n^3 - 146*n^2 + 83*n + 70)*a(n-2) + 2*(2*n - 5)*(41*n^3 - 105*n^2 + 58*n + 24)*a(n-3).

a(n) ~ 41 * 4^n / (9*sqrt(Pi)*n^(3/2)).

(End)

MAPLE

a:= proc(n) option remember; `if`(n<4, [1$2, 3, 11][n+1],

      (2*(74*n^2-69*n-110)*a(n-1)-3*(89*n^2-139*n-70)*a(n-2)+

       2*(91*n^2-204*n-52)*a(n-3)-4*(5*n+1)*(2*n-7)*a(n-4))

       /((n+2)*(23*n-43)))

    end:

seq(a(n), n=0..40);  # Alois P. Heinz, Jul 13 2017

PROG

(Python)

class Memoize:

    def __init__(self, func):

        self.func = func

        self.cache = {}

    def __call__(self, arg):

        if arg not in self.cache:

            self.cache[arg] = self.func(arg)

        return self.cache[arg]

@Memoize

def a(n): return [1, 1, 3, 11][n] if n<4 else (2*(74*n**2 - 69*n - 110)*a(n - 1) - 3*(89*n**2 - 139*n - 70)*a(n - 2) + 2*(91*n**2 - 204*n - 52)*a(n - 3) - 4*(5*n + 1)*(2*n - 7)*a(n - 4))/((n + 2)*(23*n - 43))

print map(a, xrange(27)) # Indranil Ghosh, Jul 15 2017, after Maple code

CROSSREFS

Cf. A000108, A256334.

Sequence in context: A227638 A166336 A002783 * A007482 A134760 A257290

Adjacent sequences:  A289831 A289832 A289833 * A289835 A289836 A289837

KEYWORD

nonn

AUTHOR

Kyle Goryl, Jul 13 2017

EXTENSIONS

More terms from Alois P. Heinz, Jul 13 2017

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 27 02:54 EDT 2019. Contains 323597 sequences. (Running on oeis4.)