codeforces 837E 数论

最新推荐文章于 2024-06-29 17:01:15 发布

原创最新推荐文章于 2024-06-29 17:01:15 发布 · 1k 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

数学同时被 3 个专栏收录

13 篇文章

订阅专栏

gcd

8 篇文章

订阅专栏

数论概论

3 篇文章

订阅专栏

本文介绍了一种计算特殊数论函数的方法，该函数定义为F(a, 0)=0，F(a, b)=1+F(a, b-GCD(a, b))。通过分析给出的x和y(<=1e12)，利用质因子处理来高效计算F(x, y)。

Vasya is studying number theory. He has denoted a function f(a, b) such that:

f(a, 0) = 0;
f(a, b) = 1 + f(a, b - gcd(a, b)), where gcd(a, b) is the greatest common divisor of aand b.

Vasya has two numbers x and y, and he wants to calculate f(x, y). He tried to do it by himself, but found out that calculating this function the way he wants to do that might take very long time. So he decided to ask you to implement a program that will calculate this function swiftly.

Input

The first line contains two integer numbers x and y (1 ≤ x, y ≤ 10¹²).

Output

Print f(x, y).

Example

Input

3 5

Output

Input

6 3

Output

题意：定义F（a，0） = 0，F（a，b） = 1 + F（a，b - GCD（a，b）。给定 x 和 y （<=1e12）求F（x，y）。

假设 a=A*gcd b=B*gcd 一轮后 A*gcd,(B-1)*gcd 当(B-k)%T == 0时 gcd变化 T为A的因子，由于因子数量很多，所以我们利用质因子来处理因子的问题，要找到最小的k，其实就是找最小的B%T ，这个想想就懂了，然后找到这个乘到gcd里面的因子，这个因子将会是所有y%T==0 的质因子组成，把这些质因子挑出来，其实就相当于找到了这个因子，并从y中剔除。

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

vector<ll> v;

int main()
{
    ll x,y;
    scanf("%lld%lld",&x,&y);
    ll tt=__gcd(x,y);
    x/=tt,y/=tt;
    for(ll i=2;i*i<=x;i++)
    {
        while(x%i==0)
        {
            x/=i;
            v.push_back(i);
        }
    }
    if(x>1) v.push_back(x);
    ll ans=0;
    while(y)
    {
        ll hh=y;
        for(int i=0;i<v.size();i++)
            hh=min(hh,y%v[i]);
        y-=hh;
        ans+=hh;
        vector<ll> v1;
        for(int i=0;i<v.size();i++)
        {
            if(y%v[i]==0)
            {
                y/=v[i];
            }
            else v1.push_back(v[i]);
        }
        v.swap(v1);
    }
    printf("%lld\n",ans );
}