博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
P1080 国王游戏
阅读量:4693 次
发布时间:2019-06-09

本文共 2090 字,大约阅读时间需要 6 分钟。

考虑任意一个排队方案,对于其中某两个相邻位置 $i>0,j=i+1$,如果交换更优

那么有 $max(A/r[i],Al[i]/r[j])>max(A/r[j],Al[j]/r[i])$,其中 $A=\prod_{k=0}^{i-1}l[k]$,$l[0]$ 是国王左手的数

因为 $A/r[i]<=Al[j]/r[i]$,$A/r[j]<=Al[i]/r[j]$,分类讨论一波发现上式等价于 $Al[i]/r[j]>Al[j]/r[i]$

得到 $l[i]r[i]>l[j]r[j]$ 时,交换最优,所以直接按 $l*r$ 为关键字排序即可

如果不太理解 $max$ 是怎么没的,直接按 $max(1/r[i],l[i]/r[j])>max(1/r[j],l[j]/r[i])$ 排序也行

计算答案时要用高精度,重新学了一遍压位高精

#include
#include
#include
#include
#include
using namespace std;typedef long long ll;const ll wid=100000000;const int N=2e5+7;char ch[N];struct bigint{ ll a[12507],len; bigint() { memset(a,0,sizeof(a)); len=0; } inline void read() { scanf("%s",ch+1); ll l=strlen(ch+1),t=0,k=1; len=(l+7)/8; for(int i=l;i;i--) { if(!((l-i)%8)){ k=1; t++; } a[t]+=k*(ch[i]^48); k*=10; } } inline void print() { if(!len) { printf("0\n"); return; } printf("%lld",a[len]); for(int i=len-1;i;i--) printf("%08lld",a[i]); printf("\n"); } inline bool operator < (const bigint &tmp) const { if(len!=tmp.len) return len
1) u.len--; return u; } inline bigint operator + (const bigint &tmp) const { bigint u; ll x=0; u.len= len>tmp.len ? len : tmp.len; for(int i=1;i<=u.len;i++) { u.a[i]=a[i]+tmp.a[i]+x; x=u.a[i]/wid; u.a[i]%=wid; } while(x) u.a[++u.len]=x%wid,x/=wid; return u; } inline bigint operator * (const bigint &tmp) const { bigint u; for(int i=1;i<=tmp.len;i++) for(int j=1;j<=len;j++) u.a[i+j-1]+=tmp.a[i]*a[j];//先不考虑进位 u.len=len+tmp.len; for(int i=1;i<=u.len;i++) u.a[i+1]+=u.a[i]/wid,u.a[i]%=wid;//最后统一进位 while(!u.a[u.len]&&u.len>1) u.len--; return u; } inline bigint operator / (const int &tmp) const { //整除 bigint u; u.len=len; ll x=0; for(int i=len;i;i--) { x*=wid; x+=a[i]; u.a[i]=x/tmp; x%=tmp; } while(!u.a[u.len]&&u.len>1) u.len--; return u; }}A,B,C,ans;int n;struct dat{ int x,y; inline bool operator < (const dat &tmp) const { return x*y

 

转载于:https://www.cnblogs.com/LLTYYC/p/11367907.html

你可能感兴趣的文章