题目链接:http://hihocoder.com/problemset/problem/1165
题面:
幽香今天心情不错,正在和花田里的虫子玩一个益智游戏。
这个游戏是这样的,对于一个数组A,幽香从A中选择一个数a,虫子从A中选择一个数b。a和b可以相同。她们的分数是a*b的因子的个数。
幽香和虫子当然想要获得尽可能的高的分数,你能告诉她们应该选择哪两个数吗。
由于幽香是个非常随意的人,数组A中的元素都是她随机选择的。
一行一个数n,表示A中整数的数量。
接下来一行n个数,分别表示a1,a2,...,an,为A中的元素。
n <= 100000, 1 <= ai <= 100000
保证所有的ai都是随机生成的。
一行表示最大的分数。
2 3 4
6
解题:白天刚好研究了下质数分解的问题,刚好拿代码来套一下正好。第一次拿T-shirt,刚刚好卡在最后,昨天还因为错过T-shirt伤心,今天第一次拿到了T-shirt高兴到爆啊!!接下来要努力刷有质量的题,希望能有所进步。
还是不废话了,说下这道题吧。一个数的因子数等于它完全分解为质数后,每一项质数个数加一连乘的结果。比如75=3*5*5,那么它的因子数为(1+1)*(2+1)=6。其实这个结论特别好证明,每一个因子数可以取0到它的个数个。比如,每个都取0个,那么结果即为1。如果每个都取,那么结果为那个数本身。预处理下1-100000的因子数,然后将读进来的数按因子数多少排个序,然后取其前n(当n<200)或者前200个计算,取最大值即可。
代码:
#include <iostream> #include <map> #include <cmath> #include <algorithm> using namespace std; map <int,int> store; bool is_prime(unsigned int x) { if(x==1)return false; else if(x==2)return true; int y=sqrt(1.0*x); for(int i=2;i<=y;i++) { if(x%i==0) return false; } return true; } void cal(unsigned int x) { if(is_prime(x)) store[x]++; else { int y=sqrt(1.0*x); for(int i=2;i<=y;i++) { if(x%i==0) { cal(i); cal(x/i); break; } } } } int ans[100005]; struct info { int key,val; }tempstore[100005]; bool cmp(info a,info b) { return a.val>b.val; } int main() { ans[1]=1; ans[2]=2; ans[3]=2; int tmp,res,maxx; for(int i=4;i<=100000;i++) { res=1; store.clear(); int up_limit=sqrt(1.0*i); tmp=i; for(int j=2;j<=up_limit;j++) { if(tmp%j==0) { cal(j); cal(tmp/j); break; } } if(store.size()==0) { ans[i]=2; continue; } map <int,int> ::iterator iter; for(iter=store.begin();iter!=store.end();iter++) res*=((iter->second)+1); ans[i]=res; } int n; while(cin>>n) { maxx=0; for(int i=0;i<n;i++) { cin>>tmp; tempstore[i].key=tmp; tempstore[i].val=ans[tmp]; } sort(tempstore,tempstore+n,cmp); int z=min(200,n); for(int i=0;i<z;i++) { for(int m=0;m<z;m++) { long long w=((long long int )tempstore[i].key)*tempstore[m].key; store.clear(); res=1; int up_limit=sqrt(1.0*w); //cout<<"*\n"; for(int j=2;j<=up_limit;j++) { if(w%j==0) { cal(j); cal((int)w/j); break; } } //cout<<"@\n"; if(store.size()==0) { res=2; } map <int,int> ::iterator iter; for(iter=store.begin();iter!=store.end();iter++) res*=((iter->second)+1); if(res>maxx) maxx=res; //cout<<"&\n"; } } cout<<maxx<<endl; } return 0; }
原文地址:http://blog.csdn.net/david_jett/article/details/45463043