Light OJ 1014 - Ifter Party

思路：

先求P-L的所有因子，然后用dfs搜索所有的约数，然后输出大于L的所有约数。

#include <stdio.h>
#include <iostream>
#include <string>
#include <string.h>
#include <algorithm>
#include <stdlib.h>
#include <math.h>
#include <vector>
#include <map>
#include <utility>

using namespace std;
typedef long long ll;
const int maxn = 100005;
int prime[maxn];
int isprime[maxn];


void getprime(){//筛选素数
    memset(isprime, 0,sizeof(isprime));
    int cont = 1;
    for(int i = 2;i < maxn;i++){
        if(!isprime[i])
            prime[cont++] = i;
        if(!isprime[i])
            for(int j = i + i;j < maxn;j += i)
                isprime[j] = 1;
    }
//    for(int i = 0;i < cont;i++)
//        printf("%d ",prime[i]);
}
vector< pair<int, int> >factor;
vector<int>ans;

vector< pair<int, int> >primeFactor(int n){//分解质因子
    vector< pair<int, int> >v;
    //此处需要特别注意，应该强转为long long ，因为prime[i]*prime[i]当i很大的时候非常大。
    for(int i = 1;(ll)prime[i]*prime[i] <= (ll)n;i++){
        int p = 0;
        while(n%prime[i] == 0){
            p++;
            n/=prime[i];
        }
        if(p)
            v.push_back(make_pair(prime[i], p));
    }
    if(n != 1)
        v.push_back(pair<int, int>(n, 1));
    return v;
}
int pow(int x,int y){
    int ans = 1;
    while(y){
        if(y&1) ans*=x;
        y >>= 1;
        x *= x;
    }
    return ans;
}
void dfs(int num,int fac){
    if(num == factor.size()){
        ans.push_back(fac);
        return;
    }
    pair<int, int>pp;
    pp = factor[num];
    for(int i = 0;i <= pp.second;i++){
        dfs(num+1, fac*pow(pp.first, i));
    }
}
int main(){
    getprime();
    int t;
    scanf("%d",&t);
    for(int cas = 1;cas <= t;cas++){
        int p,l;
        scanf("%d%d",&p,&l);
        factor.clear();
        ans.clear();
        factor = primeFactor(p-l);
        dfs(0, 1);
//        for(int i = 0;i < ans.size();i++)
//            printf("%d ",ans[i]);
        sort(ans.begin(), ans.end());
        int it;
        for(it = 0;it < ans.size() && ans[it] <= l; it++);
        printf("Case %d:",cas);
        if(it == ans.size())
            printf(" impossible\n");
        else{
            for(;it != ans.size();it++)
                printf(" %d",ans[it]);
            printf("\n");
        }
        
//        pair<int, int>pp;
//        for(int i = 0;i < factor.size();i++){
//            pp = factor[i];
//            printf("%d %d\n",pp.first,pp.second);
//        }
//        printf("-----------\n");
        
    }
    return 0;
}