HDU 4944 FSF’s game 规律题

据说推出了公式却没有写出来，sigh。。。

== 表示我不是很知情的补了一下。。。

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
#define mod 4294967296
#define N 500010
typedef long long ll;
ll dp[N], gc[N];
void go(){
    memset(gc, 0, sizeof gc);
    for(ll i = 1; i < N; i++)
    {
        for(ll j = i; j < N; j+=i)
        {
            gc[j] += (j/i+1)*(j/i)/2;
            gc[j] %= mod;
        }
    }
}

int main() {
    go();
    dp[1] = 1;
    for(ll i = 2; i < N; i++)dp[i] = dp[i-1] + gc[i]*i, dp[i] %= mod;
    int T, Cas = 1, n; scanf("%d", &T);
    while(T-- > 0) {
        scanf("%d",&n);
        printf("Case #%d: ",Cas++);
        cout<<dp[n]<<endl;
    }
    return 0;
}