Codeforces 837E Vasya's Function（数论、gcd）

题目原网页->   Codeforces 837E Vasya's Function

写这道题参考了->https://blog.csdn.net/calabash_boy/article/details/76652792

---

• 这道题不存在多组样例输入。一开始以为会有多组样例输入，打算预处理出1~1e6所有数的因数，非常非常麻烦。

1. 当输入的x和y不互质的时候，需要先计算出gc=gcd(x,y)。经过推理可以发现  f ( a , b ) == f ( a / gc , b / gc )   因为此时是把gc当作一个整体、一个单位看待，相当于1。对结果不会产生任何影响。这个结论应用在后面可以减少计算次数。

2. build_factor函数，用于计算出x的分解质因数后的结果。2^45=35,184,372,088,832。最多不会超过45个因数，开第二维70大小的数组绰绰有余。计算到最后如果x没有被除尽，则把x也放入数组。

3. 当y==0，f(x,y)=0,计算结束。x在之前已经被处理成质因数的乘积的形式，接下来的计算过程中不再需要x,只需要对存放质因数的factor数组进行处理即可。

4. while的每一个循环开始时，y已经与x互质（此时x表现为factor数组里不为 -1 的数的乘积的形式）此时需要不停地对y-1-1-1，直到和x产生不为1的公约数。此时简化过程，要找到最少减几个1可以和x产生公约数。就对y mod(此时x拥有的质因数)，取余数最小的那个。假设余数最小为r , 那么y减去r个1，就会和x 产生公约数。需要操作r次，所以answer+=r.

5. 接着上一步，当前x和y可以产生不为1的公约数了。那么就x/=(当前gc),y/(当前gc)。如果直接计算此时的gcd，则需要重新计算出新的x的factor数组。减少麻烦，可以利用之前的factor数组,对于既能整除x,又能整除y的factor[i]。就使y/=factor[i]。facor[i]=-1,相当于x/=factor[i]。之后不再使用factor[i]。 

---

 第一次会wa,原因：if(y%factor[i]==0)写成了if(temp%factor[i]==0),,,temp随时变化，怎么会accepted。

#include<stdio.h>
#include<math.h>
long long factor[70];
int cnt=0;
void build_factor(long long x)
{
    long long f=sqrt(x)+1;
    for(long long i=2;i<=f;i++)
    {
        while(x%i==0)
        {
            factor[cnt++]=i;
            x/=i;
        }
    }
    if(x>1)
        factor[cnt++]=x;
}

long long cut(long long y)
{
    for(int i=0;i<cnt;i++)
    {
        if(factor[i]==-1)
            continue;
        if(y%factor[i]==0)
        {
            y/=factor[i];
            factor[i]=-1;
        }
    }
    return y;
}

long long gcd(long long a,long long b)
{
    if(b==0)
        return a;
    else
        return gcd(b,a%b);
}

int main()
{
    long long x,y,gc,answer;
    scanf("%lld%lld",&x,&y);
    gc=gcd(x,y);
    x=x/gc;
    y=y/gc;
    //printf("gc=%lld x=%lld y=%lld\n",gc,x,y);
    build_factor(x);
    answer=0;
    while(y>0)
    {
        long long temp=y;
        for(int i=0;i<cnt;i++)
        {
            if(factor[i]==-1)
                continue;
            if(y%factor[i]<temp)
                temp=(y%factor[i]);
        }
        answer+=temp;
        y-=temp;
        y=cut(y);
    }
    printf("%lld\n",answer);
    return 0;
}