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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A303731 Number of noncrossing path sets on n nodes up to rotation and reflection with each path having a prime number of nodes. 6

%I

%S 1,0,1,1,1,5,6,27,53,140,649,1297,6355,18038,63226,241741,744711,

%T 3008107,10028056,37270169,138083464,488933323,1872525356,6763888465,

%U 25498771059,95467533318,355595703773,1353873044078,5077809606803,19345857682140,73533468653115

%N Number of noncrossing path sets on n nodes up to rotation and reflection with each path having a prime number of nodes.

%H Andrew Howroyd, <a href="/A303731/b303731.txt">Table of n, a(n) for n = 0..500</a>

%o (PARI) \\ number of path sets with restricted path lengths

%o NCPathSetsModDihedral(v)={ my(n=#v);

%o my(p=serreverse(x/(1 + x*v[1] + sum(k=2, #v, (k*2^(k-3))*x^k*v[k])) + O(x^2*x^n) )/x);

%o my(vars=variables(p));

%o my(h=substvec(p + O(x^(n\2+1)),vars,apply(t->t^2, vars)));

%o my(q=x*deriv(p)/p);

%o my(R=v[1]*x + sum(i=1, (#v-1)\2, v[2*i+1]*2^(i-1)*x*(x^2*h)^i), Q=sum(i=1, #v\2, v[2*i]*2^(i-1)*(x^2*h)^i), T=intformal((p - 1 + sum(d=2,n, eulerphi(d)*substvec(q + O(x^(n\d+1)), vars, apply(t->t^d, vars))))/x));

%o O(x*x^n) + (1 + T + (1 + Q + (1+R)^2*h/(1-Q) + v[2]*x^2*h)/2)/2;

%o }

%o Vec(NCPathSetsModDihedral(vector(30, k, isprime(k))))

%Y Cf. A303729, A303732.

%K nonn

%O 0,6

%A _Andrew Howroyd_, Apr 29 2018

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 July 22 02:26 EDT 2019. Contains 325210 sequences. (Running on oeis4.)