QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#796499 | #9520. Concave Hull | xixu | Compile Error | / | / | C++20 | 5.1kb | 2024-12-01 19:55:16 | 2024-12-01 19:55:16 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
// #pragma GCC optimize(2)
#define int long long
// #define int __int128
#define double long long
// #define double __int128
// #define ull __int128
#define re read()
#define pr(x) print(x)
#define fup(a, b, c, d) for(int a = (b); a <= (c); a += (d))
#define fdown(a, b, c, d) for(int a = (b); a >= (c); a -= (d))
#define sq(x) (x) * (x)
// #define read() re
#define write(x) w(x)
typedef long long ll;
// typedef unsigned long long ull;
typedef pair<int , int> PII;
typedef map<int , int> MII;
const int inf = 0x3f3f3f3f, N = 2e6 + 10, M = 4e5 + 10, mod = 1e9 + 7;
const ll INF = 0x3f3f3f3f3f3f3f3f;
char *p1,*p2,buf[100000];
#define nc() (p1==p2 && (p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++)
int read()
{
int x=0,f=1;
char ch=nc();
while(ch<48||ch>57)
{
if(ch=='-')
f=-1;
ch=nc();
}
while(ch>=48&&ch<=57)
x=x*10+ch-48,ch=nc();
return x*f;
}
void write(__int128 x)
{
if(x<0)
putchar('-'),x=-x;
if(x>9)
write(x/10);
putchar(x%10+'0');
// printf("\n");
return;
}
__int128 s(double x1, double y1, double x2, double y2, double x3, double y3)
{
return ((__int128)x1 * y2 + (__int128)x2 * y3 + (__int128)x3 * y1 - (__int128)x1 * y3 - (__int128)x2 * y1 - (__int128)x3 * y2);
}
void solve(){
struct point{
double x,y;
point operator+(point a)const{
return {x+a.x,y+a.y};
}
point operator-(point a)const{
return {x-a.x,y-a.y};
}
point operator*(double a)const{
return {a*x,a*y};
}
double operator*(point a)const{
return x*a.y-y*a.x;
}
double operator&(point a)const{
return x*a.x+y*a.y;
}
bool operator==(point a){
if(x==a.x&&y==a.y)return true;
return false;
}
bool operator<(point a)const{
if(x!=a.x)return x<a.x;
return y<a.y;
}
point rot(double t) const { // 逆时针
return {x*cos(t)-y*sin(t),x*sin(t)+y*cos(t)};
}
point rot90() const {
return {-y,x};
}
double len2() const {
return x*x+y*y;
}
double len() const {
return sqrt(x*x+y*y);
}
};
auto dis=[&](point a,point b)->double{
return sqrt(sq(a.x-b.x)+sq(a.y-b.y));
};
auto check=[&](point a,point b,point c)->bool{
return ((a.x-b.x)*(a.y-c.y)-(a.y-b.y)*(a.x-c.x))<=0;
};
auto Convex_Hull=[&](vector<point>po,int si)->stack<point>{
sort(po.begin(),po.end());
stack<point>s;
int size=0;
s.push(po[0]);
for(int i=1;i<=si-1;i++){
point c=po[i];
while(size){
point b=s.top();
s.pop();
point a=s.top();
s.push(b);
if(check(a,b,c)){
size--;
s.pop();
}else break;
}
s.push(c);
size++;
}
size=0;
for(int i=si-2;~i;i--){
point c=po[i];
while(size){
point b=s.top();
s.pop();
point a=s.top();
s.push(b);
if(check(a,b,c)){
size--;
s.pop();
}else break;
}
s.push(c);
size++;
}
if(s.size()==1)s.push(po[0]);
return s;
};
int n;
n = read();
vector<point> po(n);
for(auto &[x,y] : po) {
x = read(), y = read();
}
// cout << 0;
// return ;
stack<point> stk = Convex_Hull(po, n);
if(n <= 3) {
w(-1);
printf("\n");
return ;
}
vector<point> v1, v2;
map<point, int> mp;
while(stk.size()) {
auto t = stk.top();
stk.pop();
if(!mp[t])
v1.push_back(t);
// cout << t.x << ' ' << t.y << '\n';
mp[t] ++;
}
if(v1.size() == n) {
w(-1);
printf("\n");
return ;
}
// cout << '\n';
for(auto t : po) {
// cout << t.x << ' ' << t.y << '\n';
if(mp[t]) continue ;
v2.push_back(t);
}
// cout << '\n';
// return ;
// cout << v2.size() << '\n';
stk = Convex_Hull(v2, v2.size());
v2.clear();
// cout << stk.size() << '\n';
// return ;
while(stk.size()) {
auto t = stk.top();
stk.pop();
v2.push_back(t);
}
__int128 mi = INF;
__int128 su = 0;
point p = v1[0];
fup(i, 1, v1.size() - 1, 1) {
auto x1 = v1[i - 1], x2 = v1[i];
for(auto x3 : v2) {
mi = min(mi, (__int128)abs(s(x1.x, x1.y, x2.x, x2.y, x3.x, x3.y)));
}
if(i == 1) continue ;
su += (__int128)abs(s(x1.x, x1.y, x2.x, x2.y, p.x, p.y));
}
// cout << mi << ' ' << su << '\n';
w(su - mi);
printf("\n");
}
signed main()
{
// cin.tie(0);
// cout.tie(0);
// ios::sync_with_stdio(false);
int _ = 1;
_ = read();
while(_ --)
{
solve();
}
}
詳細信息
answer.code: In member function ‘solve()::point solve()::point::rot(long long int) const’: answer.code:85:29: warning: narrowing conversion of ‘(((__gnu_cxx::__enable_if<true, double>::__type)(long long int)((const solve::point*)this)->solve::point::x * std::cos<long long int>(t)) - ((__gnu_cxx::__enable_if<true, double>::__type)(long long int)((const solve::point*)this)->solve::point::y * std::sin<long long int>(t)))’ from ‘__gnu_cxx::__enable_if<true, double>::__type’ {aka ‘double’} to ‘long long int’ [-Wnarrowing] 85 | return {x*cos(t)-y*sin(t),x*sin(t)+y*cos(t)}; | ~~~~~~~~^~~~~~~~~ answer.code:85:47: warning: narrowing conversion of ‘(((__gnu_cxx::__enable_if<true, double>::__type)(long long int)((const solve::point*)this)->solve::point::x * std::sin<long long int>(t)) + ((__gnu_cxx::__enable_if<true, double>::__type)(long long int)((const solve::point*)this)->solve::point::y * std::cos<long long int>(t)))’ from ‘__gnu_cxx::__enable_if<true, double>::__type’ {aka ‘double’} to ‘long long int’ [-Wnarrowing] 85 | return {x*cos(t)-y*sin(t),x*sin(t)+y*cos(t)}; | ~~~~~~~~^~~~~~~~~ answer.code: In function ‘void solve()’: answer.code:202:51: error: call of overloaded ‘abs(__int128)’ is ambiguous 202 | mi = min(mi, (__int128)abs(s(x1.x, x1.y, x2.x, x2.y, x3.x, x3.y))); | ~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ In file included from /usr/include/c++/13/cstdlib:79, from /usr/include/x86_64-linux-gnu/c++/13/bits/stdc++.h:42, from answer.code:1: /usr/include/stdlib.h:840:12: note: candidate: ‘int abs(int)’ 840 | extern int abs (int __x) __THROW __attribute__ ((__const__)) __wur; | ^~~ In file included from /usr/include/c++/13/cstdlib:81: /usr/include/c++/13/bits/std_abs.h:79:3: note: candidate: ‘constexpr long double std::abs(long double)’ 79 | abs(long double __x) | ^~~ /usr/include/c++/13/bits/std_abs.h:75:3: note: candidate: ‘constexpr float std::abs(float)’ 75 | abs(float __x) | ^~~ /usr/include/c++/13/bits/std_abs.h:71:3: note: candidate: ‘constexpr double std::abs(double)’ 71 | abs(double __x) | ^~~ /usr/include/c++/13/bits/std_abs.h:61:3: note: candidate: ‘long long int std::abs(long long int)’ 61 | abs(long long __x) { return __builtin_llabs (__x); } | ^~~ /usr/include/c++/13/bits/std_abs.h:56:3: note: candidate: ‘long int std::abs(long int)’ 56 | abs(long __i) { return __builtin_labs(__i); } | ^~~ answer.code:205:36: error: call of overloaded ‘abs(__int128)’ is ambiguous 205 | su += (__int128)abs(s(x1.x, x1.y, x2.x, x2.y, p.x, p.y)); | ~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ /usr/include/stdlib.h:840:12: note: candidate: ‘int abs(int)’ 840 | extern int abs (int __x) __THROW __attribute__ ((__const__)) __wur; | ^~~ /usr/include/c++/13/bits/std_abs.h:79:3: note: candidate: ‘constexpr long double std::abs(long double)’ 79 | abs(long double __x) | ^~~ /usr/include/c++/13/bits/std_abs.h:75:3: note: candidate: ‘constexpr float std::abs(float)’ 75 | abs(float __x) | ^~~ /usr/include/c++/13/bits/std_abs.h:71:3: note: candidate: ‘constexpr double std::abs(double)’ 71 | abs(double __x) | ^~~ /usr/include/c++/13/bits/std_abs.h:61:3: note: candidate: ‘long long int std::abs(long long int)’ 61 | abs(long long __x) { return __builtin_llabs (__x); } | ^~~ /usr/include/c++/13/bits/std_abs.h:56:3: note: candidate: ‘long int std::abs(long int)’ 56 | abs(long __i) { return __builtin_labs(__i); } | ^~~