bzoj2257 [Jsoi2009]瓶子和燃料

来源：互联网发布：php全栈框架编辑：程序博客网时间：2024/04/28 00:47

题目

一道数学题。
结论：首先易知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 longusing 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.

阅读全文

0 0