Jzoj3192 球
2017-09-25 15:41
134 查看
小包有N个球和A个盒子。球标号为0到N-1,盒子编号为0到A-1。标号为x的球放在下标为x mod A的盒子里。
然后他得到了B个新盒子,编号为0到B-1。他想把所有球从旧盒子里面拿出来放到新盒子,标号为x的球放在下标为x mod B的盒子里。
对于某个球来说,如果原本放在下标为a的旧盒子,现在放在下表为b的新盒子,那么费用为|a-b|。
求出总花费。a,b<10^6 N<10^9
这个题比较有趣
首先我们发现一个性质,对于0~a-1这a个数,我们只需要知道0对应的是哪个数,就可以快速算出这a个数的贡献
比如a=3 b=5
012|012|012|012|012
012|340|123|401|234
对于第一组,我们知道贡献显然为0
第二组贡献为3*2+2*1=8
第三组类似
另外呢,我们发现对于c=lcm(a,b)来说,0~c-1和c~2c-1的答案相同,即有循环
所以我们可以将时间复杂度优化到O(a)
#include<stdio.h>
#include<algorithm>
#define L long long
using namespace std;
L gcd(L a,L b){
for(L c;b;a=b,b=c) c=a%b;
return a;
}
inline L f(L i,L j,L k){
L S1=min(i,j-k)*k;
L S2=(i-min(i,j-k))*(j-k);
return S1+S2;
}
int main(){
L i,n,a,b,S=0,r,c,S1=0;
scanf("%lld%lld%lld",&n,&a,&b);
r=gcd(a,b); c=a/r*b;
if(a>b) swap(a,b);
for(i=0;i<c;i+=a)
S1+=f(a,b,i%b);
S+=S1*(n/c); n%=c;
for(i=0;i<n;i+=a)
S+=f(a,b,i%b);
if(i>a){ i-=a; S-=f(a,b,i%b); }
for(;i<n;++i) S+=abs((i%a)-(i%b));
printf("%lld\n",S);
}
然后他得到了B个新盒子,编号为0到B-1。他想把所有球从旧盒子里面拿出来放到新盒子,标号为x的球放在下标为x mod B的盒子里。
对于某个球来说,如果原本放在下标为a的旧盒子,现在放在下表为b的新盒子,那么费用为|a-b|。
求出总花费。a,b<10^6 N<10^9
这个题比较有趣
首先我们发现一个性质,对于0~a-1这a个数,我们只需要知道0对应的是哪个数,就可以快速算出这a个数的贡献
比如a=3 b=5
012|012|012|012|012
012|340|123|401|234
对于第一组,我们知道贡献显然为0
第二组贡献为3*2+2*1=8
第三组类似
另外呢,我们发现对于c=lcm(a,b)来说,0~c-1和c~2c-1的答案相同,即有循环
所以我们可以将时间复杂度优化到O(a)
#include<stdio.h>
#include<algorithm>
#define L long long
using namespace std;
L gcd(L a,L b){
for(L c;b;a=b,b=c) c=a%b;
return a;
}
inline L f(L i,L j,L k){
L S1=min(i,j-k)*k;
L S2=(i-min(i,j-k))*(j-k);
return S1+S2;
}
int main(){
L i,n,a,b,S=0,r,c,S1=0;
scanf("%lld%lld%lld",&n,&a,&b);
r=gcd(a,b); c=a/r*b;
if(a>b) swap(a,b);
for(i=0;i<c;i+=a)
S1+=f(a,b,i%b);
S+=S1*(n/c); n%=c;
for(i=0;i<n;i+=a)
S+=f(a,b,i%b);
if(i>a){ i-=a; S-=f(a,b,i%b); }
for(;i<n;++i) S+=abs((i%a)-(i%b));
printf("%lld\n",S);
}
相关文章推荐
- jzoj3192 【中山市选2013】球
- jzoj 3192.【中山市选2013】球
- JZOJ NOIP2014模拟 8.13
- 【JZOJ】3187 的士
- JZOJ 1768. [NOI2001]炮兵阵地
- JZOJ 1756.得分
- JZOJ3975. 【NOI2015模拟1.17】串
- JZOJ4374. 【GDOI2016模拟3.9】暴走的图灵机
- JZOJ 4466【GDOI2016模拟4.22】无界单词
- {题解}[jzoj4489] GDOI2016 Day2_T1 SigemaGo
- [JZOJ 3395] Freda的传呼机
- {题解}[jzoj1299] 洗盘子(cleanup.pas/cpp)
- {题解}[jzoj1298]牛棚(graze2.pas/c/cpp)
- 【JZOJ 3463】 军训 training
- 【2011集训队出题】【BZOJ 2154】【JZOJ 1938】Crash的数字表格
- JZOJ 3453. 连通块(connect)
- JZOJ. 1217. 软件开发
- JZOJ1241. Number
- JZOJ1421. 二叉树
- jzoj 1391. 【2012.02.25普及组】危险系数