jzoj3912. 超氧化钾（B组——Day6）

jzoj3912. 超氧化钾（B组——Day6）

题目

Description

作为一个化学小考只有77分的同学，小W感觉很有压力，所以他决定搞化学。今天他开始研究超氧化钾，这是一种很神奇的离子化合物，在高温下空间结构是立方体，如下图所示。有一天教黄交给了小W一个任务：根据教黄多年淘金的经验，他发现超氧化钾有一个神奇的性质就是如果一组有序的整数(x,y,z)满足z=x mod 1+x mod 2+…+x mod y，那么对于超氧化钾来说这个位置上的离子就是一个关键离子。通过用回旋加速器得到高速粒子来攻击这个离子，他就可以获得铜铝制金法的催化剂。现在教黄想知道对于一个给定的x和y对应的z值是多少。

Input

两个整数x,y。

Output

一个整数z。

Sample Input

1 1

Sample Output

0

Data Constraint

40%的数据，y<=10^5。
100%的数据，1<=x,y<=10^9。

解析

等差数列

找规律，举例：
x=50 ，y=100时，得到x%（1~y）的值为：
0 0 2 2 0 2 1 2 5 0 6 2 11 8 5 2 16 14 12 10 8 6 4 2 0 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50

可以发现当x<y时，余数就是x，可以直接得出这一部分：z=（y-x）*x；
当x>=y时，
x/2+1~x的余数是一个公差为1的等差数列
x/3+1~x/2的余数，是一个公差为2的等差数列
…………………………………………
x/i+1~x/(i-1)的余数，是一个公差为i-1的等差数列

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

long long x,y,z;

int main(){
	scanf("%lld%lld",&x,&y);
	if (x<y) {
		z=(y-x)*x;
		y=x-1;
  	}
	long long i=1;
	while (x/i>y) i++;
	z+=((y-x/i)*(x%(x/i+1)+x%y))/2;
	while(x/i>100000){
		i++;
		long long l=x/i,r=x/(i-1);
		z+=((r-l)*(x%(l+1)+x%r))/2;
	}
	y=x/i;
	while (y>0){
		z+=x%y;
		y--;
	}
	printf("%lld",z);
	return 0;
}