# Computed 2017-08-27 - 28 in GP/PARI CALCULATOR Version 2.9.2 (released) amd64 running mingw (x86-64/GMP-6.1.2 kernel) 64-bit version, compiled: Mar 22 2017, gcc version 4.9.1 (GCC), threading engine: single # with the following PARI-program: # # ok(p, e, n) = { my(v=Vec(factorback(p, e))); (vecmin(v)>=0 && vecsum(v)==n); }; # A111588(n) = if(n<4, return(1)); my(x='x, f=factor((x^n-1)/(x-1)), p=f[, 1], e=2*f[, 2]~, u=vector(#e, i, [0, e[i]]), s, t); t=vecmax(e); for(i=1, #e, if(e[i]==t, u[i][2]\=2; break)); forvec(v=u, t=e-v; if(cmp(v, t)<=0 && ok(p, v, n) && ok(p, t, n), s++)); s \\ By _Charles R Greathouse IV_, Jan 19 2017 # for(n=1,239,write("b111588.txt", n, " ", A111588(n))); 1 1 2 1 3 1 4 2 5 1 6 2 7 1 8 4 9 2 10 2 11 1 12 8 13 1 14 2 15 2 16 10 17 1 18 8 19 1 20 8 21 2 22 2 23 1 24 33 25 2 26 2 27 4 28 8 29 1 30 13 31 1 32 26 33 2 34 2 35 2 36 57 37 1 38 2 39 2 40 33 41 1 42 13 43 1 44 8 45 8 46 2 47 1 48 140 49 2 50 8 51 2 52 8 53 1 54 33 55 2 56 33 57 2 58 2 59 1 60 125 61 1 62 2 63 8 64 71 65 2 66 13 67 1 68 8 69 2 70 13 71 1 72 348 73 1 74 2 75 8 76 8 77 2 78 13 79 1 80 140 81 10 82 2 83 1 84 122 85 2 86 2 87 2 88 33 89 1 90 118 91 2 92 8 93 2 94 2 95 2 96 583 97 1 98 8 99 8 100 59 101 1 102 13 103 1 104 33 105 13 106 2 107 1 108 349 109 1 110 13 111 2 112 140 113 1 114 13 115 2 116 8 117 8 118 2 119 2 120 956 121 2 122 2 123 2 124 8 125 4 126 123 127 1 128 197 129 2 130 13 131 1 132 120 133 2 134 2 135 33 136 33 137 1 138 13 139 1 140 118 141 2 142 2 143 2 144 2036 145 2 146 2 147 8 148 8 149 1 150 125 151 1 152 33 153 8 154 13 155 2 156 120 157 1 158 2 159 2 160 583 161 2 162 140 163 1 164 8 165 13 166 2 167 1 168 932 169 2 170 13 171 8 172 8 173 1 174 13 175 8 176 140 177 2 178 2 179 1 180 1676 181 1 182 13 183 2 184 33 185 2 186 13 187 2 188 8 189 33 190 13 191 1 192 2393 193 1 194 2 195 13 196 57 197 1 198 120 199 1 200 369 201 2 202 2 203 2 204 122 205 2 206 2 207 8 208 140 209 2 210 290 211 1 212 8 213 2 214 2 215 2 216 3050 217 2 218 2 219 2 220 120 221 2 222 13 223 1 224 583 225 57 226 2 227 1 228 122 229 1 230 13 231 13 232 33 233 1 234 119 235 2 236 8 237 2 238 13 239 1