UVA1635-唯一分解定理的基本应用2

原题：https://vjudge.net/problem/UVA-1635

这是一个极其典型的“从素因子角度出发”的题目，下面是我的代码：

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 using namespace std;

 const int maxn = 1e5 + ;
 int md[],mz[];
 int ans[maxn];
 int n,m,lenth;
 int c[];

 bool check(int n,int i)
 {
     int x=n-i;
     int y=i;
     for(int k=;k<lenth;k++)
     {
         int p=md[k];
         while(x%p==)
         {
             x/=p;
             c[k]++;
         }
         while(y%p==)
         {
             y/=p;
             c[k]--;
         }//上下两个循环不能合并,因为这是一个递增积累过程,
     }//也正因如此才需要额外开辟一个数组c.
     for(int k=;k<lenth;k++)
         if(c[k]<mz[k])return false;
     return true;
 }

 int factor()
 {//用来分解m
     memset(md,,sizeof(md));
     memset(mz,,sizeof(mz));
     int c=m,cur=;
     int p=sqrt(c);
     for(int i=;i<=p;i++)
     {
         bool ok=false;
         while(c%i==)
         {
             md[cur]=i;
             mz[cur]++;
             c/=i;
             ok=true;
         }
         if(ok)cur++;
         if(c==)break;
     }
     if(c>)md[cur]=c,mz[cur]=,cur++;
     return cur;
 }

 int main()
 {
     //freopen("input.txt","r",stdin);
     //freopen("text.txt","w",stdout);
     while(cin>>n>>m)
     {
         lenth=factor();//一开始忘记了记录lenth
         int cnt=;
         memset(c,,sizeof(c));
         //
         for(int i=;i<n-;i++)
             if(check(n,i))ans[cnt++]=i+;
         printf("%d\n",cnt);
         for(int i=;i<cnt;i++)
         {
             if(i>)printf(" ");
             printf("%d",ans[i]);
         }
         printf("\n");
     }
     return ;
 }