#include <bits/stdc++.h>

#define rep0(i, n) for(int i=0; i<(int)(n);i++)
#define rep1(i, n) for(int i=1; i<=(int)(n);i++)
#define all(x) (x).begin(), (x).end()
#define debug(a) if(flagdebug){cout << "debug (" #a "):" << endl << (a) << endl;}

using namespace std;

typedef uint64_t ulint;
typedef int64_t lint;

// #include <boost/multiprecision/cpp_int.hpp>
// using bint = boost::multiprecision::cpp_int;

template<class T> lint sum(vector<T> a);
template<class T> double ave(vector<T> a);
template<typename T> ostream& operator << (ostream& s, const vector<T>& v);
ulint gcd(lint a, lint b);
ulint lcm(lint a, lint b);
ulint digit(lint n, lint base);
ulint digit(lint n);
vector<lint> digitbreak(lint n, lint base);
vector<lint> digitbreak(lint n);
bool flagdebug = false;


const int MAX = 1000000;
const int MOD = 1000000007; // 998244353;


long long fac[MAX], finv[MAX], inv[MAX];

// テーブルを作る前処理
void COMinit() {
    fac[0] = fac[1] = 1; //累乗
    finv[0] = finv[1] = 1; 
    inv[1] = 1;
    for (int i = 2; i < MAX; i++){
        fac[i] = fac[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * inv[i] % MOD;
    }
}

// 二項係数計算
long long COM(int n, int k){
    if (n < k) return 0;
    if (n < 0 || k < 0) return 0;
    return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}




int main(){
//////////////////////////////////////////////////////

    flagdebug = true; ///本番時削除！！！！！

    //入力

    lint n;
    cin >> n;

    // vector<lint> a(n);
    //2次元配列
    //vector<vector<lint>> mat(row, vector<lint>(column, 初期値));

    // rep0(i, n)
    // {
    //     cin >> a[i];
    // }

    ulint ans3 = 0;

    rep1(i,n)
    {
        ulint inum = n / i;
        ulint sum = (inum + 1) * inum / 2;
        //debug(i);
        //debug(inum);
        //debug(sum);
        ans3 += sum * i;
    }

    // debug(ans3);


    // ulint ans2 = 0;

    // vector <lint> divnum(n+1, 2);

    // divnum[1] = 1;
    // for(lint j = 2; j * j <= n; j ++)
    // {
    //     for(lint i = j + j; i <= n; i += j)
    //     {
    //         divnum[i]++;
    //     }
    // }

    // debug(divnum);

    // rep1(i,n)
    // {
    //     ans2 += i * divnum[i];
    // }

    // debug(ans2);
    
    // ulint ans = 0;

    // //vector< int64_t > ret;
    // //愚直会
    // for(int64_t j = 1; j <= n; j++) 
    // {
    //     //debug(j);
    //     lint sum = 0;
    //     for(int64_t i = 1; i * i <= j; i++) {
    //         if(j % i == 0) {
    //             sum++;
    //             if(i * i != j)
    //             {
    //                 //debug(n/i);
    //                 sum ++;
    //             } 
    //         }

    //     }
    //     ans += j * sum;
    //     //debug(sum);
    //     //debug(ans);
    // }

    // debug(ans);

    
    cout << ans3 << endl;
    

    //2次元配列
    //vector<vector<lint>> mat(row, vector<lint>(column, 初期値));
    //ペア
    //pair<lint,lint> p = make_pair(1,2);
    //vector<pair<lint,lint>> c(100010, pair<lint,lint>(0,0));
    //タプル
    //tuple<lint,lint,lint> t  = make_tuple(1,2,3);


    //昇順ソート
    //sort(all(a));
    //降順ソート
    //sort(all(a),greater<lint>());
    //特殊なソート（ラムダ式）
    //sort(all(mat), [&i](const vector<lint> &a, const vector<lint> &b){return a[i] < b[i];});
    //特殊なソート（pair）
    //sort(all(a), [](const pair<lint,lint> &a, const pair<lint,lint> &b){return a.first < b.first;});
    // //[&i]は参照する変数，いらなければからの大括弧でよし．
    //ソート済集合（順序付き集合）
    //set<lint> st;
    //itr.insert(0);
    //auto itrmin = st.begin(); //最小値の次へアクセス
    //auto itrmax = st.end(); //最大値の次へアクセス（かならずデクリメントが必要）
    //itrmax --; //最大値はst.end()の1つ前
    //st.erase(0);
    //st.find;
    
    //////////////////////////////////////////////////////
    //キュー
    //queue<lint> que;
    //que.push(value);
    //que.front(); //トップを参照（消されないので消したい場合は直後にpop()を使う）
    //que.pop(); //追い出す
    //que.empty(); //空なら1，あれば0
    //優先度付きキュー
    //priority_queue<lint> que;
    //キューにque.push()で入れると最大値順に並べられる．que.pop()で取り出すときは常に最大値が前にくる．
    //最小値から取り出したい時はpriority_queue< lint, std::vector<lint>, std::greater<lint> >
    //////////////////////////////////////////////////////
    //スタック
    //stack<lint> stk;
    //stk.push(value);
    //stk.top(); //先頭を参照（消されないので消したい場合は直後にpop()を使う）
    //stk.pop(); //トップを削除
    //stk.empty(); //空なら1,あれば0
    //stk.size(); //数
    //
    // sort(input.begin() , input.end());
    // uint64_t max = numeric_limits<uint64_t>::max();

//////////////////////////////////////////////////////
// Ctrl + Opt + N to make
    return 0;
}

//vectorprint用
template<typename T> ostream& operator << (ostream& s, const vector<T>& v) {
    s << "[\t";
    for (auto i = v.begin(); i < v.end() - 1; i++){s << *i << "\t";}
    s << *(v.end() - 1) << "\t]\n";
	return s;
}

// 最大公約数
ulint gcd(lint a, lint b){
    if(b!=0){return gcd(b, a % b);}
    return a;    
}

//最小公倍数
ulint lcm(lint a, lint b){
    return a / gcd(a, b) * b;
}

//桁数を求める
ulint digit(lint n,lint base)
{
    if(n == 0)return 0;
    if(n < 0)n = -n;
    return log10(n)/log10(base) + 1;
}

//桁数を求める
ulint digit(lint n)
{
    if(n == 0)return 0;
    if(n < 0)n = -n;
    return log10(n) + 1;
}

//桁分解, base進数の桁ごとの数値に分解
vector<lint> digitbreak(lint n, lint base)
{ 
    vector<lint> a(base,0);
    while(n >= base)
    {
        a[n % base] ++;
        n /= base;
    }
    a[n] ++;
    return a;
}

//桁分解, 10進数の桁ごとの数値に分解
vector<lint> digitbreak(lint n)
{ 
    return digitbreak(n, 10);
}

Measured memory_bytes: 3684