[POI 2008][BZOJ 1132]Tro

这题我真是无能为力了

这题的做法还是挺简单的

枚举左下角的点做为原点，把其余点按极角排序    PS.是作为原点，如枚举到 k 时，对于所有 p[i] (包括p[k]) p[i]-=p[k] (此处为向量减法)

排序后满足 i<j 的两个向量 p[i] 和 p[j] 的叉积都是正数了

ΣΣp[i]×p[j] = ΣΣ(p[i].x*p[j].y-p[i].y*p[j].x) = Σ(p[i].x*Σp[j].y)-Σ(p[i].y*Σp[j].x)

计算叉积和的复杂度就从 O(n²) 降为了 O(n)

再加上枚举和排序的复杂度，总复杂度就是 O(n²logn)

但是我的代码似乎被 BZOJ 讨厌了%>_<%

本地测都是 A 的，校内 OJ 也 A 了，但BZOJ一交上去秒 WA 啊，55555555

求大爷指导……

这是错误的代码，求教做人

 #include <cstdio>
 #include <algorithm>
 const int size=;
 typedef struct point vector;
 typedef long long llint;

 namespace IOspace
 {
     inline int getint()
     {
         register int num=;
         register char ch;
         do ch=getchar(); while (ch<'' || ch>'');
         do num=num*+ch-'', ch=getchar(); while (ch>='' && ch<='');
         return num;
     }
     inline void putint(llint num, char ch='\n')
     {
         char stack[];
         register int top=;
         if (num==) stack[top=]='';
         for ( ;num;num/=) stack[++top]=num%+'';
         for ( ;top;top--) putchar(stack[top]);
         if (ch) putchar(ch);
     }
 }

 struct point
 {
     llint x, y;
     inline point() {}
     inline point(llint _x, llint _y):x(_x), y(_y) {}
     inline vector & operator += (vector v) {x+=v.x; y+=v.y; return *this;}
     inline vector & operator -= (vector v) {x-=v.x; y-=v.y; return *this;}
 };
 point p[size];
 inline llint operator * (vector a, vector b) {return a.x*b.y-a.y*b.x;}
 inline bool operator < (point a, point b) {return a*b>;}
 inline void swap(point & a, point & b) {point t=a; a=b; b=t;}

 int N;
 inline llint count(int);

 int main()
 {
     llint ans=;

     N=IOspace::getint();
     for (int i=;i<=N;i++) p[i].x=IOspace::getint(), p[i].y=IOspace::getint();

     for (int i=N;i>=;i--) ans+=count(i);

     bool b=ans&1LL;
     IOspace::putint(ans>>1LL, '.');
     IOspace::putint(b?:);

     return ;
 }
 inline llint count(int n)
 {
     llint ret=;

     int k=;
     for (int i=;i<=n;i++)
         if (p[i].x<p[k].x || p[i].x==p[k].x && p[i].y<p[k].y)
             k=i;
     swap(p[k], p[n]);

     for (int i=;i<=n;i++) p[i]-=p[n];

     std::sort(p+, p+n);

     vector s(, );
     for (int i=;i<=n;i++)
     {
         ret+=s*p[i];
         s+=p[i];
     }

     return ret;
 }

欲哭无泪的本傻系列

我再把正确的放上来，帮助理解题解:

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 using namespace std;

 #define maxn 30010
 typedef long long ll;
 struct xllend3
 {
     int x,y;
 } orz[maxn];

 int n,m;
 ll ans;

 bool cmp(const xllend3 &a,const xllend3 &b)
 {
     return a.x*b.y>a.y*b.x;
 }
 ll orzhzw(int n)
 {
     int i=;
     for (int j=;j<n;j++)
     if (orz[j].x<orz[i].x||(orz[j].x==orz[i].x && orz[j].y<orz[i].y))
         i=j;
     swap(orz[i],orz[n-]);
     for (int i=;i<n-;i++)
         orz[i].x-=orz[n-].x,orz[i].y-=orz[n-].y;
     sort(orz,orz+n-,cmp);
     ll gui=;
     ll sx=,sy=;
     for (int i=n-;i>=;i--)
     {
         gui+=(ll)orz[i].x*sy-(ll)orz[i].y*sx;
         sx+=orz[i].x;
         sy+=orz[i].y;
     }
     return gui;
 }

 int main()
 {
     scanf("%d",&n);
     for (int i=;i<n;i++)    scanf("%d%d",&orz[i].x,&orz[i].y);
     for (int i=n;i>;i--)    ans+=orzhzw(i);
     cout<<ans/;
     if (ans&)    cout<<".5"<<endl;else cout<<".0"<<endl;
     return ;
 }

Orz mxh1999