BZOJ-2257 瓶子和燃料 分解因数+数论方面乱搞（裴蜀定理）

一开始真没想出解法。。。后来发现那么水。。。。

2257: [Jsoi2009]瓶子和燃料
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 970 Solved: 577
[Submit][Status][Discuss]

Description
jyy就一直想着尽快回地球，可惜他飞船的燃料不够了。
有一天他又去向火星人要燃料，这次火星人答应了，要jyy用飞船上的瓶子来换。jyy
的飞船上共有 N个瓶子(1<=N<=1000) ，经过协商，火星人只要其中的K 个 。 jyy
将 K个瓶子交给火星人之后，火星人用它们装一些燃料给 jyy。所有的瓶子都没有刻度，只
在瓶口标注了容量，第i个瓶子的容量为Vi（Vi 为整数，并且满足1<=Vi<=1000000000 ） 。
火星人比较吝啬，他们并不会把所有的瓶子都装满燃料。他们拿到瓶子后，会跑到燃料
库里鼓捣一通，弄出一小点燃料来交差。jyy当然知道他们会来这一手，于是事先了解了火
星人鼓捣的具体内容。火星人在燃料库里只会做如下的3种操作：1、将某个瓶子装满燃料；
2、将某个瓶子中的燃料全部倒回燃料库；3、将燃料从瓶子a倒向瓶子b，直到瓶子b满
或者瓶子a空。燃料倾倒过程中的损耗可以忽略。火星人拿出的燃料，当然是这些操作能
得到的最小正体积。
jyy知道，对于不同的瓶子组合，火星人可能会被迫给出不同体积的燃料。jyy希望找
到最优的瓶子组合，使得火星人给出尽量多的燃料。

Input
第1行：2个整数N,K，
第2..N 行：每行1个整数，第i+1 行的整数为Vi

Output
仅1行，一个整数，表示火星人给出燃料的最大值。

Sample Input
3 2
3
4
4

Sample Output
4

HINT
选择第2 个瓶子和第 个瓶子，火星人被迫会给出4 体积的容量。

Source

首先是裴蜀定理：
若a,b是整数,且（a,b)=d，那么对于任意的整数x,y,ax+by都一定是d的倍数，特别地，一定存在整数x,y，使ax+by=d成立。
所以结果一定是gcd。那么分解因数，扫描一遍找到最大的，且数量>=k（m）的因数，即为答案

代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}

int n,m;
int yz[10000010];
int cnt=0,num;

void work(int x)
{
    for (int i=1; i<=sqrt(x); i++)
        if (x%i==0) {yz[++cnt]=i;if (i!=x/i) yz[++cnt]=x/i;} 
}

int main()
{
    n=read(),m=read();
    for (int i=1; i<=n; i++) {int x=read();work(x);}
    sort(yz+1,yz+cnt+1);
    num=1;
    for (int i=cnt; i>=1; i--)
        {
            if (yz[i]==yz[i+1]) {num++;if (num==m) {printf("%d\n",yz[i+1]);break;}}
                else  num=1;
        }
    return 0;   
}