对于每个灯求其探照范围,用向量叉积来判断是否可以照到,题目给的多边形是逆时针的,设 灯为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;
}