"Integer partitions"의 두 판 사이의 차이
37번째 줄: | 37번째 줄: | ||
− | <h5>partition and | + | <h5>partition rank and crank</h5> |
− | (*define a integer you want to investigate*) | + | (*define a integer you want to investigate*)n := 6<br> (*choose the proper moduli for the partition statistics*)<br> md := 11<br> S[n_] := IntegerPartitions[n]<br> (*define the rank of a partition with the name "pr"*)<br> pr[s_] := Max[s] - Length[s]<br> (*define the crank of a partition with the name "crank"*)<br> Om[s_] := Count[s, 1]<br> KK[s_] := Select[s, # > Om[s] &]<br> Mu[s_] := Length[KK[s]]<br> crank[s_] := If[Om[s] == 0, Max[s], Mu[s] - Om[s]]<br> (*modulus distribution of partition rank*)<br> Sort[Tally[Table[Mod[pr[s], md], {s, S[n]}]]]<br> (*modulus distribution of partition crank*)<br> Sort[Tally[Table[Mod[crank[s], md], {s, S[n]}]]]<br> (*list of paritions with rank & crank *)<br> Do[Print[s, ", rank=", pr[s], "\[Congruent]", Mod[pr[s], md], "(mod ",<br> md, ")", ", crank=", crank[s], "\[Congruent]", Mod[crank[s], md],<br> "(mod ", md, ")"], {s, S[n]}]<br> (*you will see the distribution of rank/crank modulus,the partition \<br> statistics and list of paritions with rank&crank*) |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
2010년 3월 4일 (목) 20:45 판
background
n:=9
md:=5
n:=12
md:=7
n:=6
md:=11
will be a good choice
\(p(5k+4)\equiv 0 \pmod 5\)
\(p(7k+5)\equiv 0 \pmod 7\)
\(p(11k+6)\equiv 0 \pmod {11}\)
partition rank and crank
(*define a integer you want to investigate*)n := 6
(*choose the proper moduli for the partition statistics*)
md := 11
S[n_] := IntegerPartitions[n]
(*define the rank of a partition with the name "pr"*)
pr[s_] := Max[s] - Length[s]
(*define the crank of a partition with the name "crank"*)
Om[s_] := Count[s, 1]
KK[s_] := Select[s, # > Om[s] &]
Mu[s_] := Length[KK[s]]
crank[s_] := If[Om[s] == 0, Max[s], Mu[s] - Om[s]]
(*modulus distribution of partition rank*)
Sort[Tally[Table[Mod[pr[s], md], {s, S[n]}]]]
(*modulus distribution of partition crank*)
Sort[Tally[Table[Mod[crank[s], md], {s, S[n]}]]]
(*list of paritions with rank & crank *)
Do[Print[s, ", rank=", pr[s], "\[Congruent]", Mod[pr[s], md], "(mod ",
md, ")", ", crank=", crank[s], "\[Congruent]", Mod[crank[s], md],
"(mod ", md, ")"], {s, S[n]}]
(*you will see the distribution of rank/crank modulus,the partition \
statistics and list of paritions with rank&crank*)
various partitions
(* partitions with at most 5 parts *)
IntegerPartitions[7, 5]
(* partition into exactly three parts *)
VS[n_] := IntegerPartitions[n, {3}]
VS[11]
(* number of partitions into distinct parts *)
PartitionsQ[11]
(* partition into odd parts *)
IntegerPartitions[11, All, {1, 3, 5, 7, 9, 11}]