bzoj2257 [Jsoi2009]瓶子和燃料

题目

一道数学题。
结论：首先易知k个瓶子能导出的油最小值一定是k个瓶子容量的最大公因数 于是问题转化成了在n个数中选择k个使最大公因数最大。
所以把每个数分解因数，选出公因数个数大于等于k中最大的，即可。

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cmath>
#include<map>
#define MAXV 1000000000
#define LL long long
using namespace std;
bool P[100001];
int prime[100001],size;
map <LL,int>E;
map <LL,int>::iterator id;
LL ans;
LL n,k,x;
void read(LL &x)
{
    char c;
    x=0;
    c=getchar();
    if(c==-1)return;
    while(c>'9'||c<'0')c=getchar();
    while(c<='9'&&c>='0')x=x*10+c-'0',c=getchar();
}
LL ksm(LL a,LL b,LL c)
{
    LL ans=1;
    a=a%c;
    while(b)
    {
        if(b&1)ans=ans*a%c;
        a=a*a%c;
        b=b>>1;
    }
    return ans;
}
bool is_P(LL x)
{
    if(x==1)return false;
    if(x==2)return true;
    if(x==5)return true;
    if(x==13)return true;
    if(ksm(2,x-1,x)!=1)return false;
    if(ksm(5,x-1,x)!=1)return false;
    //if(ksm(13,x-1,x)!=1)return false;
    return true;
}
void divide(LL x)
{
    if(is_P(x)){E[x]++;return;}
    for(int i=1;i<=sqrt(x);i++)
    {
        if(x%i==0)E[i]++;
        if(x%i==0&&i!=x/i)E[x/i]++;
    }
}
int main()
{
    memset(P,true,sizeof(P));
    P[1]=false;
    for(int i=2;i<=100000;i++)
    {
        if(P[i])
        {
            if(i<=100000/i)
            {
                for(int j=i*i;j<=100000;j+=i)
                    P[j]=false;
            }
        }
    }
    for(int i=1;i<=100000;i++)
        if(P[i])prime[++size]=i;
    read(n),read(k);
    for(int i=1;i<=n;i++)
    {
        read(x);
        divide(x);
    }
    ans=1;
    for(id=E.begin();id!=E.end();id++)
    {
        //cout<<id->first<<" "<<id->second<<endl;
        if(id->second>=k)ans=max(ans,id->first);
    }
    cout<<ans;
    return 0;
} 


代码写的丑233.