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;
}