对m进行分解,然后从 ( n − 1 0 ) 计算到 ( n − 1 n − 1 ) ,对每个数能否被m整除。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49
| #include<cstdio> #include<cstring> #include<cmath> using namespace std; const int maxn=100010; int p[110],e[110],pnum; int ans[maxn],pos; int main(){ int n,m; while(~scanf("%d%d",&n,&m)){ pnum=pos=0; int k=sqrt(m); for(int i=2;i<=k;++i){ if(m%i==0){ p[pnum]=i,e[pnum]=0; while(m%i==0) m/=i,++e[pnum]; ++pnum; } if(m==1) break; } if(m>1){ p[pnum]=m; e[pnum]=1; ++pnum; } int cnt=pnum; for(int i=1;i<n;++i){ int a=n-i,b=i; for(int j=0;j<pnum;++j){ while(a%p[j]==0){ a/=p[j]; --e[j]; if(e[j]==0) --cnt; } while(b%p[j]==0){ b/=p[j]; ++e[j]; if(e[j]==1) ++cnt; } } if(cnt==0) ans[pos++]=i+1; } printf("%d\n",pos); if(pos==0) printf("\n"); for(int i=0;i<pos;++i) printf("%d%c",ans[i],i==pos-1?'\n':' '); } return 0; }
|
** 本文迁移自我的CSDN博客,格式可能有所偏差。 **