题意:求sum{abs(i%a-i%b)|i=0..n-1},1<=n<=1000000000,1<=a,b<=100000
这题不难,但是我不会。
看了网上的题解,居然是根据lcm暴力模拟。
从0..n-1模拟,找到下一个最近%a=0或%b=0的断点,计算这一段的变化值。
由于这一段差值是一样的,可以一次计算一段的总差值。
为什么暴力可行呢。似乎没有人给出合理解释。
我来研究一下。
n直接取最大值1000000000。
那么当lcm(a,b)比较大时候,说明min(a,b)至少大于10000。
那么这样暴力模拟跳的断点最多只会有O(n/a+n/b),不会大于10^6。
当lcm比较小的时候,利用循环的特点,每lcm个数是一样的。
这样只要暴力跳lcm个数。
此时lcm最大就是a*b,断点最多是O(lcm/a+lcm/b),即O(b+a)。
所以复杂度是可接受的。
这题不难,但是我不会。
看了网上的题解,居然是根据lcm暴力模拟。
从0..n-1模拟,找到下一个最近%a=0或%b=0的断点,计算这一段的变化值。
由于这一段差值是一样的,可以一次计算一段的总差值。
为什么暴力可行呢。似乎没有人给出合理解释。
我来研究一下。
n直接取最大值1000000000。
那么当lcm(a,b)比较大时候,说明min(a,b)至少大于10000。
那么这样暴力模拟跳的断点最多只会有O(n/a+n/b),不会大于10^6。
当lcm比较小的时候,利用循环的特点,每lcm个数是一样的。
这样只要暴力跳lcm个数。
此时lcm最大就是a*b,断点最多是O(lcm/a+lcm/b),即O(b+a)。
所以复杂度是可接受的。
其实这是一种玩数据范围的题。。仔细分析数据之间的关系就好办啦
注意第二种的计算中,(n/lcm) 要先算,否则会爆long long。因此wa两发。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
long long gcd(long long a,long long b){
return b==0?a:gcd(b,a%b);
}
long long work(long long n,long long a,long long b){
long long au=0,bu=0;
long long ret=0,cnt=0,len;
while(cnt<n){
len=min(a-au,b-bu);
if (cnt+len>n) len=n-cnt;
ret+=abs(au-bu)*(len);
au=au+len;
bu=bu+len;
cnt+=len;
if (au>=a) au-=a;
if (bu>=b) bu-=b;
}
return ret;
}
int cas;
long long n,a,b,lcm,ans;
int Main(){
//freopen("4611in.txt","r",stdin);
scanf("%d",&cas);
while(cas--){
scanf("%I64d%I64d%I64d",&n,&a,&b);
if (a>b) swap(a,b);
lcm=a/gcd(a,b)*b;
if (lcm>n){
ans=work(n,a,b);
}
else {
ans=work(lcm,a,b)*(n/lcm)+work(n%lcm,a,b);
}
printf("%I64d\n",ans);
}
return 0;
}