#include #define ll long long #define Nmax 1000001 using namespace std; int n,P[Nmax],nrp; ll x,sav,nrb,res,y; int main() { cin>>n; for (int i=2;i>x; y = 1; sav = x; res += sav; for (int c=1;c<=nrp && x!=1;c++) { while (x%P[c]==0) { x/=P[c]; y*=P[c]; res += sav/y; } } if (x!=1) res++; } cout<