#在GP/PARI CALCULATOR版本2.9.2(已发布)和运行mingw(x86-64/GMP-6.1.2内核)64位版本的d64中计算2017-10-23,编译时间:2017年3月22日,gcc版本4.9.1(gcc),线程引擎:单#通过以下PARI计划:##v163820=readvec(“…/a163820_terms_only.txt”)#A163820(n)=v163820【n】;#A293902(n)={my(exp_combos=功率集(系数(n)[,2]),m=1);对于(i=1,#exp_comos,m*=vecproduct(exp_combos[i])!);m;};#vecproduct(v)={my(m=1);对于(i=1,#v,m*=v[i]);m;};#功率集(v)={my(siz=2^长度(v),pv=矢量(siz));对于(i=0,siz-1,pv[i+1]=选择位(v,i));pv;};#choosebit(v,m)={my(s=向量(hammingweight(m)),i=j=1);而(m>0,如果(m%2,s[j]=v[i];j++);i++;m>>=1),s;};#A293900(n)=A163820(n)/A293902(n);#对于(n=1179,写入(“b293900.txt”,n,“”,A293900(n));##文件a163820_terms_only.txt是这样准备的:#gawk“{print$2}”<b163820_upto179_from_David_Corneth.txt>a163820_terms_only.txt#1 02 13 14 15 16 27 18 19 110 211 112 913 114 215 216 117 118 919 120 921 222二23 124 4025 126 227 128 929 130 34831 132 133 234 235 236 11037 138 239 240 4041 142 34843 144 9个45 9个46 247第1页48 17549 150 951 252 953 154 4055 2个56 4057 258 259 160 13866061 162 263 9个64 165 266 34867 168 969 270 34871 172 112773 174 275 976 977 278 34879 180 17581 182 283 184 13866085 286 287 288 4089第1页90 13866091 292 993 294 2个95 2个96 75697 198 999 9100 110101 1102 348103 1104 40105 348106 2107 1108 1127109 1110 348111 2112 175113 1114 348115 2116 9117 9118 2119 2120 50107680121 1122 2123 2124 9125 1126 138660127 1128 1129 2130 348131 1132 138660133 2134 2135 40个136 40个137 1138 348139 1140 138660141 2142 2143 2144 10530145 2146 2147 9148 9149个1150 138660151 1152 40153 9154 348155 2156 138660157 1158 2159 2个160 756161 2162 175163 1164 9165 348166 2167 1168 50107680169 1170 348171 9172 9173 1174 348175 9176 175177 2个178 2179 1