1225 余数之和

找规律。。。边界要特殊考虑。。。

#include<bits/stdc++.h>
using namespace std;

const long long mod=1e9+7;

long long multimod(long long x,long long y,long long mod)
{
	long long ret=0;
	while(y)
	{
		if(y&1)
			ret=(ret+x)%mod;
		x=(x<<1)%mod;
		y>>=1;
	}
	return ret;
}

int main()
{
	long long n,ans,sq,i,inv2=500000004;
	while(~scanf("%lld",&n))
	{
		sq=sqrt(n);
		ans=multimod(n,n,mod);
		for(i=1;i<=sq;i++)
		{
			if(n/i==sq)
			{
				ans=((ans-sq*sq)%mod+mod)%mod;
				continue;
			}
			ans=((ans-multimod(n/i,i,mod))%mod+mod)%mod;
			ans=((ans-multimod(multimod(i,multimod((n/i+n/(i+1)+1)%mod,(n/i-n/(i+1))%mod,mod),mod),inv2,mod))%mod+mod)%mod;
		}
		printf("%lld\n",ans);
	}
}