CodeForcesdiv1:995C - Leaving the Bar(随机算法+贪心)
For a vector →v=(x,y)v→=(x,y), define |v|=√x2+y2|v|=x2+y2.
Allen had a bit too much to drink at the bar, which is at the origin. There are nn vectors →v1,→v2,⋯,→vnv1→,v2→,⋯,vn→. Allen will make nn moves. As Allen's sense of direction is impaired, during the ii-th move he will either move in the direction →vivi→ or −→vi−vi→. In other words, if his position is currently p=(x,y)p=(x,y), he will either move to p+→vip+vi→ or p−→vip−vi→.
Allen doesn't want to wander too far from home (which happens to also be the bar). You need to help him figure out a sequence of moves (a sequence of signs for the vectors) such that his final position pp satisfies |p|≤1.5⋅106|p|≤1.5⋅106 so that he can stay safe.
The first line contains a single integer nn (1≤n≤1051≤n≤105) — the number of moves.
Each of the following lines contains two space-separated integers xixi and yiyi, meaning that →vi=(xi,yi)vi→=(xi,yi). We have that |vi|≤106|vi|≤106 for all ii.
Output a single line containing nn integers c1,c2,⋯,cnc1,c2,⋯,cn, each of which is either 11 or −1−1. Your solution is correct if the value of p=∑ni=1ci→vip=∑i=1ncivi→, satisfies |p|≤1.5⋅106|p|≤1.5⋅106.
It can be shown that a solution always exists under the given constraints.
3
999999 0
0 999999
999999 0
1 1 -1
1
-824590 246031
1
8
-67761 603277
640586 -396671
46147 -122580
569609 -2112
400 914208
131792 309779
-850150 -486293
5272 721899
1 1 1 1 1 1 1 -1
题意:给定一些向量,你可以改变它的符号,使得这些向量之和的长度小于1.5e6。
思路:考虑到每条边的长度小于1e6,所以任意两条边的向量和长度小于1.414e6<1.5e6。任意三条边的向量和(可以改变方向)的长度也成立,下面给出证明...
所以我们贪心, 保证向量和离原点更近. 然后下面的代码AC了,然后又被FST了.是因为每次我贪心原则是一样的.最后的结果有可能大于1.5e6. 我们需要加一些随机性, 多次贪心,直到结果满足题意.(这里随机的作用就算多次贪心,直到满足)
证明:每三个向量abc中都能找到两个向量合起来 <= 1e6,然后不断合并,最后只会剩下一个或者两个向量。
对于a和b,1,若ab小于60度,其中一个转向,则变为大于120度,<1e6;2,ab角度大于120度,则其向量和的长度小于1e6。3,角度在60到120度之间,那么第三边c,不论怎么放,都或与a或与b满足上面的关系,使得其向量和长度小于1e6,那么一直这样合并,最后只剩下一边或者两边,而两边的情况最坏为1.414e6
FST代码,WA79:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=;
const ll p=;
ll ans[maxn+],x[maxn+],y[maxn+];
ll X,Y,P;
int main()
{
int N,i,j; P=p*p;scanf("%d",&N);
for(i=;i<=N;i++){
scanf("%I64d%I64d",&x[i],&y[i]);
if(X>=&&Y>=&&x[i]<=&&y[i]<=) X+=x[i],Y+=y[i],ans[i]=;
else if(X>=&&Y>=&&x[i]>=&&y[i]>=) X-=x[i],Y-=y[i],ans[i]=-;
else if(X<=&&Y<=&&x[i]<=&&y[i]<=) X-=x[i],Y-=y[i],ans[i]=-;
else if(X<=&&Y<=&&x[i]>=&&y[i]>=) X+=x[i],Y+=y[i],ans[i]=;
else if((X+x[i])*(X+x[i])+(Y+y[i])*(Y+y[i])<(X-x[i])*(X-x[i])+(Y-y[i])*(Y-y[i])) X+=x[i],Y+=y[i],ans[i]=;
else X-=x[i],Y-=y[i],ans[i]=-;
}
for(i=;i<=N;i++) printf("%d ",ans[i]);
return ;
}
AC代码:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=;
const ll p=;
ll X,Y,P; int ans[maxn];
struct in{ ll x,y,id; }s[maxn];
int main()
{
int N,i,j; P=p*p;scanf("%d",&N);
for(i=;i<=N;i++){
scanf("%I64d%I64d",&s[i].x,&s[i].y); s[i].id=i;
}
while(true){
random_shuffle(s+,s+N+); X=Y=;
for(i=;i<=N;i++){
if((X+s[i].x)*(X+s[i].x)+(Y+s[i].y)*(Y+s[i].y)>(X-s[i].x)*(X-s[i].x)+(Y-s[i].y)*(Y-s[i].y)) X-=s[i].x,Y-=s[i].y,ans[s[i].id]=-;
else X+=s[i].x,Y+=s[i].y,ans[s[i].id]=;
}
if(X*X+Y*Y<=P) {
for(i=;i<=N;i++) printf("%d ",ans[i]); return ;
}
}
return ;
}
(div1E也是需要随机算法!!!
CodeForcesdiv1:995C - Leaving the Bar(随机算法+贪心)的更多相关文章
- Codeforces 996E Leaving the Bar (随机化)
题目连接:Leaving the Bar 题意:给你n个向量,你可以加这个向量或减这个向量,使得这些向量之和的长度小于1.5e6. 题解: 按照正常的贪心方法,最后的结果有可能大于1.5e6 .这里我 ...
- 浅浅地谈一下随机算法【poj2454】【poj3318】
随机算法我也只是稍微接触了一下,就是想写篇博客自己稍微总结一下 其实随机算法也算是一个玄学吧,运气不好还是会wa.但是我们知道,计算机可以在短时间内计算大量的数据,所以碰到正确答案的概率还是挺大的. ...
- 微信红包中使用的技术:AA收款+随机算法
除夕夜你领到红包了吗?有的说“我领了好几K!”“我领了几W!” 土豪何其多,苦逼也不少!有的说“我出来工作了,没压岁钱了,还要发红包”.那您有去抢微信红包吗?微信群中抢“新年红包”春节爆红.618微信 ...
- POJ 3318 Matrix Multiplication(随机算法)
题目链接 随机算法使劲水...srand((unsigned)time(0))比srand(NULL)靠谱很多,可能是更加随机. #include <cstdio> #include &l ...
- 抽奖随机算法的技术探讨与C#实现
一.模拟客户需求 1.1 客户A需求:要求每次都按照下图的概率随机,数量不限,每个用户只能抽一次,抽奖结果的分布与抽奖概率近似. 1.2 客户B需求:固定奖项10个,抽奖次数不限,每个用户只能抽一次, ...
- hdu 4712 (随机算法)
第一次听说随机算法,在给的n组数据间随机取两个组比较,当随机次数达到一定量时,答案就出来了. #include<stdio.h> #include<stdlib.h> #inc ...
- 权重随机算法的java实现
一.概述 平时,经常会遇到权重随机算法,从不同权重的N个元素中随机选择一个,并使得总体选择结果是按照权重分布的.如广告投放.负载均衡等. 如有4个元素A.B.C.D,权重分别为1.2.3.4,随机结果 ...
- hdu 4712 Hamming Distance ( 随机算法混过了 )
Hamming Distance Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others) ...
- HDU4712+随机算法
随机算法 求n个20位的2进制串的MinDist. Dist:两个串的异或结果中1的个数 /* 随机算法 */ #include<algorithm> #include<iostre ...
随机推荐
- Shell脚本笔记 1
函数别名 设置别名 alias name="command" alias ll="ls -laS" 取消别名 unalias name 求取数学表达式 valu ...
- 一起学android之怎样设置TextView中不同字段的字体颜色(22)
在这里先看看效果图: OK,有时候,在我们的项目中会要求TextView中文本有一部分的字体颜色不一样.这时我们应该使用 SpannableStringBuilder这个工具类,当然这个类的功能非常强 ...
- 微服务之旅:从Netflix OSS到 Istio Service Mesh
在这篇文章中,我们从Netflix开始,通过Envoy和Istio的崛起,快速浏览微服务的历史. 微服务是具有边界上下文的松散耦合服务,使您能够独立开发,部署和扩展服务.它还可以定义为构建独立开发和部 ...
- 图像处理之基础---卷积及其快速算法的C++实现
头文件: /* * Copyright (c) 2008-2011 Zhang Ming (M. Zhang), zmjerry@163.com * * This program is free so ...
- freopen - C/C++文件输入输出利器
freopen以前经常使用,比较方便,可以当作模板,在中间替换为自己的代码即可使用. ? 1 2 3 4 5 6 7 8 9 10 11 12 13 14 #include <stdio.h&g ...
- win7激活附带激活软件
链接: https://pan.baidu.com/s/1i46yoHR 密码: 7k6y
- 【BZOJ1018】[SHOI2008]堵塞的交通traffic 线段树
[BZOJ1018][SHOI2008]堵塞的交通traffic Description 有一天,由于某种穿越现象作用,你来到了传说中的小人国.小人国的布局非常奇特,整个国家的交通系统可以被看成是一个 ...
- rtsp转rtmp、hls网页直播服务器EasyNVR前端兼容性调试:ie下的 pointer-events- none
发现问题: 之前在做EasyNVR 的web页面开发过程中,力求的都是一个播放效果的.功能的展示.对于兼容性也有注意,但有些细节还是难免有所疏忽. 内部测试发现:由于我们是流媒体的实时视频直播,在we ...
- Vue中表单校验
1.安装校验插件vee-validate npm install vee-validate --save 2.在main.js中引用插件 // 表单校验 import VeeValidate, { V ...
- sql server charindex函数和patindex函数详解(转)
charindex和patindex函数常常用来在一段字符中搜索字符或字符串.假如被搜索的字符中包含有要搜索的字符,那么这两个函数返回一个非零的整数,这个整数是要搜索的字符在被搜索的字符中的开始位数. ...