背景 Background
NOIP 2012 提高组 题2描述 Description
恰逢 H 国国庆,国王邀请 n 位大臣来玩一个有奖游戏。首先,他让每个大臣在左、右手上面分别写下一个整数,国王自己也在左、右手上各写一个整数。然 后,让这 n 位大臣排成一排,国王站在队伍的最前面。排好队后,所有的大臣都会获得国王奖赏的若干金币,每位大臣获得的金币数分别是:排在该大臣前面的 所有人的左手上的数的乘积除以他自己右手上的数,然后向下取整得到的结果。国王不希望某一个大臣获得特别多的奖赏,所以他想请你帮他重新安排一下队伍的顺序,使得获得奖赏最多的大臣,所获奖赏尽可能的少。注意,国王的位置始终在队伍的最前面。
输入格式 InputFormat
第一行包含一个整数 n,表示大臣的人数。第二行包含两个整数 a 和 b,之间用一个空格隔开,分别表示国王左手和右手上的整数。
接下来 n 行,每行包含两个整数 a 和 b,之间用一个空格隔开,分别表示每个大臣左手和右手上的整数。
输出格式 OutputFormat
输出只有一行,包含一个整数,表示重新排列后的队伍中获奖赏最多的大臣所获得的金币数。样例输入 SampleInput [复制数据]
3 1 1 2 3 7 4 4 6
样例输出 SampleOutput [复制数据]
2
数据范围和注释 Hint
【输入输出样例说明】按 1、2、3 号大臣这样排列队伍,获得奖赏最多的大臣所获得金币数为 2;
按 1、3、2 这样排列队伍,获得奖赏最多的大臣所获得金币数为 2;
按 2、1、3 这样排列队伍,获得奖赏最多的大臣所获得金币数为 2;
按 2、3、1 这样排列队伍,获得奖赏最多的大臣所获得金币数为 9;
按 3、1、2 这样排列队伍,获得奖赏最多的大臣所获得金币数为 2;
按 3、2、1 这样排列队伍,获得奖赏最多的大臣所获得金币数为 9。
因此,奖赏最多的大臣最少获得 2 个金币,答案输出 2。
【数据范围】
对于 20%的数据,有 1≤ n≤ 10,0 < a、b < 8;
对于 40%的数据,有 1≤ n≤20,0 < a、b < 8;
对于 60%的数据,有 1≤ n≤100;
对于 60%的数据,保证答案不超过 109;
对于 100%的数据,有 1 ≤ n ≤1,000,0 < a、b < 10000。
思路:当时比赛时我只有50分。当然没想到正解,好像是N^3的伪DP
现在知道了该怎么做了。。讲讲为什么要按照a*b排序
我们假设一个方案
a b
c d
e f
g h
第三项的值:a*c/f
第四项的值:a*c*e/h
假设我们交换第三项和第四项
变成
a b
c d
g h
e f
第三项的值:a*c/h
第四项的值:a*c*g/f
如果我们交换能获得更优的值,那么我们需要做比较。。第一次的第四项的值显然>第二次第三项的值
那么主要是比较第一次第四项的值和第二次第四项的值
需要满足:a*c*e/h>a*c*g/f 整理一下 e*f>h*g
那么也就是说:如果前一项的两个数的乘积>后一项的两个数的乘积,那么总能交换它们得到更优的值,故应该按照乘积大小排序
注意要高精!!!压8位的话,数组开到600就够了
#include <cmath>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <cstdlib>
using namespace std;
typedef long long ll;
const int maxn=2210,ml=10200,maxl=600;
const ll mod=100000000;
struct qq
{
int a,b;
friend bool operator < (qq a,qq b)
{
return a.a*a.b<b.a*b.b;
}
} q[maxn];
ll c[ml],d[ml],add,ans[ml],t;
int vv,n;
void close()
{
exit(0);
}
bool judge()
{
for (int i=1;i<=maxl;i++)
if (d[i]>ans[i])
return true;
else
if (ans[i]>d[i])
return false;
return true;
}
void compare()
{
if (judge())
memcpy(ans,d,sizeof(d));
}
void mul(ll *a,int v)
{
add=0;
for (int i=maxl;i>=1;i--)
{
t=a[i]*v+add;
a[i]=t % mod;
add=t / mod;
}
}
void divide(ll *a,int div)
{
add=0;
for (int i=1;i<=maxl;i++)
{
t=add*mod+a[i];
a[i]=t / div;
add=t % div;
}
}
void print(ll *a)
{
int i;
for (i=1;i<maxl && a[i]==0;i++);
printf("%lld",a[i]);
if (i==maxl) printf("\n");
for (int j=i+1;j<=maxl;j++)
{
printf("%08lld",a[j]);
if (j==maxl)
printf("\n");
}
}
void work()
{
ans[maxl]=vv / q[1].b;
c[maxl]=vv;
for (int i=1;i<n;i++)
{
mul(c,q[i].a);
memcpy(d,c,sizeof(c));
divide(d,q[i+1].b);
/*
printf("c:");
printf("d:");
printf("--------------------\n");
print(c);
print(d);
*/
compare();
}
print(ans);
}
void init()
{
scanf("%d",&n);
int useless;
scanf("%d %d",&vv,&useless);
for (int i=1;i<=n;i++)
scanf("%d %d",&q[i].a,&q[i].b);
sort(q+1,q+n+1);
ans[maxl]=vv;
work();
}
int main ()
{
init();
close();
return 0;
}