#include <bits/stdc++.h>
using namespace std;
const int maxn = 1E5 + 10;
bool prim[maxn];
vector<pair<int, int> > Fta;
vector<int>primes, ans;
int n, m, kase, Irrelevant[maxn];
void init()
{
for (int i = 2; i <= maxn; i++)
if (!prim[i])
{
primes.push_back(i);
for (int j = i + i; j <= maxn; j++)
prim[j] = 1;
}
}
void prime_factors(int n)
{
int m = floor(sqrt(n) + 0.5);
for (int i = 2; i <= m; i++)
{
int cnt = 0, p = i;
if (n % p == 0)
{
while (n % p == 0) n /= p, cnt++;
Fta.push_back(make_pair(p, cnt));
}
}
if (n > 1) Fta.push_back(make_pair(n, 1));
}
int main(int argc, char const *argv[])
{
init();
while (cin >> n >> m)
{
n--;//计算的是C(0~n-1,n-1)!;
ans.clear();
Fta.clear();
prime_factors(m);
memset(Irrelevant, 0, sizeof(Irrelevant));
for (int i = 0; i < Fta.size(); i++)
{
int p = Fta[i].first, cnt = 0, x;
for (int k = 1; k < n; k++)
{
x = n - k + 1;
while (x % p == 0) {x /= p; cnt++;}
x = k;
while (x % p == 0) {x /= p; cnt--;}
if (cnt < Fta[i].second) Irrelevant[k] = 1;
}
}
for (int i = 1; i < n; i++) if (!Irrelevant[i]) ans.push_back(i + 1);
cout << ans.size() << endl;
for (int i = 0; i < ans.size(); i++)
cout << ans[i] << (i == ans.size() - 1 ? "\n" : " ");
if (!ans.size()) cout << endl;//千万别忘了没有答案,也要换行!
}
return 0;
}
唯一分解定理的模板题,顺带附带线性打表求素数
虽然是模板题,但是解法不好想,可以利用公式C(k,n) = (n-k+1)/k*C(n,k-1) 来计算
但还是绕不过去大数(显然上不了大树) 于是就有了很有意思的观点,对于每一个m的质数考察,考察对于每一个c(),是否能被整除,对于下一个c()由公式可得,这样上一个计算出来的指数可以被下一个继承,而本来需要高精度表示的数字,就只需要指数来表示就可以了。十分精妙