题目大意:输入I和J,输出n,满足:1.I<=n<=J。2.在1的所有数中,n=(a1+a2+a3+...)/n最小,其中:ai|n,且:ai!=n。
Time Limit:2000MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64u
数据规模:1<=I<=J<=1000000。
理论基础:无。
题目分析:规模相对于时间来说,直接枚举O(n^2)可能会超时(我没有试过),大家可以去试一下,我在这里想说的是根据题目描述,如果区间内包含1,那么0/1=0,而其他的数都大于0,故直接输出1。
我们讨论没有1的情况,可分为两种:
1.区间内存在素数,那么解必为区间内最大的素数。
2.区间内不存在素数,这时我们再进行枚举即可。
代码如下:
#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
typedef double db;
#define DBG 1
#define maa (1<<31)
#define mii ((1<<31)-1)
#define sl(c) ((int)(c).size()) //取字符串长度;
#define forl(i, a, b) for(int i = (a); i < (b); ++i) //不带边界值循环,正序
#define forle(i, a, b) for(int i = (a); i <= (b); ++i) //带边界值循环,正序
#define forh(i, a, b) for(int i = (a); i > (b); --i) //不带边界值,逆序
#define forhe(i, a, b) for(int i = (a); i >= (b); --i) //带边界值,逆序
#define forlc(i, a, b) for(int i##_b = (b), i = (a); i < i##_b; ++i) //带别名的循环,用于不可改变值
#define forlec(i, a, b) for(int i##_b = (b), i = (a); i <= i##_b; ++i)
#define forgc(i, a, b) for(int i##_b = (b), i = (a); i > i##_b; --i)
#define forgec(i, a, b) for(int i##_b = (b), i = (a); i >= i##_b; --i)
#define forall(i, v ) forl(i, 0, sz(v)) //循环所有
#define forallc(i, v ) forlc(i, 0, sz(v))
#define forlla(i, v ) forhe(i, sz(v)-1, 0)
#define forls(i, n, a, b) for(int i = a; i != b; i = n[i]) //搜表用
#define rep(n) for(int repp = 0; repp < (n); ++repp)
#define repc(n) for(int repp_b = (n), repp = 0; repp < repp_b; ++repp)
#define rst(a, v) memset(a, v, sizeof a) //把字符v填充到a reset 重置
#define cpy(a, b) memcpy(a, b, sizeof a) //copy b 的sizeof(a)个字符到a
#define rstn(a, v, n) memset(a, v, (n)*sizeof((a)[0])) //把字符v填充到a[n]之前的字节
#define cpyn(a, b, n) memcpy(a, b, (n)*sizeof((a)[0])) //copy b 的 n 个字符到a
#define ast(b) if(DBG && !(b)) { printf("%d!!|\n", __LINE__); while(1) getchar(); } //调试
#define dout DBG && cout << __LINE__ << ">>| "
#define pr(x) #x"=" << (x) << " | "
#define mk(x) DBG && cout << __LINE__ << "**| "#x << endl
#define pra(arr, a, b) if(DBG) {\
dout<<#arr"[] |" <<endl; \
forlec(i, a, b) cout<<"["<<i<<"]="<<arr[i]<<" |"<<((i-(a)+1)%8?" ":"\n"); \
if((b-a+1)%8) puts("");\
} //数列查看
#define rd(type, x) type x; cin >> x //读数
inline int rdi() { int d; scanf("%d", &d); return d; }
inline char rdc() { scanf(" "); return getchar(); }
inline string rds() { rd(string, s); return s; }
inline db rddb() { db d; scanf("%lf", &d); return d; }
template<class T> inline bool updateMin(T& a, T b) { return a>b? a=b, true: false; }
template<class T> inline bool updateMax(T& a, T b) { return a<b? a=b, true: false; }
bool nprime[1000001]={true,true};
void psieve(int n)
{
int temp;
for(int i=2;i*i<=n;i++)
{
if(!nprime[i])
{
for(int j=i;(temp=j*i)<=n;j++)nprime[temp]=true;
}
}
}
int main()
{
int m,n,flag=1,ans=mii;
float Min=mii;
scanf("%d%d",&m,&n);
if(m==1)
{
printf("1\n");
exit(0);
}
else
{
psieve(n);
forle(i,m,n)
{
if(!nprime[i])
{
ans=i;
flag=0;
}
}
if(flag==0)
{
printf("%d\n",ans);
exit(0);
}
else
{
forle(j,m,n)
{
float sum=0;
for(int k=2;k*k<=j;k++)
{
if(j%k==0&&j/k!=k)sum=sum+k+j/k;
else sum=sum+k;
}
if(sum/j<Min)
{
Min=sum/j;
ans=j;
}
}
printf("%d\n",ans);
exit(0);
}
}
return 0;
}
其中:exit(0)<==>return 0;
by : Jsun_moon http://blog.csdn.net/Jsun_moon