【SinGuLaRiTy-1028】 Copyright (c) SinGuLaRiTy 2017. All Rights Reserved.
对于所有题目：Time Limit: 1s | Memory Limit: 256MB
[USACO2016 Jan] 愤怒的奶牛 (Angry Cows)
题目描述
在数轴x上摆放有n（2<=n<=50000）捆干草堆,没有任何两堆在同样的位置，所有的位置均为整数。你可以用弹弓射击射击数轴上的任意地点。如果你用弹弓以R的力度射击x处，那么该处会发生爆炸，爆炸的范围是以R为半径的圆形区域，所以它会使得[x-R,x+R]的所有干草堆同时发生爆炸。这些干草堆的爆炸半径是R-1。它们又会触发连锁反应，第三轮的爆炸的半径为R-2，依次递减。请选择最小的力度射击，使得所有的干草堆全部爆炸。
输入
第一行包含N。接下来N个整数，表示干草堆的位置。所有位置在[0,1000000000]内。
输出
输出最小的力度R，使得所有的干草堆发生爆炸。四舍五入保留一位小数。
样例数据
| 样例输入 | 样例输出 |
|
5 |
3.0 |
<样例解释>
如果以力度3射击坐标5，则坐标3，坐标8处的干草堆会发生爆炸，然后又会引爆坐标1和坐标10的干草堆，最后引爆坐标11处的干草堆。
解析
f[i]表示要炸掉i之前的所有干草堆，在i点最少需要的半径。
g[i]表示要炸掉i之后的所有干草堆，在i点最少需要的半径。
找到满足(j<i&&a[i]-a[j]>f[j]+1)的最后一个j，f[i]=min(a[i]-a[j],f[j+1]+1)，g数组同理。
最后枚举起始的射击区间，计算答案。
由于答案的小数部分只可能是0或者0.5，所以将所有数乘以2就可以避免小数问题，最后答案再除以2就可以了。
Code
#include<iostream> #include<cstdio> #include<cstdlib> #include<cmath> #include<cstring> #include<algorithm> #define ll long long #define N 50005 #define inf 2000000000 using namespace std; int n,a[N],f[N],g[N]; int main() { scanf("%d",&n); for(int i=1;i<=n;i++) { scanf("%d",&a[i]); a[i]*=2; } sort(a+1,a+n+1); for(int i=1;i<=n;i++) f[i]=g[i]=inf; int t=1; f[1]=0; for(int i=2;i<=n;i++) { while(t+1<i&&a[i]-a[t+1]>f[t+1]+2) t++; f[i]=min(a[i]-a[t],f[t+1]+2); } t=n; g[n]=0; for(int i=n-1;i>=1;i--) { while(t-1>i&&a[t-1]-a[i]>g[t-1]+2) t--; g[i]=min(a[t]-a[i],g[t-1]+2); } int ans=inf; for(int i=1,j=n;i<j;) { ans=min(ans,max((a[j]-a[i])/2,max(f[i],g[j])+2)); if(f[i+1]<g[j-1]) i++; else j--; } printf("%.1lf",(double)ans/2); return 0; }
[USACO2016 Jan] 无线电通信 (Radio Contact)
题目描述
农夫约翰和奶牛贝西要去寻找丢失的奶牛，为了彼此能联系对方，他们带着无线电通讯设备。不幸的是电池快没有电了。所以它们要尽量节省电能。农夫从位置（fx，fy）出发，一共走N步，贝西从位置（bx,by）出发，一共走M步。农夫的路线是由一个长度为N的字符串限制，字符串只出现’E’或’S’或’W’或’N’中，表示东南西北四个方向。农夫每一单位时间可以选择不动，或者按照限制走出一步。奶牛贝西也是如此。现在无线电设备每一单位时间消耗的电量等于他们的距离的平方。请问，他们走到终点，最少消耗多少电量？
输入
第一行两个整数n，m（1<=n,m<=1000）,第二行为fx，fy表示农夫的起始位置，第三行为bx,by，表示贝西的起始位置。
接下来有两个字符串，第一个字符串表示农夫的路线，第二个字符串表示贝西的路线。
他们的坐标总是在（0<=x,y<=1000）。北是y的正方向，东为x的正方向。
输出
输出一个整数，表示最少消耗的电量。
样例数据
| 样例输入 | 样例输出 |
|
2 7 |
28 |
解析
其实，这是一道水题。
设f[i][j]表示当前第一个人走了i步第二个人走了j步，那么就有f[i+1][j]=min(f[i+1][j],f[i][j]+dis())。然后f[i][j+1]和f[i+1][j+1]用同样的方法推出来。
Code
#include<cstdio> #include<iostream> #include<cstring> #define N 1024 typedef long long ll; using namespace std; int n,m,a[N][2],b[N][2]; ll f[N][N]; char s[N]; inline int dis(int a,int b,int c,int d) { return (c-a)*(c-a)+(b-d)*(b-d); } int main() { scanf("%d%d",&n,&m); cin>>a[1][0]>>a[1][1]>>b[1][0]>>b[1][1]; scanf("%s",s+1); for(int i=1;i<=n;i++) { if(s[i]=='N') a[i+1][0]=a[i][0],a[i+1][1]=a[i][1]+1; else if(s[i]=='S') a[i+1][0]=a[i][0],a[i+1][1]=a[i][1]-1; else if(s[i]=='W') a[i+1][0]=a[i][0]-1,a[i+1][1]=a[i][1]; else if(s[i]=='E') a[i+1][0]=a[i][0]+1,a[i+1][1]=a[i][1]; } scanf("%s",s+1); for(int i=1;i<=m;i++) { if(s[i]=='N') b[i+1][0]=b[i][0],b[i+1][1]=b[i][1]+1; else if(s[i]=='S') b[i+1][0]=b[i][0],b[i+1][1]=b[i][1]-1; else if(s[i]=='W') b[i+1][0]=b[i][0]-1,b[i+1][1]=b[i][1]; else if(s[i]=='E') b[i+1][0]=b[i][0]+1,b[i+1][1]=b[i][1]; } for(int i=1;i<=n+1;i++) for(int j=1;j<=m+1;j++) f[i][j]=2147483647; f[1][1]=0; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { f[i+1][j]=min(f[i+1][j],f[i][j]+dis(a[i+1][0],a[i+1][1],b[j][0],b[j][1])); f[i][j+1]=min(f[i][j+1],f[i][j]+dis(a[i][0],a[i][1],b[j+1][0],b[j+1][1])); f[i+1][j+1]=min(f[i+1][j+1],f[i][j]+dis(a[i+1][0],a[i+1][1],b[j+1][0],b[j+1][1])); } printf("%I64d",f[n+1][m+1]); return 0; }
[USACO2016 Jan] 熄灯 (Lights Out)
题目描述
谷仓是一个简单的多边形。它的每一条边或平行于x轴，或平行于y轴。从任意一个点沿顺时针方向走，横边和纵边是交替出现的。多边形一共有n个点，编号为1到n。其中1号点为出口。奶牛贝西已经完全记住了谷仓的地图。它从某一个点出发，只能沿着边界走，在开灯的情况，它能很快知道怎么走才能更快到达出口。可是有一天，灯熄灭了。它一下慌了神，然后它忘了自己在哪个点了。但是幸运的是它仍然记得整个谷仓的地图，而且它能够凭触觉知道当前点的内角有多大，它也能感觉到当前点是不是出口，而且经过一条边后也能精确地计算出该边的长度。现在，为了寻找出口，它决定采取这样的策略：沿着顺时针方向继续走下去，直到它能够判断出自己的位置，然后再选择距离最短的方向（顺时针或逆时针）走到出口。求最坏的情况下，贝西走到出口要比灯没坏的情况下多花多少时间。
输入
第一行包含N（4<=N<=200）,接下来N行每行包含两个整数，表示多边形的n个顶点（xi,yi），按顺时针的方向给出。这些点的范围均在[-100000,100000]
输出
输出贝西按照上述策略在最坏情况下比正常情况下要多花的时间。
样例数据
| 样例输入 | 样例输出 |
|
4 |
2 |
<样例解释>
贝西如果在出口处(Point 1)，肯定可以感知到，这绝对不会是最坏情况。现在考虑它在其他3个点时的情形：它当前可以感知到内角的大小。但因为内角都是90度，所以他无法确定自己位置。于是它按照既定的策略，顺时针走：
1.如果它在开始在点2，它需要走到点3，此时，它知道自己在哪里了。于是它找最短的路径，不管哪边，都是11.所以，它一共需要走12个单位。如果是开灯的情况下，只需要走10个单位。所以，它要多走两个单位。
2.如果它在点3，它要走11个单位。开灯的情况下也要走11个单位。
3.如果它在点4，它要走1个单位。开灯的情况下也要走1个单位。
于是最坏的情况下，它要多走两个单位。
解析
这道题，其实最重要的就是贝西判定当前的位置的方法：通过走过一段独一无二的路线。那么，我们怎么用代码进行实现呢?在这里，我用的是vector套pair，这样写也显得有些混乱；有一种用hash来判断的方法，相比之下还是好理解一点。(hash是一个很重要的东西，以后的字符串处理常常会用到) 大体思路就是通过用角度的大小(90度和270度)和边的长度来表达区分不同的路线，如果这个特征值一样，我们就认为贝西无法区分，也就找不到自己的位置；如果特征值不同，贝西就可以辨别出方向了。需要注意的是：数据的处理以及表达方式一定要确保能描述每一条路线的特征，而不会出现不同的路线有相同特征值的情况。
Code
#include<iostream> #include<vector> #include<set> #include<algorithm> #include<cstdio> #define MAXN 210 using namespace std; int opt[MAXN]; int main() { int N; scanf("%d",&N); vector<pair<long long,long long> > A(N); for(int i=0;i<N;i++) { cin>>A[i].first>>A[i].second; } vector<int> S(1,0); for(int i=0;i<N;i++) { int j=(i+1)%N; int k=(i+2)%N; S.push_back(abs(A[i].first-A[j].first)+abs(A[i].second-A[j].second)); if((A[i].first-A[j].first)*(A[k].second-A[j].second)-(A[k].first-A[j].first)*(A[i].second-A[j].second)>0) S.push_back(-1); else S.push_back(-2); } S.back()=0; for(int i=0;i<N;i++) { opt[i+1]=opt[i]+S[i*2+1]; } opt[N]=0; for(int i=N-1;i>=0;i--) { opt[i]=min(opt[i],opt[i+1]+S[i*2+1]); } multiset<vector<int> > st; for(unsigned int i=0;i<S.size();i+=2) { for(unsigned int ln=1;i+ln<=S.size();ln+=2) { st.insert(vector<int>(S.begin()+i,S.begin()+i+ln)); } } int result=0; for(unsigned int i=2;i+2<S.size();i+=2) { int ln; int cost=0; for(ln=1;;ln+=2) { if(st.count(vector<int>(S.begin()+i,S.begin()+i+ln))==1) break; cost+=S[i+ln]; } result=max(result,cost+opt[(i+ln)/2]-opt[i/2]); } cout<<result; return 0; }
[Vijos 1157] 分梨子
题目描述
Finley家的院子里有棵梨树，最近收获了许多梨子。于是，Finley决定挑出一些梨子，分给幼稚园的宝宝们。可是梨子大小味道都不太一样，一定要尽量挑选那些差不多的梨子分给孩子们，那些分到小梨子的宝宝才不会哭闹。每个梨子都具有两个属性值，Ai和Bi，本别表示梨子的大小和甜度情况。假设在选出的梨子中，两个属性的最小值分别是A0和B0。只要对于所有被选出的梨子i，都满足C1*(Ai-A0)+C2*(Bi-B0)≤C3（其中，C1、C2和C3都是已知的常数），就可以认为这些梨子是相差不多的，可以用来分给小朋友们。那么，作为幼稚园园长的你，能算出最多可以挑选出多少个梨子吗？
输入
第一行一个整数N（1≤N≤2000），表示梨子的总个数。 第二行三个正整数，依次为C1,C2和C3（C1,C2≤2000，C3≤10^9）。 接下来的N行，每行两个整数。第i行的两个整数依次为Ai和Bi。
输出
只有一个整数，表示最多可以选出的梨子个数。
样例数据
| 样例输入 | 样例输出 |
|
3 |
2 |
解析
神奇的题~
对于这道题，我们要采用"数形结合"的思想：每一个梨子不是有两个特征值Ai和Bi吗，我们要把它们看做一组坐标(Ai,Bi)，这样一来，每个梨子就变为了坐标轴上的一个点。我们来继续处理给出的题目信息：在选出的梨子里，有两个最小值，也就是说，对于所有我们选中的梨子，都满足：Ai>A0,Bi>B0，也就是说，所有的我们选中的梨子所代表的的点，都应该位于直线x=A0的右边，直线y=B0的上方。<是不是很"神奇"?>我们继续，题目又说了：每一个梨子(点)都满足C1*(Ai-A0)+C2*(Bi-B0)≤C3，我们用更"数学"的语言翻译一下，就变成：“每一个点(X,Y)都满足M*(X-a)+N(Y-b)≤C (M,N,C均为常数，我们在这里也暂且把a和b当做常数)”，嘿，这不是高中数学教材里的线性规划不等式吗？下面我们来画一画图(图中的点代表梨子)：
如图，三条直线围成的这个直角三角形内的点就是我们能够选择的梨子了。欢呼了？雀跃啦？不行，现在直接枚举A0、B0，就已经是O(n^2)了，然后再在三角形里面扫描满足条件的梨子，O(n^3)直接超时。于是乎，我们要想一个法子，让扫描梨子这一项的工程量减小一点。我们先将原来的不等式变为：C1Ai+C2Bi-C3<=C1A0+C2B0 。然后我们可以这么搞：先在外层来一个for循环确定A0，接下来我们再由低到高枚举B0，这时看一看我们刚刚得出的不等式，发现等式右边的值变大了，也就是说会有新的梨子满足这个不等式。——哎，有法子了！我们在枚举A0、B0之前预处理一下：1>记录每一个y轴(y=1,y=2,y=3......)上的梨子(点)的数量，2>将梨子(点)按C1Ai+C2Bi-C3排序。我们在枚举B0的过程中，先删去由于B0增大而从三角形下方剔除出去的梨子(由于有了<1>预处理，这个过程就是O(1)的了)，再加上新的满足不等式的值(<2>预处理也大大见减少了时间)，不断更新，然后取最大值，最后就能得到答案了。
Code
#include<iostream> #include<cstdio> #include<algorithm> #include<cstring>
using namespace std; const int MAXN=2005;
struct node { int a,b,s; }p[MAXN];
int c1,c2,c3; int ans; int n; inline bool cmp(const node&a,const node&b) { return a.s<b.s; } void init() { scanf("%d",&n); scanf("%d%d%d",&c1,&c2,&c3); for(int i=1;i<=n;i++) { scanf("%d%d",&p[i].a,&p[i].b); p[i].s=p[i].a*c1+p[i].b*c2; } } int main() { init(); sort(p+1,p+n+1,cmp); for(int i=1;i<=n;i++) { int a0=p[i].a; int b0=p[i].b; int cur=i; int tot=0; while(p[cur++].s<=c1*a0+c2*b0+c3&&cur<=n) tot++; ans=max(ans,tot); } cout<<ans<<endl; return 0; }
Time: 2017-07-21