2021-03-29から1日間の記事一覧
maspypy.com これの解読をする aが乗法的な時のDirichlet 積の計算 c=b; for(int p:primes) for(int i=N/p;i;i--){ int n=p*i,q=p,m=i; while(true){ c[n]+=a[q]*c[m]; if(m%p)break; q*=p;m/=p; } } pについてのループが終わると、 が成立する様になること…
maspypy.com これの解読をする aが乗法的な時のDirichlet 積の計算 c=b; for(int p:primes) for(int i=N/p;i;i--){ int n=p*i,q=p,m=i; while(true){ c[n]+=a[q]*c[m]; if(m%p)break; q*=p;m/=p; } } pについてのループが終わると、 が成立する様になること…