降雨量java_P4273 [NOI2004] 降雨量

把马路竖起来，横轴是时间

那么显然扫过的面积就是遮挡的水量

之后就是计算几何的事了

#include

#include

#include

#include

#include

using namespace std;

int n,k;

int xx[505],yy[505];

int x,y;

struct b{

int f;

int t;

double l;

} ed[2500];

int fa[100001];

int bl[2500];

struct re{

int lx,rx,uy,dy;

}rec[2500];

re tem;

re recc[5];

int find(int x){

return fa[x]==x?x:fa[x]=find(fa[x]);

}

int f(int x){

return x*x;

}

int di(int x,int y){

return f(xx[x]-xx[y])+f(yy[x]-yy[y]);

}

bool cmp(b x,b y){

return x.l

}

int p;

void add(int f,int to,int ne){

ed[++p].f=f;

ed[p].t=to;

ed[p].l=ne;

return ;

}

int cnt;

void com(int x,int y){

recc[y].lx=min(recc[y].lx,rec[x].lx);

recc[y].uy=min(recc[y].uy,rec[x].uy);

recc[y].rx=max(recc[y].rx,rec[x].rx);

recc[y].dy=max(recc[y].dy,rec[x].dy);

}

void can(int x,int y,re tem){

//cnt--;

recc[y]=tem;

}

int Aimee;

int znx=0x3fffff;

int ss(){

Aimee=0;

for(int i=1;i<=cnt;++i){

Aimee+=(recc[i].rx-recc[i].lx)*(recc[i].dy-recc[i].uy);

}

return Aimee;

}

void endd(){

Aimee=0;

for(int i=1;i<=cnt;++i){

Aimee+=(recc[i].rx-recc[i].lx)*(recc[i].dy-recc[i].uy);

}

znx=min(znx,Aimee);

}

bool ch(int i,int xx,int yy){

return (xx<=recc[i].rx&&xx>=recc[i].lx&&recc[i].dy>=yy&&recc[i].uy<=yy);

}

bool che(int x,int xx,int yy){

for(int i=1;i<=cnt;++i){

if(i==x)

continue;

if(ch(x,recc[i].lx,recc[i].dy)||

ch(x,recc[i].lx,recc[i].uy)||

ch(x,recc[i].rx,recc[i].dy)||

ch(x,recc[i].rx,recc[i].uy))

return 0;

}

return 1;

}

void kr(int i,int s){

if(s>znx){

return ;

}

if(i==n+1){

if(cnt==k)

znx=s;

return ;

}

if(cnt

re tem=recc[cnt+1];

cnt++;

recc[cnt]=rec[i];

if(che(cnt,rec[i].lx,rec[i].dy))

kr(i+1,ss());

recc[cnt]=tem;

cnt--;

}

for(int j=1;j<=cnt;++j){

re tem=recc[j];

com(i,j);

if(che(j,rec[i].lx,rec[i].dy))

kr(i+1,ss());

can(i,j,tem);

}

}

int main(){

scanf("%d%d",&n,&k);

for(int i=1;i<=n;++i){

scanf("%d%d",&x,&y);

rec[i].lx=rec[i].rx=x;

rec[i].dy=rec[i].uy=y;

}

kr(1,0);

cout<

return 0;

}

```

```

`那么显然扫过的面积就是遮挡的水量

之后就是计算几何的事了

````cpp

#include

#include

#include

#include

#include

using namespace std;

int n,k;

int xx[505],yy[505];

int x,y;

struct b{

int f;

int t;

double l;

} ed[2500];

int fa[100001];

int bl[2500];

struct re{

int lx,rx,uy,dy;

}rec[2500];

re tem;

re recc[5];

int find(int x){

return fa[x]==x?x:fa[x]=find(fa[x]);

}

int f(int x){

return x*x;

}

int di(int x,int y){

return f(xx[x]-xx[y])+f(yy[x]-yy[y]);

}

bool cmp(b x,b y){

return x.l

}

int p;

void add(int f,int to,int ne){

ed[++p].f=f;

ed[p].t=to;

ed[p].l=ne;

return ;

}

int cnt;

void com(int x,int y){

recc[y].lx=min(recc[y].lx,rec[x].lx);

recc[y].uy=min(recc[y].uy,rec[x].uy);

recc[y].rx=max(recc[y].rx,rec[x].rx);

recc[y].dy=max(recc[y].dy,rec[x].dy);

}

void can(int x,int y,re tem){

//cnt--;

recc[y]=tem;

}

int Aimee;

int znx=0x3fffff;

int ss(){

Aimee=0;

for(int i=1;i<=cnt;++i){

Aimee+=(recc[i].rx-recc[i].lx)*(recc[i].dy-recc[i].uy);

}

return Aimee;

}

void endd(){

Aimee=0;

for(int i=1;i<=cnt;++i){

Aimee+=(recc[i].rx-recc[i].lx)*(recc[i].dy-recc[i].uy);

}

znx=min(znx,Aimee);

}

bool ch(int i,int xx,int yy){

return (xx<=recc[i].rx&&xx>=recc[i].lx&&recc[i].dy>=yy&&recc[i].uy<=yy);

}

bool che(int x,int xx,int yy){

for(int i=1;i<=cnt;++i){

if(i==x)

continue;

if(ch(x,recc[i].lx,recc[i].dy)||

ch(x,recc[i].lx,recc[i].uy)||

ch(x,recc[i].rx,recc[i].dy)||

ch(x,recc[i].rx,recc[i].uy))

return 0;

}

return 1;

}

void kr(int i,int s){

if(s>znx){

return ;

}

if(i==n+1){

if(cnt==k)

znx=s;

return ;

}

if(cnt

re tem=recc[cnt+1];

cnt++;

recc[cnt]=rec[i];

if(che(cnt,rec[i].lx,rec[i].dy))

kr(i+1,ss());

recc[cnt]=tem;

cnt--;

}

for(int j=1;j<=cnt;++j){

re tem=recc[j];

com(i,j);

if(che(j,rec[i].lx,rec[i].dy))

kr(i+1,ss());

can(i,j,tem);

}

}

int main(){

scanf("%d%d",&n,&k);

for(int i=1;i<=n;++i){

scanf("%d%d",&x,&y);

rec[i].lx=rec[i].rx=x;

rec[i].dy=rec[i].uy=y;

}

kr(1,0);

cout<

return 0;

}