数字对 C组模拟赛

题目大意：

将(1,1)变为(n,i)(i>0)变化形式
1.(a+b,a)
2.(a,a+b)

---

解题思路：

30分

对于30%的数据， 1 <= n <= 1000
f[i][j]表示从（1, 1）变为（i，j）所需最少步数
使用n$^2$dp轻松解决
用 f[i][j]表示a=i，b=j时所使用的最小代价
这样空间时间都会炸
所以只能拿30分

---

60分

对于60%的数据，1 <= n <= 20000
反着来推
给定（a,b），将其还原
唯一还原方法
若 a > b，则（a,b）-> (a-b, b)
若 b >= a，则（a,b) -> (a, b-a)
就是还原(n,i)为(1,1)，这时枚举i
(没打60分代码)
时间复杂度
O（n*logn）

---

100分

对于100%的数据，1 <= n <= 10$^6$
对于给定的n，枚举所有的i，模拟还原（n,i）
取其中最少步数为解
如何加速模拟过程保证复杂度？
类比求gcd的辗转相除法
若 a > b，则（a,b）-> (a%b, b)
若 b >= a，则（a,b) -> (a, b%a)

---

源程序：

30分（dp）：

#include<cstdio>
#include<cstring>
using namespace std;
int f[6001][6001],n,ans,k;
long long min(int x,int y){return x>y?y:x;}
int main()
{
    scanf("%d",&n);
    memset(f,127/3,sizeof(f));
    f[1][1]=0;
    f[1][2]=1;
    f[2][1]=1;
    ans=2147483647;
    for (int i=1;i<=n;i++)
     for (int j=1;j<=n;j++)
      f[i+j][j]=f[i][i+j]=min(f[i][j]+1,min(f[i+j][j],f[i][i+j]));
    for (int i=1;i<=n;i++)
     ans=min(ans,min(f[i][n],f[n][i]));
    printf("%d",ans);
}

AC（GCD优化）：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,ans;
int min(int x,int y){return x>y?y:x;}
int Ngcd(int a,int b)
{
    if (b==1) return a-1;
    if (a==b||a==0||b==0) return 1e9;
    if (a<b) return Ngcd(b%a,a)+b/a;
    else return Ngcd(a%b,b)+a/b;
}
int main()
{
    scanf("%d",&n); 
    ans=2147483647;
    for (int i=1;i<=n;i++)
     ans=min(ans,Ngcd(n,i));
    printf("%d",ans);
}