CH Round #53 -【Nescafé 32】杯NOIP模拟赛 GCD Path 题解

GCD Path CH Round #53 -【Nescafé 32】杯NOIP模拟赛

描述

给定一张N个点的有向图，点i到点j有一条长度为 i/(gcd(i,j))的边。有Q个询问，每个询问包含两个数x和y，求x到y的最短距离。

输入格式

第一行包含两个用空格隔开的整数，N和Q。
接下来Q行，每行两个数x和y。

输出格式

输出Q行整数，表示从x到y的最短距离。

样例输入

6 2
4 6
2 5

样例输出

2
2

数据范围与约定

• 对于30%的数据，1<=N<=100。
• 对于70%的数据，1<=N<=10^5。
• 对于100%的数据，1<=N<=10^7，1<=x,y<=N，Q<=10^5。

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++

这道题刚开始看用Floyd写了一遍，可是只过30……
然后看别人的题解，发现原来用数学的知识…
假设求x->y的最短路，则直接走 这条路  长度为x/gcd（x，y）设这个长度的质因数分解为a1*a2*a3*a4……(两项可以相等）
然后要用到一个结论：
若a>=2,b>=2,则a+b<=a*b
移项就是 （1-a）*（1-b）>=1  这是显然的。
所以我们不妨把这个长度分开来走，
因为每拆一项都会使答案减小或不边，那我们不妨直接将该数全部分解为质数，一个一个质数来走。   
举个例子
100-1，则100/gcd（100，1）的质因数分解为2*2*5*5
我们不妨使每次走的长度为2 2 5 5，而2+2+5+5=14<100 这样使长度之和达到最小。
所以我们可以这样走 100->50->25->5->1->1  
所以就有了下面的代码……可是发现会有错误答案……
经过仔细的研究我发现有一种情况：自己到自己……按题目的要求应该是1，可是很明显答案是0……
然后我又加了一行if，最后过了……
哎……题目有歧义啊！！

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
int n,q,x,y;


int is_prime(int x)
{
for(int i=2;i<=(int)sqrt(x);i++)
if (x%i==0)return 0;
return 1; 
}


int pr(int x)  
{
int sum=0;
if (is_prime(x))return x;
    int a=2; 
    while(a*a<=x)
    {
        while(x%a==0)
        {
            sum+=a;
            x=x/a;
        }
        a++;
        
    }
    if(x>1)sum+=x; 
    return sum;
} 


int i_gcd(int x,int y)
{
int ii=x;
if (!x||!y)return x>y?x:y;
for(int i;i=x%y;x=y,y=i);
return pr(ii/y); 
}


int main()
{
scanf("%d%d",&n,&q);
for (int i=1;i<=q;i++) 
{ 
scanf("%d%d",&x,&y);
if(x==y)printf("0\n");
else printf("%d\n",i_gcd(x,y));
}
return 0; 
}