【NOIP2016提高A組模擬9.15】Osu
阿新 • • 發佈:2018-05-17
oid std lib 最大 algorithm ios d+ -a gin
題目
分析
考慮二分答案,
二分小數顯然是不可取的,那麽我們將所有可能的答案求出來,記錄在一個數組上,排個序(C++調用函數很容易超時,手打快排,時間復雜度約為\(O(>8*10^7)\),但相信夢想的力量)。
剩下就簡單了,將二分出的值判斷是否可以獲得k分以上,
這裏可以用多種方法,spfa、dp
dp:
\(dp_i\)表示移動到了第i個點的最大分數
#include <cmath> #include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <queue> const long long maxlongint=2147483647; const long long mo=1000000007; const long long N=2005; using namespace std; long long ans,a[N][3],n,m,tot,f[N]; double sum[N][N]; struct ddx { long long a,b,c; double ans; }d[N*N]; double cale(ddx x) { return 1.0*x.a*sqrt(x.b)/x.c; } bool cmp(ddx x,ddx y) { return x.ans<y.ans; } void q(int l,int r) { int i=l,j=r,e1; double mid=d[(l+r)/2].ans; double e; while(i<j) { while(d[i].ans<mid) i++; while(d[j].ans>mid) j--; if(i<=j) { e1=d[i].a; d[i].a=d[j].a; d[j].a=e1; e1=d[i].b; d[i].b=d[j].b; d[j].b=e1; e1=d[i].c; d[i].c=d[j].c; d[j].c=e1; e=d[i].ans; d[i].ans=d[j].ans; d[j].ans=e; i++; j--; } } if(i<r) q(i,r); if(l<j) q(l,j); } bool ok(double x) { memset(f,0,sizeof(f)); for(long long i=1;i<=n;i++) { for(long long j=0;j<=i-1;j++) { if(sum[i][j]<=x) { f[i]=max(f[i],f[j]+1); } } if(f[i]>=m) return true; } return false; } long long gcd(long long x,long long y) { return x==0?y:gcd(y%x,x); } int work(int x) { for(long long i=sqrt(d[x].b);i>=1;i--) { if(d[x].b%(i*i)==0) { d[x].b/=i*i; d[x].a*=i; break; } } long long k=gcd(d[x].a,d[x].c); if(!k) return 0; d[x].a/=k; d[x].c/=k; } double rf(long long l,long long r) { while(l<r) { long long mid=(l+r)/2; if(ok(d[mid].ans)) r=mid; else l=mid+1; } work(l); printf("%lld %lld %lld",d[l].a,d[l].b,d[l].c); } int main() { scanf("%lld%lld",&n,&m); for(long long i=1;i<=n;i++) for(long long j=0;j<=2;j++) { scanf("%lld",&a[i][j]); } for(long long i=0;i<=n-1;i++) for(long long j=i+1;j<=n;j++) { d[++tot].a=1; d[tot].b=(a[j][1]-a[i][1])*(a[j][1]-a[i][1])+(a[j][2]-a[i][2])*(a[j][2]-a[i][2]); d[tot].c=a[j][0]-a[i][0]; d[tot].ans=cale(d[tot]); sum[i][j]=sum[j][i]=d[tot].ans; } q(1,tot); rf(1,tot); }
【NOIP2016提高A組模擬9.15】Osu