[ARC082-C][几何]ConvexScore

Description

题解

栋老师的神仙题…
话说栋老师的题现在我似乎还是基本不会做啊
这么菜拿头去考试啊…
考虑一个可以构成凸包的点集$S$，其内部包括边界上的点为$T$
那么贡献给出是$2^{|T|-|S|}$
实际上就是要计算$(S,U)$这样的有序对数量，其中$S$表示凸包点集，$U$表示子集，允许空集
然后再想一想，突然就会发现
$S\cup U$的集合互不相同
于是就是要计算有面积的点集个数
只会有空集和共线存在，$n^3$枚举每条线一下就可以了
upd:
操没写这个凸包dp直接被thu打自闭了
补了一个$n^4$的做法
然后似乎有一个很套路的可以变成$n^3$的做法
你在枚举$i,j,k$转移的时候，枚举中间那个点
然后枚举需要转移的那个点
如果你把其他所有点按到中间那个点极角排序一下
能转移过来的区间是连续一段
双指针优化一下就能$n^3$了

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#include<set>
#define LL long long
#define mp(x,y) make_pair(x,y)
#define pll pair<long long,long long>
#define pii pair<int,int>
using namespace std;
inline int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int stack[20];
inline void write(LL x)
{
	if(x<0){putchar('-');x=-x;}
    if(!x){putchar('0');return;}
    int top=0;
    while(x)stack[++top]=x%10,x/=10;
    while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);putchar(' ');}
inline void pr2(LL x){write(x);putchar('\n');}
const LL mod=998244353;
const int MAXN=205;
LL pow_mod(LL a,int b)
{
	LL ret=1;
	while(b)
	{
		if(b&1)ret=ret*a%mod;
		a=a*a%mod;b>>=1;
	}
	return ret;
}
void dl(LL &x,LL y){x-=y;if(x<0)x+=mod;}
void ad(LL &x,LL y){x+=y;if(x>=mod)x-=mod;}
struct pt{int x,y;}w[MAXN];
int n;
int multi(int p0,int p1,int p2)
{
	int u1=w[p0].x-w[p1].x;
	int u2=w[p0].x-w[p2].x;
	int v1=w[p0].y-w[p1].y;
	int v2=w[p0].y-w[p2].y;;
	return u1*v2-u2*v1;
}
bool ok(int u1,int u2)
{
	for(int i=1;i<u2;i++)if(i!=u1){if(!multi(u1,u2,i))return false;}
	return true;
}
int main()
{
	n=read();
	for(int i=1;i<=n;i++)w[i].x=read(),w[i].y=read();
	LL ans=(pow_mod(2LL,n)-1+mod)%mod;
	for(int i=1;i<=n;i++)for(int j=i+1;j<=n;j++)if(ok(i,j))
	{
		int cnt=2;
		for(int k=j+1;k<=n;k++)if(!multi(i,j,k))cnt++;
		LL sum=(pow_mod(2,cnt)-(cnt+1)%mod+mod)%mod;
		dl(ans,sum);
	}
	pr2((ans-n+mod)%mod);
	return 0;
}

没卡过的$n^4$

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#include<set>
#define LL long long
#define mp(x,y) make_pair(x,y)
#define pll pair<long long,long long>
#define pii pair<int,int>
using namespace std;
inline int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int stack[20];
inline void write(LL x)
{
	if(x<0){putchar('-');x=-x;}
    if(!x){putchar('0');return;}
    int top=0;
    while(x)stack[++top]=x%10,x/=10;
    while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);putchar(' ');}
inline void pr2(LL x){write(x);putchar('\n');}
const LL mod=998244353;
const int MAXN=205;
struct pt{int x,y;double pp;}w[MAXN],g[MAXN];int tot;
int f[MAXN][MAXN],n,pre[MAXN][MAXN],l1[MAXN];
int multi(pt p0,pt p1,pt p2)
{
	int ux,uy,vx,vy;
	ux=p0.x-p1.x;uy=p0.y-p1.y;
	vx=p0.x-p2.x;vy=p0.y-p2.y;
	return ux*vy-vx*uy;
}
bool cmp(pt u1,pt u2){return u1.pp!=u2.pp?u1.pp<u2.pp:u1.y<u2.y;}
bool check(pt u1,pt u2,pt u3,pt u4)
{
	int g1=multi(u1,u2,u3),g2=multi(u1,u2,u4);
	if((g1<=0&&g2<=0)||(g1>=0&&g2>=0))return true;
	return false;
}
bool ok(pt p0,pt p1,pt p2,pt p3)//p0~p1不能取 
{
	if(multi(p0,p1,p3)==0)return false;
	if(!check(p0,p1,p2,p3)||!check(p1,p2,p0,p3)||!check(p2,p0,p1,p3))return false;
	return true;
}
void ad(int &x,int y){x+=y;if(x>=mod)x-=mod;}
int temp[MAXN];
int ans;
bool in(pt u1,pt u2,pt u3)
{
	if(u2.x>u3.x)swap(u2,u3);
	else if(u2.x==u3.x&&u2.y>u3.y)swap(u2,u3);
	if(u1.x>=u2.x&&u1.x<=u3.x&&(u1.y>=u2.y&&u1.y<=u3.y||u1.y>=u3.y&&u1.y<=u2.y))return true;
	return false;
}
void solve(int x)
{
	tot=0;
	for(int i=1;i<=n;i++)if(w[i].y>w[x].y||(w[i].y==w[x].y&&w[i].x>w[x].x))g[++tot]=w[i];
	for(int i=1;i<=tot;i++)g[i].pp=atan2(g[i].y-w[x].y,g[i].x-w[x].x);
	sort(g+1,g+1+tot,cmp);
	memset(pre,0,sizeof(pre));
	memset(l1,0,sizeof(l1));
	for(int i=1;i<=tot;i++)
		for(int j=1;j<=n;j++)
			if(multi(w[x],g[i],w[j])==0&&in(w[j],w[x],g[i]))
				l1[i]++;
	for(int i=1;i<=tot;i++)for(int j=1;j<=tot;j++)if(i!=j)
	{
		for(int k=1;k<=n;k++)if(ok(w[x],g[i],g[j],w[k]))
			pre[i][j]++;
	}
	memset(f,0,sizeof(f));g[0]=w[x];
	for(int i=1;i<=tot;i++)
		for(int j=i+1;j<=tot;j++)
		{
			for(int k=0;k<i;k++)if(multi(g[k],g[i],g[j])>0)
			{
				int pt=pre[j][i];
				if(k==0)ad(f[i][j],temp[pt-1]);
				else ad(f[i][j],(LL)f[k][i]*temp[pt-1]%mod);
			}
			int pt=l1[j];
			ad(ans,(LL)f[i][j]*temp[pt-2]%mod);
		}
}
int main()
{
	temp[0]=1;for(int i=1;i<MAXN;i++)temp[i]=temp[i-1]*2%mod;
	n=read();
	for(int i=1;i<=n;i++)w[i].x=read(),w[i].y=read();
	for(int i=1;i<=n;i++)
		solve(i);
	pr2(ans);
	return 0;
}