codeforces 490D D. Chocolate(数论)

题目链接：

codeforces 490D

---

题目大意：

给出两块长方形的巧克力，每次可以选择其中一块吃掉三分之一，或者吃掉一半，问最少多少次可以可以得到两块同样面积的巧克力。

---

题目分析：

• 首先求取两块巧克力面积的gcd,如果$\frac{area_1}{gcd} 和\frac{area_2}{gcd}$的质因数只有2和3的情况下才能有解
• 构建一组解的方法，就是记录两块巧克力质因数2和3的个数，首先去掉多余的3，这样会导致2更多，然后再去掉多余的2就得到了最终结果。
  

  ---

AC代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>

using namespace std;

typedef long long LL;
LL a1,b1,a2,b2,s1,s2;

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

int main ( )
{
    while ( ~scanf ( "%I64d%I64d%I64d%I64d" , &a1 , &b1 , &a2 , &b2 ) )
    {
        s1 = a1*b1;
        s2 = a2*b2;
        LL d = gcd ( s1 , s2 );
        s1 /= d;
        s2 /= d;
        LL temp = s1;
        LL p1,p2,q1,q2;
        p1 = p2 = q1 = q2 = 0;
        bool flag = true;
        while ( temp%2LL == 0LL )
            temp /= 2LL,p1++;
        while ( temp%3LL == 0LL )
            temp /= 3LL,p2++;
        if ( temp != 1LL ) flag = false;
        temp = s2;
        while ( temp%2LL == 0LL )
            temp /= 2LL,q1++;
        while ( temp%3LL == 0LL )
            temp /= 3LL,q2++;
        if ( temp != 1LL ) flag = false;
        if ( flag )
        {
            LL ans = 0;
            while ( p2 != q2 )
            {
                ans++;
                if ( p2 > q2 )
                {
                    p2--;p1++;
                    if ( a1%3LL == 0LL )
                    {
                        a1 /= 3LL;
                        a1 *= 2LL;
                    }
                    else
                    {
                        b1 /= 3LL;
                        b1 *= 2LL;
                    }
                }
                else if ( p2 < q2 )
                {
                    q2--;q1++;
                    if ( a2%3LL == 0LL )
                    {
                        a2 /= 3LL;
                        a2 *= 2LL;
                    }
                    else
                    {
                        b2 /= 3LL;
                        b2 *= 2LL;
                    }
                }
            }
            while ( p1 != q1 )
            {
                ans++;
                if ( p1 > q1 )
                {
                    p1--;
                    if ( a1%2LL == 0LL ) a1/= 2LL;
                    else b1 /= 2LL;
                }
                else if ( p1 < q1 )
                {
                    q1--;
                    if ( a2%2LL == 0LL ) a2 /= 2LL;
                    else b2 /= 2LL;
                }
            }
            printf ( "%I64d\n" , ans );
            printf ( "%I64d %I64d\n" , a1 , b1 );
            printf ( "%I64d %I64d\n" , a2 , b2 );
        }
        else puts ( "-1" );
    }
}