%%%

%%%

UVA 1635 Irrelevant Elements (唯一分解定理 + 组合数)

这题需要分解的m为1e9不能先打表后分解
直接分解外层循环为sqrt(m)

然后再利用组合数的递推式进行求解

#include <iostream>
#include <vector>
#include <cstring>
#include <cmath>

using namespace std;
const int maxn = 100000;

vector<int> primes;
bool flag[maxn];

void factors(int n) {
    int m = floor(sqrt(n) + 0.5);
    for(int i = 2; i <= m; ++i) {
        if(n % i == 0) {
            primes.push_back(i);
            while(n % i == 0) n /= i;
        }
    }
    if(n > 1) primes.push_back(n);
}

int main() {
    int n, m;
    while(cin >> n >> m) {
        memset(flag, 0, sizeof(flag));
        primes.clear();
        factors(m);
        n--;
        for(int i = 0; i < primes.size(); ++i) {
            int p = primes[i], e = 0;
            int min_e = 0, x = m;
            while(x % p == 0) { x /= p; min_e++; }
            for(int j = 1; j < n; ++j) {
                x = n - j + 1;
                while(x % p == 0) x /= p, e++;
                x = j;
                while(x % p == 0) x /= p, e--;
                if(e < min_e) flag[j] = true;
            }
        }
        vector<int> ans;
        for(int i = 1; i < n; ++i)
            if(!flag[i]) ans.push_back(i + 1);
        cout << ans.size() << endl;
        if(ans.size()) cout << ans.front();
        for(int i = 1; i < ans.size(); ++i)
            cout << ' ' << ans[i];
        cout << endl;
    }
    return 0;
}

阅读更多
上一篇Codeforces 954C - Matrix Walk(思维 + 模拟)
下一篇UVA 1262 Password (编码与解码)
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭