【其他】计蒜客 ICPC Arab Collegiate Programming Contest 2014 A题 Game of Peace

题意：
有 2 个容器，初始其中一个容器中有 x 个大理石，另一个有 0 个
克隆机可以将数量较多的容器中的大理石复制后倒入数量较少的容器
当复制第 n 次时，数量较多的容器大理石数量加 y，
复制 m 次结束，求出结束时两容器大理石数目的最大公约数

刚开始按照纯模拟来做的，结果超时了
后来看了题解才知道，有可以简化运算的方法。

• 由更相减损术可以知道，gcd(a, b) == gcd(a+b, b)，所以当加上 y 之后就不用算了，此时的最大公约数与结果一样，后续的操作并不会影响最大公约数。
• 由于操作数可以达到 10w 次，加数会很大，注意数据类型。

模拟代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

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

int main()
{
    int T; scanf("%d", &T);
    for(int p=1; p<=T; p++){
        int x, n, y, m; scanf("%d%d%d%d", &x, &n, &y, &m);
        ll a=x, b=0, k=0;
        for(;;){
            if(k==n) { if(a<b) swap(a, b); a+=y; break;}
            // if(k==m) break;  只需要在模拟的基础上去掉这一行就对了
            k++;
            if(a<b) swap(a, b);  // a>=b
            b+=a;
        }

        cout<<"Case "<<p<<": "<<gcd(a, b)<<endl;
    }

    //system("pause");

    return 0;
}

递推思路：

• 令 f(1) = 0, f(2) = 1，这是初始状态。递推方程为 f(i) = f(i-1) + f(i-2) .
• 此时第 1 次为 f(2), f(3)，第 2 次为 f(3), f(4)，…，第 n 次为 f(n+1), f(n+2)；由递推方程可知，f(n+2) > f(n+1)，所以 f(n+2) + y .
• 综上，gcd( f(n+1), f(n+2)+y ) 就是最后要求的最大公约数。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

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

ll f[100];

int main()
{
    int T; scanf("%d", &T);
    for(int p=1; p<=T; p++){
        int x, n, y, m; scanf("%d%d%d%d", &x, &n, &y, &m);
        f[1]=0; f[2]=x;
        for(int i=3; i<=n+2; i++)
            f[i]=f[i-1]+f[i-2];
        printf("Case %d: %lld\n", p, gcd(f[n+1], f[n+2]+y));
    }

    //system("pause");

    return 0;
}