技巧 筛1~n的所有因子

从 i : 1~n, 是i的倍数， 则计入该数

复杂度 n*(1/1+1/2+1/3+...1/n)=nlogn

ll d[N]; // 计每个数的因子数
set<ll> s[N]; // 存入每个数的所有因子
for(i=1; i<=n; i++) 
    for(j=i; j<=n; j+=i) // <n的i的倍数  n/i个
        d[j]++, s[j].insert(i);