POJ 2657 Comfort（一元线性同余方程）

Description
有n个点顺时针围成一圈（编号1~n），初始时在1点，每次可以顺时针移动k点，n个点中有一些点有障碍物，求最小的k使得在不经过任何一个有障碍物的点就可以到达z点
Input
第一行三个整数n,z,m分别表示点数，目标点以及障碍物数，之后m个整数表示哪些点有障碍物
Output
求最小的k使得在不经过任何一个有障碍物的点就可以到达z点
Sample Input
9 7 2
2 3
Sample Output
3
Solution
为方便取模，先将所有点编号减一变成0~n-1，从1到z枚举k，对于每个k，如果一元线性同余方程k*x=z(mod n)有解g（此处g为最小正整数解），说明每次跳k步可以到达z点，之后枚举每个障碍物判断是否会经过障碍物，对于某障碍物ob[i]，如果有解且最小正整数解y< g，说明会经过这个障碍物，否则不会经过，以此来判断当前k是否满足条件，第一次遇到满足条件的k一定为最小的k
Code

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define maxn 1111
int n,z,k,m,ob[maxn];
int extend_gcd(int a,int b,int &x,int &y)
{ 
    int d=a;
    if(b!=0)
    {
        d=extend_gcd(b,a%b,y,x);
        y-=(a/b)*x;
    }
    else 
    {
        x=1;
        y=0;
    }
    return d;
}
int linear(int a,int b,int c)//求解一元线性同余方程ax=b(mod c)
{
    int x,y;
    int g=extend_gcd(a,c,x,y);
    if(b%g)//无解 
        return -1;
    x=x*(b/g);
    int mod=c/g;
    x=(x%mod+mod)%mod;
    return x;
}
int main()
{
    while(~scanf("%d%d%d",&n,&z,&m))
    {
        z--;
        for(int i=0;i<m;i++)scanf("%d",&ob[i]),ob[i]--;
        for(k=1;k<=z;k++)
        {
            int g=linear(k,z,n);
            if(g!=-1)
            {
                int flag=1;
                for(int i=0;i<m;i++)
                {
                    int x=linear(k,ob[i],n);
                    if(x!=-1&&x<g)flag=0;
                }
                if(flag)break;
            }
        }
        printf("%d\n",k);
    }
    return 0;
}