BZOJ1766: [Ceoi2009]photo

给出平面上n<=100个点，求最少用多少面积≤S<=200000的矩形能全覆盖。

第一想法是区间DP，f[i,j]表示覆盖第i到j的点最少矩形数，$f[i,j]=∑f[i,x]+f[x+1,j]，x∈[i,j)$，但好在样例良心，这样的方案过不了样例。怎么把两个矩形相交的方案处理出来呢？

f[i,j,h]表示第i到第j个点（按横坐标排序）中，把高度>=h的点都覆盖了的最小矩形数，$f[i,j,k]=min(f[i,x,k]+f[x+1,j,k])$，若S/(Xj-Xi)>h，表明这个区间可以用两个矩形相交的方式放置，$f[i,j,k]<?=f[i,j,S/(Xj-Xi)+1]+1$

上式成立的前提是i,j两个点高度都>=h，在考虑范围内。如果不是，$f[i,j,k]=min(f[xl,xr,k])$,xl,xr表示在第i到第j个点中>=h的最左最右的点，xl>i,xr<j。

高度范围可以到S，因此先离散化。

 1 #include<stdio.h>
 2 #include<stdlib.h>
 3 #include<string.h>
 4 #include<algorithm>
 5 //#include<iostream>
 6 using namespace std;
 7 
 8 int n,S;
 9 #define maxn 111
10 #define maxs 200011
11 int h[maxs],rh[maxs],lh=0;
12 struct Point
13 {
14     int x,y;
15     bool operator < (const Point &b) const
16     {return x<b.x || (x==b.x && y<b.y);}
17 }tmp[maxn],a[maxn];int la=0;
18 int f[maxn][maxn][maxn];
19 const int inf=0x3f3f3f3f;
20 int main()
21 {
22     scanf("%d%d",&n,&S);
23     for (int i=1;i<=n;i++) scanf("%d%d",&tmp[i].x,&tmp[i].y);
24     sort(tmp+1,tmp+1+n);
25     for (int i=1;i<n;i++)
26         if (tmp[i].x!=tmp[i+1].x) a[++la].x=tmp[i].x,a[la].y=tmp[i].y;
27     a[++la].x=tmp[n].x,a[la].y=tmp[n].y;n=la;
28     for (int i=1;i<=n;i++) h[i]=a[i].y;
29     sort(h+1,h+1+n);
30     for (int i=1;i<=n;i++) if (h[i]!=h[i-1])
31     {
32         rh[h[i]]=++lh;
33         for (int j=h[i-1];j<h[i];j++) rh[j]=rh[h[i-1]];
34     }
35     for (int i=h[n];i<=200000;i++) rh[i]=rh[h[n]];
36     for (int i=1;i<=n;i++)
37     {
38         for (int j=rh[a[i].y]+1;j<=rh[h[n]];j++) f[i][i][j]=0;
39         for (int j=0;j<=rh[a[i].y];j++) f[i][i][j]=1;
40     }
41     for (int k=rh[h[n]];k>=0;k--)
42         for (int i=2;i<=n;i++)
43             for (int j=1;j<=n-i+1;j++)
44             {
45                 int xl=inf,xr=-inf,lid,rid;
46                 for (int x=j;x<=j+i-1;x++) if (a[x].y>=k)
47                 {
48                     if (a[x].x<xl)
49                     {
50                         xl=a[x].x;
51                         lid=x;
52                     }
53                     if (a[x].x>xr)
54                     {
55                         xr=a[x].x;
56                         rid=x;
57                     }
58                 }
59                 if (xl==a[j].x && xr==a[j+i-1].x)
60                 {
61                     f[j][j+i-1][k]=inf;
62                     int h2=rh[S/(a[j+i-1].x-a[j].x)]+1;
63                     if (h2>k) f[j][j+i-1][k]=min(f[j][j+i-1][k],f[j][j+i-1][h2]+1);
64                     for (int x=j;x<j+i-1;x++)
65                         f[j][j+i-1][k]=min(f[j][x][k]+f[x+1][j+i-1][k],f[j][j+i-1][k]);
66                 }
67                 else if (xl!=inf) f[j][j+i-1][k]=f[lid][rid][k];
68                 else f[j][j+i-1][k]=0;
69             }
70     printf("%d\n",f[1][n][0]);
71     return 0;
72 }

View Code

 

转载于:https://www.cnblogs.com/Blue233333/p/7304852.html