# Computed 2018-09-11 in GP/PARI CALCULATOR Version 2.9.4 (released) amd64 running linux (x86-64/GMP-6.1.2 kernel) 64-bit version, compiled: Dec 19 2017, gcc version 7.3.0 (Ubuntu 7.3.0-1ubuntu1), threading engine: pthread # with the following program: # # allocatemem(2^30); # default(parisizemax,2^31); # strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4); # partitions_into(n, parts, from=1) = if(!n, 1, my(k = #parts, s=0); for(i=from, k, if(parts[i]<=n, s += partitions_into(n-parts[i], parts, i))); (s)); # toplevel_starting_sets(orgn, n, parts, from=1, ss=List([])) = { my(k = #parts, s=0, newss); if(lcm(Vec(ss))==orgn, s += partitions_into(n, ss)); for(i=from, k, if(parts[i]<=n, newss = List(ss); listput(newss, parts[i]); s += toplevel_starting_sets(orgn, n-parts[i], parts, i+1, newss))); (s) }; # A317624(n) = if(n<=1, 0, toplevel_starting_sets(n, n, strong_divisors_reversed(n))); # for(n=0,359,write("b317624.txt", n, " ", A317624(n))); # # Note that using the following version of partitions_into would have been somewhat faster: # # partitions_into(n, parts, from=1) = if(!n, 1, if(#parts==from, (0==(n%parts[from])), my(s=0); for(i=from, #parts, if(parts[i]<=n, s += partitions_into(n-parts[i], parts, i))); (s))); # 0 0 1 0 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 3 13 1 14 1 15 1 16 1 17 1 18 3 19 1 20 5 21 1 22 1 23 1 24 17 25 1 26 1 27 1 28 7 29 1 30 60 31 1 32 1 33 1 34 1 35 1 36 76 37 1 38 1 39 1 40 55 41 1 42 105 43 1 44 11 45 10 46 1 47 1 48 187 49 1 50 6 51 1 52 13 53 1 54 30 55 1 56 111 57 1 58 1 59 1 60 5043 61 1 62 1 63 15 64 1 65 1 66 230 67 1 68 17 69 1 70 242 71 1 72 4173 73 1 74 1 75 12 76 19 77 1 78 310 79 1 80 1007 81 1 82 1 83 1 84 15214 85 1 86 1 87 1 88 285 89 1 90 14989 91 1 92 23 93 1 94 1 95 1 96 3833 97 1 98 10 99 24 100 607 101 1 102 505 103 1 104 403 105 421 106 1 107 1 108 9550 109 1 110 533 111 1 112 2893 113 1 114 620 115 1 116 29 117 29 118 1 119 1 120 988874 121 1 122 1 123 1 124 31 125 1 126 47954 127 1 128 1 129 1 130 722 131 1 132 71946 133 1 134 1 135 341 136 697 137 1 138 885 139 1 140 70719 141 1 142 1 143 1 144 542689 145 1 146 1 147 20 148 37 149 1 150 57001 151 1 152 873 153 38 154 947 155 1 156 130169 157 1 158 1 159 1 160 36169 161 1 162 1017 163 1 164 41 165 920 166 1 167 1 168 5237011 169 1 170 1175 171 43 172 43 173 1 174 1370 175 36 176 11641 177 1 178 1 179 1 180 38271145 181 1 182 1282 183 1 184 1287 185 1 186 1555 187 1 188 47 189 721 190 1448 191 1 192 172067 193 1 194 1 195 1242 196 2557 197 1 198 238892 199 1 200 131232 201 1 202 1 203 1 204 343523 205 1 206 1 207 52 208 19369 209 1 210 17243357 211 1 212 53 213 1 214 1 215 1 216 8280086 217 1 218 1 219 1 220 332487 221 1 222 2180 223 1 224 147227 225 4233 226 1 227 1 228 516770 229 1 230 2069 231 1618 232 2059 233 1 234 438033 235 1 236 59 237 1 238 2093 239 1 240 681648139 241 1 242 21 243 1 244 61 245 42 246 2655 247 1 248 2355 249 1 250 287 251 1 252 253532171 253 1 254 1 255 2032 256 1 257 1 258 2910 259 1 260 601199 261 66 262 1 263 1 264 54754885 265 1 266 2575 267 1 268 67 269 1 270 32297774 271 1 272 43733 273 2178 274 1 275 61 276 1049712 277 1 278 1 279 71 280 50537861 281 1 282 3455 283 1 284 71 285 2498 286 2840 287 1 288 201622291 289 1 290 3206 291 1 292 73 293 1 294 527416 295 1 296 3367 297 1922 298 1 299 1 300 584480167 301 1 302 1 303 1 304 61239 305 1 306 1174588 307 1 308 984284 309 1 310 3632 311 1 312 133912123 313 1 314 1 315 685303 316 79 317 1 318 4360 319 1 320 2946237 321 1 322 3671 323 1 324 8161103 325 72 326 1 327 1 328 4141 329 1 330 164125832 331 1 332 83 333 85 334 1 335 1 336 6445165993 337 1 338 28 339 1 340 1586892 341 1 342 1776517 343 1 344 4557 345 3561 346 1 347 1 348 2507589 349 1 350 1168418 351 2728 352 933221 353 1 354 5370 355 1 356 89 357 3561 358 1 359 1