[luogu4571 JSOI2009] 瓶子和燃料 (数论)

传送门

Solution

题目说的很迷，但可以发现两个瓶子互相倒最少是容积的gcd
那么题目就转化为求其中选k个瓶子gcd的最大值，这个可以分解因数，枚举因数得到

Code

//By Menteur_Hxy
#include <map>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define Re register
#define Fo(i,a,b) for(Re int i=(a),_=(b);i<=_;i++)
#define Ro(i,a,b) for(Re int i=(b),_=(a);i>=_;i--)
//#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin)),p1==p2?EOF:*p1++)
using namespace std;
typedef long long LL;

char buf[1<<21],*p1,*p2;
inline int read() {
	int x=0,f=1;char c=getchar();
	while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
	while(isdigit(c)) x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}

const int MAXN=1010,MAX=1e7;
int n,k,tot,ans;
int vis[MAX],a[MAX];
map<int,int> id;

int main() {
	n=read(),k=read();
	Fo(i,1,n) {
		int x=read(),sqr=sqrt(x);
		Fo(i,1,sqr) if(x%i==0) {
			if(!id[i]) id[i]=++tot,a[tot]=i;vis[id[i]]++;
			if(!id[x/i]) id[x/i]=++tot,a[tot]=x/i;vis[id[x/i]]++;
		}
	}
	Fo(i,1,tot) if(vis[i]>=k) ans=max(ans,a[i]);
	printf("%d",ans);
	return 0;
}