hihoCoder挑战赛14 A,B,C题解

转载请注明出处： http://www.cnblogs.com/fraud/ ——by fraud

题目1 : 不等式

时间限制:10000ms
单点时限:1000ms
内存限制:256MB

描述

给定n个关于X的不等式，问最多有多少个成立。

每个不等式为如下的形式之一：

X < C

X <= C

X = C

X > C

X >= C

输入

第一行一个整数n。

以下n行，每行一个不等式。

数据范围:

1<=N<=50,0<=C<=1000

输出

一行一个整数，表示最多可以同时成立的不等式个数。

样例输入

4
X = 1
X = 2
X = 3
X > 0

样例输出

2

枚举所有值，每隔0.5一个，这样就能取到所有情况了。少了小数部分，WA了两下才发现。。。真是惨

/**
* code generated by JHelper
* More info: https://github.com/AlexeyDmitriev/JHelper * @author xyiyy @https://github.com/xyiyy
*/

#include <iostream>
#include <fstream>

//#####################
//Author:fraud
//Blog: http://www.cnblogs.com/fraud/ //#####################
//#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <sstream>
#include <ios>
#include <iomanip>
#include <functional>
#include <algorithm>
#include <vector>
#include <string>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <climits>
#include <cctype>

using namespace std;
#define XINF INT_MAX
#define INF 0x3FFFFFFF
#define mp(X, Y) make_pair(X,Y)
#define pb(X) push_back(X)
#define rep(X, N) for(int X=0;X<N;X++)
#define rep2(X, L, R) for(int X=L;X<=R;X++)
#define dep(X, R, L) for(int X=R;X>=L;X--)
#define clr(A, X) memset(A,X,sizeof(A))
#define IT iterator
#define ALL(X) (X).begin(),(X).end()
#define PQ std::priority_queue
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
typedef vector<PII> VII;
typedef vector<int> VI;

int getTheAnswer() {
return 42;
}
//
// Created by xyiyy on 2015/8/10.
//

#ifndef ICPC_P_HPP
#define ICPC_P_HPP

const double EPS = 1e-10;

int sig(double x) {
if (fabs(x) < EPS)return 0;
return x < 0 ? -1 : 1;
}

double add(double a, double b) {
if (fabs(a + b) < EPS * (fabs(a) + fabs(b)))return 0;
return a + b;
}

class P {
public:
double x, y;

P() { }

P(double x, double y) : x(x), y(y) { }

P  operator+(const P &p) {
return P(add(x, p.x), add(y, p.y));
}

P operator-(const P &p) {
return P(add(x, -p.x), add(y, -p.y));
}

P operator*(const double &d) {
return P(x * d, y * d);
}

P operator/(const double &d) {
return P(x / d, y / d);
}

double Dot(P p) {
return x * p.x + y * p.y;
}

double dot(P p) {
return add(x * p.x, y * p.y);
}

double Det(P p) {
return x * p.y - y * p.x;
}

double det(P p) {
return add(x * p.y, -y * p.x);
}

double abs() {
return sqrt(abs2());
}

double abs2() {
return dot(*this);
}

//绕原点旋转角度B（弧度值）产生新的点
P rot(double rad) {
return P(add(x * cos(rad), -y * sin(rad)), add(x * sin(rad), y * cos(rad)));
}

P rot90() {
return P(-y, x);
}

bool equals(P p) {
return compareTo(p) == 0;
}

int compareTo(P p) {
int b = sig(x - p.x);
if (b != 0)return b;
return sig(y - p.y);
}

};

//判断点是否在线段上
bool on_seg(P p1, P p2, P q) {
return (p1 - q).det(p2 - q) == 0 && (p1 - q).dot(p2 - q) <= 0;
}

//求两条线段的交点
P intersection(P p1, P p2, P q1, P q2) {
return p1 + (p2 - p1) * ((q2 - q1).det(q1 - p1) / (q2 - q1).det(p2 - p1));
}

//线段相交判定
bool crsSS(P p1, P p2, P q1, P q2) {
if (max(p1.x, p2.x) + EPS < min(q1.x, q2.x))return false;
if (max(q1.x, q2.x) + EPS < min(p1.x, p2.x))return false;
if (max(p1.y, p2.y) + EPS < min(q1.y, q2.y))return false;
if (max(q1.y, q2.y) + EPS < min(p1.y, p2.y))return false;
/*(if((p1 - p2).det(q1 - q2) == 0){
return (on_seg(p1,p2,q1) || on_seg(p1,p2,q2) || on_seg(q1,q2,p1) || on_seg(q1,q2,p2));
}else{
P r = intersection(p1,p2,q1,q2);
return on_seg(p1,p2,r) && on_seg(q1,q2,r);

}*/
return (p2 - p1).det(q1 - p1) * (p2 - p1).det(q2 - p1) <= 0
&& (q2 - q1).det(p1 - q1) * (q2 - q1).det(p2 - q1) <= 0;
}

//直线和线段相交判断
bool crsLS(P l1, P l2, P s1, P s2) {
return (s1 - l2).det(l1 - l2) * (s2 - l2).det(l1 - l2) <= 0;
}

//直线相交判断
//返回-1表示重合，0表示平行，1表示相交
int crsLL(P p1, P p2, P q1, P q2) {
if (sig((p1 - p2).det(q1 - q2)) != 0)return 1;
if (sig((p1 - q2).det(q1 - p2)) != 0)return 0;
return -1;
}

//直线和直线的交点
/*P isLL(P p1,P p2,P q1,P q2){
double d = (q2 - q1).det(p2 - p1);
if(sig(d)==0)return NULL;
return intersection(p1,p2,q1,q2);
}*/

//点到直线的垂足
P proj(P p1, P p2, P q) {
return p1 + ((p2 - p1) * ((p2 - p1).dot(q - p1) / (p2 - p1).abs2()));
}

//直线到点的距离
double disLP(P p1, P p2, P q) {
return fabs((p2 - p1).det(q - p1)) / (p2 - p1).abs();
}

//线段到点的距离
double disSP(P p1, P p2, P q) {
if ((p2 - p1).dot(q - p1) <= 0)return (q - p1).abs();
if ((p1 - p2).dot(q - p2) <= 0)return (q - p2).abs();
return disLP(p1, p2, q);
}

//圆和线段相交的判定
bool crsCS(P c, double r, P p1, P p2) {
return disSP(p1, p2, c) < r + EPS && (r < (c - p1).abs() + EPS || r < (c - p2).abs() + EPS);
}

//圆与圆相交的判定
bool crsCC(P c1, double r1, P c2, double r2) {
double dis = (c1 - c2).abs();
return dis < r1 + r2 + EPS && fabs(r1 - r2) < dis + EPS;
}

//四点共圆判定
/*bool onC(P p1,P p2,P p3,P p4){
P c = CCenter(p1,p2,p3);
if(c == NULL) return false;
return add((c - p1).abs2(), -(c - p4).abs2()) == 0;
}*/

//三点共圆的圆心
/*P CCenter(P p1,P p2,P p3){
//if(disLP(p1, p2, p3) < EPS)return NULL;//三点共线
P q1 = (p1 + p2) * 0.5;
P q2 = q1 + ((p1 - p2).rot90());
P s1 = (p3 + p2) * 0.5;
P s2 = s1 + ((p3 - p2).rot90());
return isLL(q1,q2,s1,s2);
}*/

//求两圆的极角 以p为中心
double polarangle(P p, P q) {
return atan2(q.y - p.y, q.x - p.x);
}

bool cmp_x(const P &p, const P &q) {
if (p.x != q.x) return p.x < q.x;
return p.y < q.y;
}

vector<P> qs;

void convex_hull(P *ps, int n) {
sort(ps, ps + n, cmp_x);
int k = 0;
qs.resize(2 * n);
for (int i = 0; i < n; qs[k++] = ps[i++]) {
while (k > 1 && (qs[k - 1] - qs[k - 2]).det(ps[i] - qs[k - 1]) < EPS)k--;
}
for (int i = n - 2, t = k; i >= 0; qs[k++] = ps[i--]) {
while (k > t && (qs[k - 1] - qs[k - 2]).det(ps[i] - qs[k - 1]) < EPS)k--;
}
qs.resize(k - 1);
}

//求凸包的直径
double convexDiameter() {
int qsz = qs.size();
if (qsz == 1)return 0;
if (qsz == 2) {
return (qs[0] - qs[1]).abs();
}
int i = 0, j = 0;
rep(k, qsz) {
if (!cmp_x(qs[i], qs[k]))i = k;
if (cmp_x(qs[j], qs[k])) j = k;
}
double res = 0;
int si = i, sj = j;
while (i != sj || j != si) {
res = max(res, (qs[i] - qs[j]).abs());
if ((qs[(i + 1) % qsz] - qs[i]).det(qs[(j + 1) % qsz] - qs[j]) < 0) i = (i + 1) % qsz;
else j = (j + 1) % qsz;
}
return res;
}

#endif //ICPC_P_HPP

P p[1010][2];

class TaskC {
public:
void solve(std::istream &in, std::ostream &out) {
int n;
in >> n;
rep(i, n)
rep(j, 2)in >> p[i][j].x >> p[i][j].y;
double ans = 0;
rep(i, n)
rep(l1, 2) {
P p1 = p[i][l1];
rep(j, n) {
if (i == j)continue;
rep(l2, 2) {
P p2 = p[j][l2];
double sum = p1.Det(p2) * 0.125;
p2 = p2 - p1;
rep(k, n) {
if (i == k || j == k)continue;
if (sig(sum) == 0)break;
int num = 0;
rep(l3, 2) {
if (p2.Det(p[k][l3] - p1) > 0);
else num++;
}
sum *= num/2.0;
}
ans += sum;
}
}
}
out << fixed << setprecision(10) << fabs(ans) << endl;
}
};

int main() {
std::ios::sync_with_stdio(false);
std::cin.tie(0);
TaskC solver;
std::istream &in(std::cin);
std::ostream &out(std::cout);
solver.solve(in, out);
return 0;
}

代码君
当然，对于这题你也可以对于每个点先求出所有点的极角，排个序然后来做。这样就可以从O(n^3)变成O(n^2logn)了

还剩下D题还是不会做，待以后再补