[NOIP2012][贪心][高精度]国王游戏

恰逢H国国庆,国王邀请n位大臣来玩一个有奖游戏。首先,他让每个大臣在左、右手上面分别写下一个整数,国王自己也在左、右手上各写一个整数。然后,让这n位大臣排成一排,国王站在队伍的最前面。排好队后,所有的大臣都会获得国王奖赏的若干金币,每位大臣获得的金币数分别是:排在该大臣前面的所有人的左手上的数的乘积除以他自己右手上的数,然后向下取整得到的结果。
国王不希望某一个大臣获得特别多的奖赏,所以他想请你帮他重新安排一下队伍的顺序,使得获得奖赏最多的大臣,所获奖赏尽可能的少。注意,国王的位置始终在队伍的最前面。
格式
输入格式

第一行包含一个整数n,表示大臣的人数。
第二行包含两个整数a和b,之间用一个空格隔开,分别表示国王左手和右手上的整数。接下来n行,每行包含两个整数a和b,之间用一个空格隔开,分别表示每个大臣左手和右手上的整数。
输出格式

输出只有一行,包含一个整数,表示重新排列后的队伍中获奖赏最多的大臣所获得的金币数。
样例1
样例输入1[复制]

3
1 1
2 3
7 4
4 6
样例输出1[复制]

2
限制
每个测试点1s
提示
对于20%的数据,有1≤ n≤ 10,0 < a、b < 8;
对于40%的数据,有1≤ n≤20,0 < a、b < 8;
对于60%的数据,有1≤ n≤100;
对于60%的数据,保证答案不超过10^9;
对于100%的数据,有1 ≤ n ≤1,000,0 < a、b < 10000。
来源
Noip2012提高组复赛Day1T2

临近NOIP,来打一打高精度
看到最大最小会以为是二分吧,但幸好这个明显没有单调性,所以就可以直接丢了
考虑贪心的证明
首先对于相邻的i,j他们任意交换位置不会对前面的大臣和后面的大臣产生影响
设i前面大臣累乘结果为sum,i(Li,Ri),j(Lj,Rj)
①i在前面
ans=max(sum/Ri,sum*Li/Rj)
②j在前面
ans=max(sum/Rj,sum*Lj/Ri)
显然的,①的后一项比②的前一项大,②的后一项比①的前一项大,么只要②的后一项比①的后一项大,一定选择i在前面更优,这个条件就可以推出,Li*Ri

#include
#include
#include
#include
using namespace std;
const int N=1e3;
int n;
int readint()
{
    int x=0; char ch=getchar();
    while (!isdigit(ch)) ch=getchar();
    while (isdigit(ch)) {x=x*10+ch-'0'; ch=getchar();}
    return x;
}
struct Minister
{
    int A,B;
    bool operator<(const Minister &x)const
    {
        return A*B1];
#define A(x) minister[x].A
#define B(x) minister[x].B
struct hp
{
    int w[50000];
}ans,temp,sum;
bool compare(hp a,hp b)
{
    if (a.w[0]>b.w[0]) return false;
    if (a.w[0]0]) return true;
    for (int i=a.w[0];i;i--) 
      if (a.w[i]return true;
      else if (a.w[i]>b.w[i]) return false;
    return false;
}
hp multi(hp a,int b)
{
    hp ret;
    memset(ret.w,0,sizeof(ret.w));
    ret.w[0]=a.w[0]+5;
    for (int i=1;i<=a.w[0];i++)
    {
        ret.w[i]+=a.w[i]*b;
        ret.w[i+1]+=ret.w[i]/10;
        ret.w[i]%=10;
    }
    for (int i=a.w[0]+1;i<=ret.w[0];i++)
    {
        ret.w[i+1]+=ret.w[i]/10;
        ret.w[i]%=10;
    }
    while (ret.w[0]>1&&!ret.w[ret.w[0]]) ret.w[0]--;
    return ret;
}
hp division(hp a,int b)
{
    hp ret;
    memset(ret.w,0,sizeof(ret.w));
    ret.w[0]=a.w[0];
    int x=0;
    for (int i=a.w[0];i;i--)
    {
        x=x*10+a.w[i];
        ret.w[i]=x/b;
        x%=b;
    }
    while (ret.w[0]>1&&!ret.w[ret.w[0]]) ret.w[0]--;
    return ret;
}
int main()
{
//  freopen("king.in","r",stdin);
    n=readint();
    for (int i=0;i<=n;i++) A(i)=readint(),B(i)=readint();
    sort(minister+1,minister+n+1);
    ans.w[0]=1; ans.w[1]=0;
    sum.w[0]=1; sum.w[1]=1;
    sum=multi(sum,A(0));
    for (int i=1;i<=n;i++)
    {
        temp=division(sum,B(i));
        if (compare(ans,temp)) ans=temp;
        sum=multi(sum,A(i));
    }
    for (int i=ans.w[0];i;i--) printf("%d",ans.w[i]);
    return 0; 
}

你可能感兴趣的:(贪心,高精度)