2018徐州icpc M

對於每個燈求其探照範圍，用向量叉積來判斷是否可以照到，題目給的多邊形是逆時針的，設 燈爲A點，多邊形其中一邊按逆時針方向來講起始的點爲B結尾的點爲C
那麼（向量）AB x （向量）AC是負數就說明這邊能被這個燈照到，所有的燈照的邊都找到之後，就是貪心了，對於每個燈把他看做答案的起點，從這個燈照的範圍找到一個燈泡，使當前光照範圍最大，再從當前光照的區間找到下一個燈泡，使光照範圍最大，以此類推的貪心下去。
做了一些優化，只跑了717ms還挺快

#include <algorithm>
#include <iostream>
#include <stdio.h>
#include <string>
#include <cstring>
#include <vector>
#include <set>
#include <math.h>
#define eps 1e-14
#define pi acos(-1)
#define ll long long
#define RD T*(rand()*2-RAND_MAX)
#define Drand (long double)rand()/RAND_MAX
#define LINF 0x7f7f7f7f7f7f7f7f
#define INF 0x3f3f3f3f
using namespace std;
const int maxn=1e6+100;
const long long mod=1e18;
ll MOD(ll a,ll m){return a>m?a%m+m:a;}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}

struct nn
{
    int l,r,id;
}area[1005];

int cmp(nn a,nn b)
{
    return a.l<b.l;
}

struct Point
{
    ll x,y;
    Point(){}
    Point(ll x,ll y):x(x),y(y){}
}pol[1005],lig[1005];

struct Vt
{
    ll x,y;
    Vt(){}
    Vt(Point beg,Point ed)
    {
        x=ed.x-beg.x;
        y=ed.y-beg.y;
    }
};

ll cross(Vt a,Vt b)
{
    return a.x*b.y-a.y*b.x;
}

int cal(Point light,Point polb,Point pole)
{
    Vt vt1(light,polb),vt2(light,pole);
    ll flag=cross(vt1,vt2);
    return flag<0?1:0;
}

int main()
{
//    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
    int t;
    scanf("%d",&t);
    while(t--){
        memset(area,0,sizeof area);
        int n,m;
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)scanf("%lld%lld",&pol[i].x,&pol[i].y);
        for(int i=1;i<=m;i++)scanf("%lld%lld",&lig[i].x,&lig[i].y);
        pol[n+1]=pol[1];

        for(int i=1;i<=m;i++){
            int st=cal(lig[i],pol[1],pol[2]);
            int ed=cal(lig[i],pol[n],pol[n+1]);
            if(st==1 && ed==1){
                for(int j=1;j<=n;j++){
                    if(cal(lig[i],pol[j+1],pol[j+2])==0){
                        area[i].r=j+n;
                        break;
                    }
                }
                for(int j=n;j>=1;j--){
                    if(cal(lig[i],pol[j-1],pol[j])==0){
                        area[i].l=j;
                        break;
                    }
                }
            }
            else{
                for(int j=1;j<=n;j++){
                    if(area[i].l==0 && cal(lig[i],pol[j],pol[j+1])){
                        area[i].l=j;
                    }
                    if(area[i].l!=0 && cal(lig[i],pol[j],pol[j+1])){
                        area[i].r=area[i].r==0?j:max(area[i].r,j);
                    }
                }
            }
        }

        for(int i=1;i<=m;i++)area[i].id=i;

        sort(area+1,area+1+m,cmp);

//        for(int i=1;i<=m;i++){
//            cout<<" area "<<area[i].id<<" "<<area[i].l<<" "<<area[i].r<<endl;
//        }

        vector<int>tmp,ans;
        for(int i=1;i<=m;i++){
            int r=area[i].r;
            int j=i;
            int pos=i;
            int flag=0;
            tmp.clear();
            tmp.push_back(area[i].id);
            while(r-n<area[i].l-1 && j<=m){

                for(;j<=m && area[j].l<=r+1;j++){
                    if(area[j].r>area[pos].r){
                        pos=j;
                    }
                }
//                cout<<" i = "<<i<<"  pos = "<<pos<<endl;
                if(area[pos].r<=r){
                    flag=1;
                    break;
                }
                tmp.push_back(area[pos].id);
                r=area[pos].r;
            }
            if(area[pos].r-n<area[i].l-1)flag=1;
            if(flag==1)continue;
            if(ans.empty() || ans.size()>tmp.size())ans=tmp;
        }
        if(ans.empty()){
            printf("-1\n");
        }
        else{
            int len=ans.size();
            printf("%d\n%d",len,ans[0]);
            for(int i=1;i<len;i++){
                printf(" %d",ans[i]);
            }
            printf("\n");
        }
    }
    return 0;
}