bzoj 3778: 共鳴【計算幾何+dp】
阿新 • • 發佈:2019-05-03
con b- 三角形面積 cmp getchar space for cpp std
枚舉起點,然後設f[i][j]為上凸殼上一個點是i當前點是j的最大面積,g是下凸殼,然後合並的時候枚舉結束點t合並上下凸殼即可
這樣的好處是每次轉移都是往凸多邊形裏加一個三角形(s,i,j),所以判斷轉移合法只要預處理出所有三角形是否合法即可,同時預處理出三角形面積,轉移就是f[j][k]=max(f[j][k],f[i][j]+c[s][j][k]);
#include<iostream> #include<cstdio> #include<algorithm> #include<cmath> using namespace std; const int N=105; int n,m; double c[N][N][N],f[N][N],g[N][N],ans; bool v[N][N][N]; struct dian { double x,y; dian(double X=0,double Y=0) { x=X,y=Y; } dian operator + (const dian &a) const { return dian(x+a.x,y+a.y); } dian operator - (const dian &a) const { return dian(x-a.x,y-a.y); } }a[N],b[N]; bool cmp(const dian &a,const dian &b) { return a.x<b.x||(a.x==b.x&&a.y<b.y); } int read() { int r=0,f=1; char p=getchar(); while(p>'9'||p<'0') { if(p=='-') f=-1; p=getchar(); } while(p>='0'&&p<='9') { r=r*10+p-48; p=getchar(); } return r*f; } double cj(dian a,dian b) { return a.x*b.y-a.y*b.x; } bool ok(dian x,dian y,dian z) { for(int i=1; i<=m; i++) { if(cj(b[i]-x,y-x)>=0&&cj(b[i]-y,z-y)>=0&&cj(b[i]-z,x-z)>=0) return 0; if(cj(b[i]-x,y-x)<=0&&cj(b[i]-y,z-y)<=0&&cj(b[i]-z,x-z)<=0) return 0; } return 1; } double mj(dian a,dian b,dian c) { return abs(cj(b-a,c-a))/2; } int main() { n=read(),m=read(); for(int i=1;i<=n;i++) a[i].x=read(),a[i].y=read(); for(int i=1;i<=m;i++) b[i].x=read(),b[i].y=read(); sort(a+1,a+1+n,cmp); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) for(int k=1;k<=n;k++) v[i][j][k]=ok(a[i],a[j],a[k]),c[i][j][k]=mj(a[i],a[j],a[k]); for(int s=1;s<=n;s++) { for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) f[i][j]=g[i][j]=-1e9; for(int i=s+1;i<=n;i++) f[s][i]=g[s][i]=0; for(int i=s;i<=n;i++) for(int j=i+1;j<=n;j++) for(int k=j+1;k<=n;k++) if(v[s][j][k]&&cj(a[i]-a[j],a[k]-a[j])>=0) f[j][k]=max(f[j][k],f[i][j]+c[s][j][k]); for(int i=s;i<=n;i++) for(int j=i+1;j<=n;j++) for(int k=j+1;k<=n;k++) if(v[s][j][k]&&cj(a[i]-a[j],a[k]-a[j])<=0) g[j][k]=max(g[j][k],g[i][j]+c[s][j][k]); for(int t=s+1;t<=n;t++) for(int i=s;i<=t;i++) for(int j=s;j<=t;j++) ans=max(ans,f[i][t]+g[j][t]); } printf("%.2f\n",ans); return 0; }
bzoj 3778: 共鳴【計算幾何+dp】