Jzoj P3912 超氧化钾___数学+推理

题目大意：

给出数$x,y$，求${\sum }_{i=1}^{y}x$ $mod$ $i$。

$1\le x,y\le 1e9$

分析：

$x$ $mod$ $i=$ $x-\lfloor \frac{x}{i}\rfloor \ast i$
那么${\sum }_{i=1}^{y}x$ $mod$ $i=$ ${\sum }_{i=1}^{y}x-\lfloor \frac{x}{i}\rfloor \ast i=$ $xy-{\sum }_{i=1}^y\lfloor \frac{x}{i}\rfloor \ast i$
因为$\lfloor \frac{x}{i}\rfloor$取值至多有$2\sqrt{x}$
所以${\sum }_{i=1}^y\lfloor \frac{x}{i}\rfloor \ast i$中的$\lfloor \frac{x}{i}\rfloor$必定是相同的一段一段的，$\sqrt{n}$的求一下分界点即可，然后对于后面的乘就是等差数列了。

代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>
#include <algorithm>

#define N 100005

using namespace std;

typedef long long ll;

ll Num[N], A, B, Answer = 0;
int cnt;

ll Get_Sum(ll x, ll y)
{
	if (x > y) return 0;
	return (x + y) * (y - x + 1) / 2; 
} 

int main()
{
	scanf("%lld %lld", &A, &B);
	for (ll i = 1; i * i <= A; i++) Num[++cnt] = A / i;       			
	for (ll i = 2; i <= sqrt(A) + 1; i++) 
	    if (A / i != A / (i - 1)) Num[++cnt] = i - 1;
	Num[++cnt] = 0;
	sort(Num + 1, Num + cnt + 1);
	Answer = (ll)A * B; 
	for (int i = 1; i < cnt; i++)
		if (Num[i + 1] > B) { Answer = Answer - (ll)Get_Sum(Num[i] + 1, B) * (A / B); break; }
		               else Answer = Answer - (ll)Get_Sum(Num[i] + 1, Num[i + 1]) * (A / Num[i + 1]); 
	printf("%lld\n", Answer);
	return 0;
}