HDU 4710 Balls Rearrangement

这道题太搞笑了，估计是HDU的测试数据太弱了，我直接暴力加一点小技巧都过了，sigh。。。

不过有一点值得注意的是，64位数据最好用cin/cout，否则可能会产生输入和输出错误。

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4710

好吧，就是求a,b 的最小公倍数Q，然后如果Q>n，就直接进行上限为Q的模拟，当然，要加一点小技巧。

技巧：设现在要往A，B箱子里放球，箱子编号分别为la，lb，那么剩余量就是A-la，和B-lb。用剩余量最小的那个直接进行加法操作，相见代码。。。

貌似这个技巧叫做跳跃求值，我不知道。。。

如果Q<n，就进行上限为Q的模拟，模拟n/Q次，再模拟一次，上限为n%Q。

总量就是n/Q*calcu(Q,a,b) + calcu(n%Q,a,b) .

上代码：

#include<stdio.h>
#include<iostream>
using namespace std;
#include<math.h>
#include<string>
#include<map>
#include<queue>
#include<stdlib.h>
#include<algorithm>

#define repA(p,q,i)  for( int (i)=(p); (i)!=(q); ++(i) )
#define repAE(p,q,i)  for( int (i)=(p); (i)<=(q); ++(i) )
#define repD(p,q,i)  for( int (i)=(p); (i)!=(q); --(i) )
#define repDE(p,q,i)  for( int (i)=(p); (i)>=(q); --(i) )
#define ll long long

ll gcd(ll a, ll b);
ll calcu( ll round, ll ca, ll cb );

int main()
{
    ll a,b,n;
    int test;  scanf("%d",&test);
    while(test--)
    {
        cin>>n>>a>>b ;
        ll GCD = gcd(a, b) ;
        ll circle = a*b/GCD ;
        if(a == b)
        {
             printf("0\n");
             continue;
        }
        
        if( circle > n )
        {
            cout<< calcu(n,a,b) << endl ;
            continue;
        }
        ll ans = n/circle * calcu(circle,a,b);
    
        ans += calcu(n%circle,a,b);
        cout<<ans<<endl; 
        
    }
    return 0;
}

ll gcd(ll a, ll b)
{  return b==0?a:gcd(b,a%b);  }


ll calcu( ll round, ll ca, ll cb )
{
     ll la=0, lb=0, gone, ans=0;
     int i=0;
     while( i < round )
     {
         gone = min( ca-la , cb-lb ) ;
         if( i+gone < round )
         {
             i += gone;
             ans += gone * abs(la-lb) ;
             la = (la+gone) % ca ;
             lb = (lb+gone) % cb ; 
         }
         else 
         {
              ans += (round-i) * abs(la-lb) ;
              break;
         }
     }
     return ans ;
} 

当然我以为我的代码应该是野鸡代码，水过去的，不过后来搜了下题解，发现别人的好像跟我差不多思路，代码还没有我简洁。。。。

当然，简洁换来的是速度的下降，本质没什么区别。

贴一下吧，摘自：http://www.cnblogs.com/kuangbin/p/3214791.html

#include <stdio.h>
#include <algorithm>
#include <iostream>
#include <string.h>
#include <set>
#include <map>
#include <vector>
#include <queue>
#include <string>
#include <math.h>
using namespace std;
long long gcd(long long a,long long b)
{
    if(b==0)return a;
    else return gcd(b,a%b);
}
long long lcm(long long a,long long b)
{
    return a/gcd(a,b)*b;
}
long long calc(int n,int a,int b)
{
    long long ans = 0;
    int i = 0;
    int ta=0,tb=0;
    int p = 0;
    while(i < n)
    {
        if(ta+a >= n && tb+b >= n)
        {
            ans += (long long)(n-i)*p;
            i = n;
            continue;
        }
        if(ta+a < tb+b)
        {
            ans += (long long)p*(ta+a-i);
            i = ta+a;
            p = i - tb;
            ta+=a;
        }
        else if(ta+a==tb+b)
        {
            ans+= (long long)p*(ta+a-i);
            i = ta+a;
            ta+=a;
            tb+=b;
            p = 0;
        }
        else
        {
            ans += (long long)p*(tb+b-i);
            i = tb+b;
            tb+= b;
            p = i-ta;
        }
    }
    return ans;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int T;
    int n,a,b;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&a,&b);
        if(a==b)
        {
            printf("0\n");
            continue;
        }
        if(a < b)swap(a,b);
        long long LCM = lcm(a,b);
        if(LCM >= n)
        {
            printf("%I64d\n",calc(n,a,b));
            continue;
        }
        long long tmp = calc(LCM,a,b);
        long long ans = tmp * (n/LCM)+calc(n%LCM,a,b);
        printf("%I64d\n",ans);
    }
    return 0;
}