Ralf Hemmecke <[EMAIL PROTECTED]> writes:
> I made it two functions for some reasons.
Thanks for all the good reasons.
I'm off to Lyon now, I just checked in what I have so far (it's not much,
though).
I checked Subst from combstruct, but I couldn't really find out what it
does. It seems that the arguments are in different order, and the results do
not always coincide. Maybe you can find out more:
with(combstruct);
f := {B=Union(Z, Prod(B, B)), F = Subst(Set(Z, card >= 1), B), T = Subst(B,
Set(Z, card >= 1))};
> seq(count([F, f, labeled], size=i), i=1..10);
1, 3, 19, 207, 3211, 64383, 1581259, 45948927, 1541641771, 58645296063
-- this agrees with Compose(B, NonEmpty SetSpecies)
> seq(count([F, f, unlabeled], size=i), i=1..10);
1, 2, 5, 15, 51, 188, 731, 2950, 12235, 51822
-- this agrees with Compose(B, NonEmpty SetSpecies)
> seq(count([T, f, labeled], size=i), i=1..10);
1, 3, 19, 193, 2721, 49171, 1084483, 28245729, 848456353, 28875761731
-- this agrees with Compose(SetSpecies, B)
> seq(count([T, f, unlabeled], size=i), i=1..10);
1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796
-- surprise! I expected
-- 1,1,2,4,10,26,77,235,758,2504
Martin
-------------------------------------------------------------------------
Take Surveys. Earn Cash. Influence the Future of IT
Join SourceForge.net's Techsay panel and you'll get the chance to share your
opinions on IT & business topics through brief surveys-and earn cash
http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV
_______________________________________________
Aldor-combinat-devel mailing list
Aldor-combinat-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/aldor-combinat-devel
