NOIP2017提高组模拟赛 9 (总结)##
第一题 星星
天空中有N(1≤N≤400)颗星,每颗星有一个唯一的坐标(x,y),(1≤x,y ≤N)。请计算可以覆盖至少K(1≤K≤N)颗星的矩形的最小面积。矩形的边必须平行于X轴或Y轴,长度必须为正整数。星如果在矩形的边上,也认为它是属于矩形内的。
思路:枚举矩形左上角的x,以及在x轴上的长度(矩形的长),然后单调的往下扫。
#include<cstdio>
#include<algorithm>
#include<cmath>
#define imax(a,b) ((a>b)?(a):(b))
#define imin(a,b) ((a<b)?(a):(b))
typedef long long ll;
using namespace std;
int n,ne,ans;
int sum[410][410];
bool d[410][410];
int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
scanf("%d%d",&n,&ne);
ans=n*n;
for(int i=1;i<=n;i++)
{
int a,b;
scanf("%d%d",&a,&b); d[a][b]=1;
}
for(int i=1;i<=n;i++)
{
sum[i][0]=0;
for(int j=1;j<=n;j++) sum[i][j]=sum[i][j-1]+(d[i][j]?1:0);
}
for(int ll=1;ll<=n;ll++)
{
for(int i=1;i<=n-ll+1;i++)
{
int fr=1,la=1,ss=sum[fr][i+ll-1]-sum[fr][i-1];
while(la<=n)
{
while(ss<ne)
{
la++; ss+=sum[la][i+ll-1]-sum[la][i-1];
if(la>n) break;
}
if(la>n) break;
ans=imin(ans,(la-fr+1)*ll);
ss-=sum[fr][i+ll-1]-sum[fr][i-1]; fr++;
}
}
}
printf("%d\n",ans);
return 0;
}
第二题 战争
X国和Y国是死对头,X国有N个炮台, Y国有M个基地和K个发电站,炮台、基地、发电站所在的位置的坐标都是整数。Y国的每个发电站能对距离其L以内(包括L)的基地供电。X国的每个炮台都可以发射无限次,每次发射可以击毁一个基地或者一个发电站,消耗的能量为两者距离的平方,这里的距离是欧几里德距离。X国决定要摧毁Y国的所有基地,我们说Y国的某个基地被摧毁的条件是:基地本身被直接摧毁或者对其供电的所有发电站被击。
问X国摧毁Y国所有基地需要消耗的总能量最小是多少。
提示:点(X1,Y1)和点(X2,Y2)的欧几里德距离是:
dis = sqrt((X2-X1)*(X2-X1)+(Y2-Y1)*(Y2-Y1)).
经典的网络流01模型,直接构图,跑一次dinic就行了。
S向各个基地基地bi连边,代价为距离基地bi最近的炮台之间的距离。
各个电厂向T连边,代价为距离电厂ci最近的炮台之间的距离。
若电厂ci可以供电给基地dj,那么基地dj连向电厂ci,代价为∞。
图的S到T的最小割即为answer。
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define imax(a,b) ((a>b)?(a):(b))
#define imin(a,b) ((a<b)?(a):(b))
typedef long long ll;
using namespace std;
const int N=50;
const int M=15000;
const int oo=1e9;
struct data { int x,y; } ad[N+10],bd[N+10],cd[N+10];
int to[M],h[M],ne[M],val[M],tt,S,T;
int vis[M],lev[M],q[M];
int an,bn,cn,L,bfstime;
void read(int n,data* A)
{
for(int i=1;i<=n;i++) scanf("%d",&A[i].x);
for(int i=1;i<=n;i++) scanf("%d",&A[i].y);
}
void add(int a,int b,int c)
{
to[++tt]=b; ne[tt]=h[a]; val[tt]=c; h[a]=tt;
to[++tt]=a; ne[tt]=h[b]; val[tt]=0; h[b]=tt;
}
int dis(data A,data B) { return ((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y)); }
void init()
{
scanf("%d%d%d%d",&an,&bn,&cn,&L);
read(an,ad); read(bn,bd); read(cn,cd);
tt=1; S=110; T=111;
for(int i=1;i<=bn;i++)
{
int dd=oo;
for(int j=1;j<=an;j++) dd=imin(dd,dis(ad[j],bd[i]));
add(S,i,dd);
}
for(int i=1;i<=cn;i++)
{
int dd=oo;
for(int j=1;j<=an;j++) dd=imin(dd,dis(ad[j],cd[i]));
add(50+i,T,dd);
}
int ml=L*L;
for(int i=1;i<=bn;i++)
{
for(int j=1;j<=cn;j++)
if(dis(bd[i],cd[j])<=ml) add(i,50+j,oo);
}
}
bool bfs()
{
lev[S]=1; vis[S]=++bfstime;
int he=1,ta=1; q[1]=S;
while(he<=ta)
{
int u=q[he++];
for(int p=h[u];p;p=ne[p])
{
int v=to[p];
if(vis[v]==bfstime || !val[p]) continue;
q[++ta]=v; vis[v]=bfstime; lev[v]=lev[u]+1;
}
}
return (vis[T]==bfstime);
}
int dfs(int now,int maxf)
{
int ret=0,t=0;
if(!maxf || now==T) return maxf;
for(int p=h[now];p;p=ne[p])
{
int v=to[p];
if(lev[v]!=lev[now]+1 || !val[p]) continue;
t=dfs(v,imin(maxf,val[p]));
ret+=t; maxf-=t;
val[p]-=t; val[p^1]+=t;
}
if(!maxf) lev[now]=-1;
return ret;
}
int dinic()
{
int ret=0;
while(bfs()) ret+=dfs(S,oo);
return ret;
}
int main()
{
freopen("b.in","r",stdin);
freopen("b.out","w",stdout);
init();
printf("%d\n",dinic());
return 0;
}
第三题 染色树
一棵共含有X个结点的树,结点编号1至X,根结点编号是1。有Y种不同的颜色,颜色编号从1至Y。
现在给每个结点都染上一种颜色,整颗树染色后满足:
1、对于编号是i的颜色,整颗树当中,至少有一个结点被染成了颜色i。
2、根结点必须被染成1号颜色,而且整颗树当中,恰好要有Z个结点被染成1号颜色。
染色过程结束后,现在要计算染色的总代价,总代价等于每一条边的代价之和,那么怎么计算一条边的代价呢?
不妨设结点a与结点b之间有一条边,该边的权重是W。
1、如果结点a和结点b最终被染成不同的颜色,那么a与b之间的那条边的代价等于0。
2、如果结点a和结点b最终被染成相同的颜色,那么a与b之间的那条边的代价等于W。
现在的问题是:在满足上述要求的前提下,染色树最小的总代价是多少?如果无法完成染色的任务,输出-1。
显然的,这是一道树形DP题(DP做的有点多-_-||),X,Y,Z好像都很大。直接做无疑是超时的。
那么想考虑一下,answer=-1的情况。每个节点都可以是任意颜色(只是会有代价而已),唯一不满足的情况就是Z>X-Y+1(极限情况)。
这是一棵树(这很重要),那么只需要两种颜色就可以使整棵树的任意x和x的父亲的颜色不同。
又多了一个颜色1恰好为Z的限制条件,那么也只需要3种颜色就行了(多余的颜色是没任何作用的)。
时间复杂度:O(9N³)
#include<cstdio>
#include<algorithm>
#include<cmath>
#define imax(a,b) ((a>b)?(a):(b))
#define imin(a,b) ((a<b)?(a):(b))
typedef long long ll;
using namespace std;
const int oo=1e9;
const int M=700;
int X,Y,Z,ans;
int to[M],ne[M],h[M],val[M],tt;
int f[302][4][302];
void add(int a,int b,int c)
{
to[++tt]=b; ne[tt]=h[a]; h[a]=tt; val[tt]=c;
to[++tt]=a; ne[tt]=h[b]; h[b]=tt; val[tt]=c;
}
void dfs(int x,int fa)
{
for(int i=1;i<=Y;i++)
{
for(int j=0;j<=Z;j++) f[x][i][j]=oo;
if(i==1) f[x][i][1]=0; else f[x][i][0]=0;
}
for(int p=h[x];p;p=ne[p])
{
int v=to[p];
if(v==fa) continue;
dfs(v,x);
for(int i=1;i<=Y;i++)
for(int j=Z;j>=0;j--)
{
int bes=oo;
for(int g=1;g<=Y;g++) bes=imin(bes,f[v][g][0]+(i==g?val[p]:0));
f[x][i][j]=imin(oo,f[x][i][j]+bes);
for(int k=1;k<=j;k++)
{
int yu=oo;
for(int g=1;g<=Y;g++) yu=imin(yu,f[v][g][k]+(i==g?val[p]:0));
f[x][i][j]=imin(f[x][i][j],f[x][i][j-k]+yu);
}
}
}
}
int main()
{
freopen("c.in","r",stdin);
freopen("c.out","w",stdout);
scanf("%d%d%d",&X,&Y,&Z); tt=1;
if(Z>X-Y+1) printf("-1\n"); else
{
Y=imin(Y,3);
for(int i=1;i<X;i++) { int a,b,c; scanf("%d%d%d",&a,&b,&c); add(a,b,c); }
dfs(1,1);
printf("%d\n",f[1][1][Z]);
}
return 0;
}
终于拿了一次AK了(虽然这次的题目有点水……)。