GYM 100090 B.Birthday Cake（尺取）

Description
给出一个n个点的凸多边形，先要选择其中两个不相邻的顶点把这个凸多边形分成两部分，要求这两部分面积差值最小
Input
第一行一整数n表示点数，之后顺时针输入这n个点的坐标x,y(4<=n<=200000,-1e8<=x,y<=1e8)
Output
输出所选的两个点的编号
Sample Input
4
0 0
1 0
1 1
0 1
Sample Output
1 3
Solution
尺取，预处理1~i点这部分的面积之和，每次对于一个左端点l，只要l~r点这部分面积没有超过整个凸多边形面积的一半r就一直右移，然后更新最小差值，r往后移一下再更新下最小差值（因为不知道在最优条件下两部分面积哪边大）
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 222222
ll s[maxn],x[maxn],y[maxn];
ll mul(ll x1,ll y1,ll x2,ll y2)
{
    return x1*y2-x2*y1;
}
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)scanf("%I64d%I64d",&x[i],&y[i]);
        s[1]=0;
        for(int i=2;i<=n;i++)s[i]=s[i-1]+mul(x[i-1],y[i-1],x[i],y[i]);
        ll sum=abs(s[n]+mul(x[n],y[n],x[1],y[1]));
        ll ans=sum;
        int j=2,l=1,r=1;
        for(int i=1;i<=n&&j<=n;i++)
        {
            ll temp=2ll*abs(s[j]-s[i]+mul(x[j],y[j],x[i],y[i]));
            while(j<=n&&sum>=temp)
            {
                if(sum-temp<ans)
                {
                    ans=sum-temp;
                    l=i,r=j;
                }
                j++;
                temp=2ll*abs(s[j]-s[i]+mul(x[j],y[j],x[i],y[i]));
            }
            if(j<=n)
            {
                temp=2ll*abs(s[j]-s[i]+mul(x[j],y[j],x[i],y[i]));
                if(temp-sum<ans)
                {
                    ans=temp-sum;
                    l=i,r=j;
                }
            }
            j--;
        }
        printf("%d %d\n",l,r);
    }
    return 0;
}