从 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);