【NOIP2016提高A组模拟9.15】Osu

news/2024/6/27 8:54:34 标签: NOIP2016提高, 9-15模拟, osu, 二分, DP

Description

在平面直角坐标系中有n个点会依次出现,其中第i个点在第ti秒出现在(xi,yi)。
你第0秒在(0,0),求至少踩k个点的情况下你移动的速度的最大值最小是多少。
答案输出a,b,c,表示 abc
n,k<=2000

Solution

既然是双最值,那么就二分喽。。。
输出答案肿么鬼畜肿么办,不怕,我们只有N^2个可能的答案。
二分在这里面进行就可以了,不过4*10^6的快排还真是有梦想。。。
不过,我们也可以二分实数,然后记录路径
最后在路径上找最大值输出
各种卡常任你选喽~

Code

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 2005
#define sqr(x) (x)*(x)
using namespace std;
typedef double db;
db dis[N][N],l,r,ans;
int n,k,a,b,c,p,q,t[N],x[N],y[N],f[N],g[N];
bool check(db x) {
    memset(f,0,sizeof(f));ans=0;
    fo(i,1,n) {
        fo(j,0,i-1) if (dis[j][i]<=x&&f[j]+1>f[i]) 
            f[i]=f[j]+1,g[i]=j;
        if (f[i]>=k) {
            for(int x=i;g[x];x=g[x]) 
                if (dis[g[x]][x]>ans) ans=dis[g[x]][x],p=g[x],q=x;
            return 1;
        }
    }
    return 0;
}
int gcd(int x,int y) {return y?gcd(y,x%y):x;}
int main() {
    scanf("%d%d",&n,&k);
    fo(i,1,n) scanf("%d%d%d",&t[i],&x[i],&y[i]);
    fo(i,0,n-1) fo(j,i+1,n) {
        dis[i][j]=sqrt(sqr(x[i]-x[j])+sqr(y[i]-y[j]))*1.0/(t[j]-t[i]);
        r=max(r,dis[i][j]);
    }r+=1e-4;
    while (r-l>1e-4) {
        db mid=(l+r)/2.0;
        if (check(mid)) r=mid;else l=mid;
    }
    a=1;b=sqr(x[p]-x[q])+sqr(y[p]-y[q]);c=t[q]-t[p];
    fo(i,2,sqrt(b)) while (!(b%(i*i))) a*=i,b/=i*i;
    int g=gcd(a,c);a/=g;c/=g;
    printf("%d %d %d",a,b,c);
}

http://www.niftyadmin.cn/n/942356.html

相关文章

BZOJ 1088 [SCOI2005]扫雷Mine 动态规划

Description 相信大家都玩过扫雷的游戏。那是在一个n*m的矩阵里面有一些雷&#xff0c;要你根据一些信息找出雷来。万圣节到了 &#xff0c;“余”人国流行起了一种简单的扫雷游戏&#xff0c;这个游戏规则和扫雷一样&#xff0c;如果某个格子没有雷&#xff0c;那么它里面的数…

2016.9.17NOIPA组模拟赛总结

前言 中秋模拟赛后续。。。 为什么在异端的下午&#xff1f; 2:10 迟到了10分钟。。。 然而还有许多更晚的老司机在路上。。。 一看题边认为这套题不可做 不是不可做&#xff0c;1,2题放原题有意义吗&#xff1f; T3又是什么辣鸡东西&#xff1f;随手贪心一下应该可以吧…

BZOJ 3036 绿豆蛙的归宿 期望动规

Description 随着新版百度空间的下线&#xff0c;Blog宠物绿豆蛙完成了它的使命&#xff0c;去寻找它新的归宿。 给出一个有向无环的连通图&#xff0c;起点为1终点为N&#xff0c;每条边都有一个长度。绿豆蛙从起点出发&#xff0c;走向终点。 到达每一个顶点时&#xff0c;如…

【NOIP2016提高A组模拟9.17】序列

Description 给出两个长度为N的序列A和B&#xff0c;每次可以选择一个i和j(1<i<j<N),使Ai~Aj都1然后模4&#xff0c;求最少的操作次数把A变成B。 N<10^5 Solution 显然我们可以求出一个a数组ai表示Ai要用多少次才能变成Bi 然后就相当于每次把一个区间的a全部-…

【GDOI2017模拟8.20】准备食物2

Description 此生无悔入东方&#xff0c;来世愿生幻想乡&#xff01; 我没入教&#xff0c;WorldWide_D&#xff08;出题人&#xff09;入了 古明地觉拥有的第三只眼&#xff0c;可以读取别人的内心想法&#xff0c;于是无论是妖怪&#xff0c;还是怨灵都为之感到恐惧&#…

BZOJ 1044 [HAOI2008]木棍分割 二分+动态规划

Description 有n根木棍, 第i根木棍的长度为Li,n根木棍依次连结了一起, 总共有n-1个连接处. 现在允许你最多砍断m个连 接处, 砍完后n根木棍被分成了很多段,要求满足总长度最大的一段长度最小, 并且输出有多少种砍的方法使得总长 度最大的一段长度最小. 并将结果mod 10007。。。 …

2016.9.24GDOI组模拟总结

前言 周六英语课不敢逃啊&#xff01;&#xff01; 于是就光荣地少掉了40分钟。。。 后记 好吧我承认这是一次很无聊的比赛。。。 写完前言就到后记了。。。 原因是T1大水&#xff0c;随手费用流就切了。。。 &#xff08;然而WerKeyTom_FTD说20000个点的费用流玄学&…

BZOJ 1876 [SDOI2009]SuperGCD 高精度 更相减损术

Description Sheng bill有着惊人的心算能力&#xff0c;甚至能用大脑计算出两个巨大的数的GCD&#xff08;最大公约 数&#xff09;&#xff01;因此他经常和别人比赛计算GCD。有一天Sheng bill很嚣张地找到了你&#xff0c;并要求和你比 赛&#xff0c;但是输给Sheng bill岂不…