题意:见体面,就是求f(x,y)
题解:假如我们模拟递推肯定炸,因为x和y有1e12这么大,因此我们要看看这里边是否存在着无用功,首先因为每次b都会减掉一个gcd(a,b),我们可以将a和b分别写成x1*gcd(a,b)和x2*gcd(a,b),当然这里x1和x2一定是互质的(废话,都提取最大公约数了),到这里你可能恍然大悟,你也可能误入错误的思想无法自拔,因为你一定会觉得肯定要减掉min(x1,x2)次呀,直接除掉不就好了?...... 当然不是了,因为假如你减了m次,此时x1还是x1,而x2是x2-m,但是假如gcd(x1,x2-m)>1呢不就炸了,因为此时减的gcd不再是之前的gcd了,想到这里,这道题就已经可以过了,对于刚才说的情况,我们可以通过暴力a的所有质因数,呢每次该gcd对答案的贡献一定是min(x2,b/(x1的某个质因数));然后我们重复过程直到b==0停止即可。
#include<stdio.h>
#include<algorithm>
using namespace std;
#define ll long long
ll x,y;
ll gcd(ll a,ll b)
{
if(a%b==0)
return b;
return gcd(b,a%b);
}
ll work(ll a,ll b)
{
if(b==0) return 0;
ll g=gcd(a,b);a=a/g;b=b/g;
ll aa=a,bb=b; ll mn=b;
for(ll i=2;i*i<=a;i++)
{
if(a%i) continue;
mn=min(mn,b%i);
while(a%i==0) a/=i;
}
if(a!=1) mn=min(mn,b%a);
return work(aa,bb-mn)+mn;
}
int main(void)
{
scanf("%lld%lld",&x,&y);
printf("%lld\n",work(x,y));
return 0;
}