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

 

终于滚回来填坑了，高精度数组要开大点，我压位开500都还太小，查了好久的错。。。

重新更新了一遍高精度模板

 

国王游戏贪心证明：

设两人左右手的数字为(a,b)和(x,y)，S为两人之前所有人左手数字的乘积

ans=max{(S*a)/y,(S*x)/b}

设(S*a)/y<(S*x)/b

显然a/y<x/b

即a*b<x*y

证毕

#include<iostream>
#include<algorithm>
#include<cstring>
#include<iomanip>
using namespace std;

const int MAXN=1000;
const int DLEN=4;
const int WIDE=10000;
class BigNum
{
public:
    int NUM[MAXN];
    int L;
    bool flag;
    BigNum(){memset(NUM,0,sizeof(NUM));L=1;flag=0;}
    BigNum(const BigNum &T){memcpy(NUM,T.NUM,sizeof(NUM));L=T.L;flag=T.flag;}
    BigNum(int n){memset(NUM,0,sizeof(NUM));NUM[0]=n;L=1;while(NUM[L-1]>=WIDE){NUM[L]+=NUM[L-1]/WIDE;NUM[L-1]%=WIDE;L++;}flag=0;}
};

void Input(string s,BigNum &T)
{
    int k=1,num=0;
    memset(T.NUM,0,sizeof(int)*MAXN);T.L=0;
    for(int i=s.size()-1;i>=0;i--)
    {
        if(k==WIDE) T.NUM[T.L++]=num,num=0,k=1;
        num+=k*(s[i]-'0');
        k*=10;
    }
    if(num>0) T.NUM[T.L++]=num;
}

void Output(const BigNum T)
{
    if(T.flag==1) cout<<'-';
    cout<<T.NUM[T.L-1];
    for(int i=T.L-2;i>=0;i--)
    {
        cout.width(DLEN);
        cout.fill('0');
        cout<<T.NUM[i];
    }
}

bool cmp(const BigNum A,const BigNum B)
{
    if(A.L!=B.L) return A.L<B.L;
    for(int i=A.L-1;i>=0;i--)
        if(A.NUM[i]!=B.NUM[i])
            return A.NUM[i]<B.NUM[i];
    return 0;
}

/*BigNum Add(const BigNum A,const BigNum B)
{
    BigNum C;
    int L=max(A.L,B.L);C.L=L;
    for(int i=0;i<L;i++)
    {
        C.NUM[i]+=A.NUM[i]+B.NUM[i];
        if(C.NUM[i]>=WIDE)
            C.NUM[i]-=WIDE,C.NUM[i+1]++;
    }
    if(C.NUM[L]>0) C.L++;
    return C;
}

BigNum Add(const BigNum A,int B)
{
    BigNum C(A);
    C.NUM[0]+=B;
    for(int i=0;i<C.L;i++)
        if(C.NUM[i]>=WIDE)
            C.NUM[i]-=WIDE,C.NUM[i+1]++;
    if(C.NUM[C.L]>0) C.L++;
    return C;
}

BigNum Dec(BigNum A,BigNum B)
{
    BigNum *X=&A,*Y=&B,C;
    int L=max(X->L,Y->L);
    C.L=L;
    if(cmp(*X,*Y)) swap(X,Y),C.flag=1;
    for(int i=0;i<L;i++)
    {
        C.NUM[i]+=X->NUM[i]-Y->NUM[i];
        if(C.NUM[i]<0)
            C.NUM[i]+=WIDE,C.NUM[i+1]--;
    }
    while(C.NUM[C.L-1]==0) C.L--;
    return C;
}

BigNum Mult(const BigNum A,const BigNum B)
{
    BigNum C;
    for(int i=0;i<A.L;i++)
        for(int j=0;j<B.L;j++)
        {
            C.L=i+j;
            C.NUM[C.L]+=A.NUM[i]*B.NUM[j];
            if(C.NUM[C.L]>=WIDE)
                C.NUM[C.L+1]+=C.NUM[C.L]/WIDE,C.NUM[C.L]%=WIDE;
        }
    C.L=A.L+B.L;
    if(C.NUM[C.L-1]==0) C.L--;
    return C;
}*/

BigNum Mult(const BigNum A,int B)
{
    BigNum C(A);
    int i,tmp,k=0;
    for(i=0;i<C.L||k;i++)
    {
        tmp=C.NUM[i]*B+k;
        k=tmp/WIDE;
        C.NUM[i]=tmp%WIDE;
    }
    C.L=i;
    return C;
}

/*BigNum Div2(const BigNum A)
{
    BigNum C(A);
    for(int i=C.L-1;i>=0;i--)
    {
        if(C.NUM[i]%2&&i>0)
            C.NUM[i-1]+=WIDE;
        C.NUM[i]/=2;
    }
    if(C.NUM[C.L-1]==0) C.L--;
    return C;
}

BigNum Div(const BigNum A,const BigNum B)
{
    BigNum L(1),R(A),T;
    while(cmp(Add(L,1),R))
    {
        T=Div2(Add(L,R));
        if(cmp(A,Mult(T,B)))
            R=T;
        else 
            L=T;
    }
    return L;
}*/

BigNum Div(const BigNum A,int B)
{
    BigNum C(A);
    int k=0;
    for(int i=C.L-1;i>=0;i--)
    {
        k=k*WIDE+C.NUM[i];
        C.NUM[i]=k/B;
        k%=B;
    }
    while(C.NUM[C.L-1]==0) C.L--;
    return C;
}

const int MAXN_=1000+10;
struct Num
{
    int L,R;
}T[MAXN_];
int n;
BigNum x=1,maxn=1;
BigNum t;

bool cmp_(Num A,Num B)
{
    return A.L*A.R<B.L*B.R;
}

int main()
{
    cin>>n;
    cin>>T[0].L>>T[0].R;
    for(int i=1;i<=n;i++)
        cin>>T[i].L>>T[i].R;
    sort(&T[1],&T[n+1],cmp_);
    for(int i=1;i<=n;i++)
    {
        x=Mult(x,T[i-1].L);
        t=Div(x,T[i].R);
        if(cmp(maxn,t))
            maxn=t;
    }
    Output(maxn);
    return 0;
}