- num = 500; a = Array[20 &, num];
- a[[1]] = a[[2]] = 1;
- Do[Do[s = a[[k]] + a[[n - k]]; If[s < a[[n]], a[[n]] = s];
- If[Mod[n, k] == 0, s = a[[k]] + a[[n/k]];
- If[s < a[[n]], a[[n]] = s]], {k, n/2}], {n, 3, num}];
- Table[{k, Select[Range@num, a[[#]] == k &]}, {k, 16}]
复制代码
{{1,{1,2}},{2,{3,4}},{3,{5,6,8}},{4,{7,9,10,12,16}},{5,{11,13,14,15,17,18,20,24,32}},{6,{19,21,22,25,26,27,28,30,33,34,36,40,48,64}},{7,{23,29,31,35,37,38,39,41,42,44,45,49,50,51,52,54,56,60,65,66,68,72,80,96,128}},{8,{43,46,47,53,55,57,58,61,62,63,67,69,70,73,74,75,76,78,81,82,84,85,88,90,97,98,99,100,102,104,108,112,120,129,130,132,136,144,160,192,256}},{9,{59,71,77,79,83,86,87,89,91,92,93,94,95,101,103,105,106,109,110,111,113,114,116,117,119,121,122,123,124,125,126,131,133,134,135,137,138,140,145,146,147,148,150,152,153,156,161,162,164,165,168,170,176,180,193,194,195,196,198,200,204,208,216,224,240,257,258,260,264,272,288,320,384}},{10,{107,115,118,127,139,141,142,143,149,151,154,155,157,158,159,163,166,167,169,171,172,174,175,177,178,181,182,183,184,185,186,187,188,189,190,197,199,201,202,205,206,207,209,210,212,217,218,219,220,221,222,225,226,228,231,232,234,238,241,242,243,244,245,246,248,250,252,255,259,261,262,265,266,268,270,273,274,276,280,289,290,291,292,294,296,297,300,304,306,312,321,322,324,325,328,330,336,340,352,360,385,386,387,388,390,392,396,400,408,416,432,448,480}},{11,{173,179,191,203,211,213,214,215,223,227,229,230,233,235,236,237,239,247,249,251,253,254,263,267,269,271,275,277,278,279,281,282,284,285,286,287,293,295,298,299,301,302,303,305,307,308,309,310,313,314,315,316,318,323,326,327,329,331,332,333,334,335,337,338,339,341,342,343,344,345,348,350,351,353,354,356,357,361,362,363,364,365,366,368,369,370,372,374,375,376,378,380,389,391,393,394,397,398,399,401,402,404,405,409,410,411,412,414,417,418,420,424,425,429,433,434,435,436,438,440,441,442,444,449,450,452,455,456,459,462,464,468,476,481,482,483,484,485,486,488,490,492,495,496,500}},{12,{283,311,317,319,346,347,349,355,358,359,367,371,373,377,379,381,382,395,403,406,407,413,415,419,421,422,423,426,427,428,430,431,437,439,443,445,446,447,451,453,454,457,458,460,461,463,465,466,469,470,471,472,474,475,477,478,487,489,491,493,494,497,498}},{13,{383,467,473,479,499}},{14,{}},{15,{}},{16,{}}} |