正文

Gym -102007 ：Benelux Algorithm Programming Contest (BAPC 18) (寒假自训第5场)

IP云
IP云V管理员/4 阅读
09/05
提示：扫一扫查出行
【扫一扫了解最新限行尾号】
复制提示

A .A Prize No One Can Win

题意：给定N，S，你要从N个数中选最多是数，使得任意两个之和不大于S。

思路：排序，然后贪心的选即可。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
ll a[maxn];
int main()
{
int N,ans; ll M;
scanf("%d%lld",&N,&M);
rep(i,,N) scanf("%lld",&a[i]);
sort(a+,a+N+);
ans=;
rep(i,,N){
if(a[i]+a[i-]<=M) ans=i;
}
printf("%d\n",ans);
return ;
}

B .Birthday Boy

题意：给出N个人的生日，让你选择一天，使得这一天的前一个生日距离它最远，如果有多个一样的，有点选择10月27之后的

思路：模拟即可。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int Y[]={,,,,,,,,,,,,};
int vis[],tot;
int mp[][],C[][],A[],B[];
int main()
{
int N,y,d,ansx=,ansy=,ans=-;char c[];
scanf("%d",&N);
rep(i,,N){
scanf("%s%d-%d",c+,&y,&d);
mp[y][d]=;
}
rep(i,,)
rep(j,,Y[i]){
tot++; A[tot]=i; B[tot]=j; C[i][j]=tot;
if(mp[i][j]) vis[tot]=;
}
rep(i,,tot) vis[i+tot]=vis[i],A[i+tot]=A[i],B[i+tot]=B[i];
rep(i,tot,tot+tot-) {
if(vis[i]) continue;
int j=i; while(j->=i-tot+&&!vis[j-]) j--;
if(i-j>ans){
ans=i-j,ansx=A[i],ansy=B[i];
}
else if(i-j==ans&&C[ansx][ansy]<=C[][]&&C[A[i]][B[i]]>C[][]) ansx=A[i],ansy=B[i];
}
printf("%02d-%02d\n",ansx,ansy);
return ;
}

C .Cardboard Container

题意: 给定N个1*1*1的盒子,让你用最小的表面积把它包起来.

思路: 枚举长和宽和高(因子)即可.

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=;
const int maxn=+;
const int inf=;
int main()
{
int v;
cin>>v;
ll ans=inf;
for(int i=;i<=v;i++){
if(v%i==){
int l=v/i;
for(int j=;j<=l;j++){
if(l%j==){
int k=l/j;
ans=min(ans,2LL*(i*j+k*j+k*i));
}
}
}
}
cout<<ans<<endl;
return ;
}

E .Entirely Unsorted Sequences

https://blog.csdn.net/ccsu_cat/article/details/86754931

G .Game Night

模拟几种排列即可

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int sum[maxn][],N; char c[maxn];
int get(int A,int B,int C)
{
int res=;
rep(i,,N){
int X=sum[i+sum[N][A]-][A]-sum[i-][A];
int Y=sum[i+sum[N][A]+sum[N][B]-][B]-sum[i+sum[N][A]-][B];
int Z=sum[i+N-][C]-sum[i+sum[N][A]+sum[N][B]-][C];
res=max(res,X+Y+Z);
//cout<<A<<" "<<B<<" "<<C<<" : "<<i<<" "<<res<<endl;
}
return N-res;
}
int main()
{
scanf("%d%s",&N,c+);
int ans=N;
rep(i,,N) c[N+i]=c[i];
rep(i,,N+N){
rep(j,,) sum[i][j]=sum[i-][j];
if(c[i]=='A') sum[i][]++;
if(c[i]=='B') sum[i][]++;
if(c[i]=='C') sum[i][]++;
}
rep(i,,)
rep(j,,)
rep(k,,)
if(i!=j&&i!=k&&j!=k)
ans=min(ans,get(i,j,k));
printf("%d\n",ans);
return ;
}

I. In Case of an Invasion, Please...

题意：给定N城市每个城市人数是Pi，M双向路，其中S个点(S<=10)是避难所，避难所容量是Ci，(N<1e5; Ci,Pi<1e9)最小化最大的转移时间，使得所有人都转移到避难所。

思路：二分+最大流。由于图比较大，此题需要一些优化才能过。

由于上限是1e14； 所以二分次数是50次； 我们可以把距离离散化上线变为了1e6，二分次数是20；这是一个优化。

第二个优化，对于二分到的mid，如果一个城市可以去的避难所容量之和<P[i]，那么显然不用跑网络流也知道不合法。

还有就是sap跑二分图比较慢，要用dinic。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e5+;
ll inf=1ll<<;
struct Edge
{
int from,to;
ll cap,flow;
};
struct Dinic
{
int n,m,s,t;
vector<Edge>edges;
vector<int>G[maxn];
bool vis[maxn];
int d[maxn],cur[maxn];
void init(int n)
{
this->n=n;
for(int i=;i<n;i++)
G[i].clear();
edges.clear();
}
void AddEdge(int from,int to,ll cap)
{
edges.push_back((Edge){from,to,cap,});
edges.push_back((Edge){to,from,,});
m=edges.size();
G[from].push_back(m-);
G[to].push_back(m-);
}
bool bfs()
{
memset(vis,,sizeof(vis));
queue<int>Q;
Q.push(s);
d[s]=;
vis[s]=;
while(!Q.empty())
{
int x=Q.front();Q.pop();
for(int i=;i<G[x].size();i++)
{
Edge& e=edges[G[x][i]];
if(!vis[e.to]&&e.cap>e.flow)
{
vis[e.to]=;
d[e.to]=d[x]+;
Q.push(e.to);
}
}
}
return vis[t];
}
ll dfs(int x,ll a)
{
if(x==t||a==)return a;
ll flow=,f;
for(int& i=cur[x];i<G[x].size();i++)
{
Edge& e=edges[G[x][i]];
if(d[x]+==d[e.to]&&(f=dfs(e.to,min(a,e.cap-e.flow)))>)
{
e.flow+=f;
edges[G[x][i]^].flow-=f;
flow+=f;
a-=f;
if(a==)break;
}
}
return flow;
}
ll Maxflow(int s,int t)
{
this->s=s,this->t=t;
ll flow=;
while(bfs())
{
memset(cur,,sizeof(cur));
flow+=dfs(s,inf);
}
return flow;
}
}solve;
ll d[][maxn],sum,vis[maxn];
vector<int>G[maxn],dis[maxn];
struct node
{
int u;
ll w;
node(int a,ll b)
{
u=a,w=b;
}
bool operator<(const node&t)const
{
return w>t.w;
}
};
priority_queue<node>q;
void dij(int k,int s,int n)
{
for(int i=;i<=n;i++)d[k][i]=inf,vis[i]=;
d[k][s]=;
q.push(node(s,));
while(!q.empty())
{
node e=q.top();q.pop();
int u=e.u; vis[u]=;
for(int i=;i<G[u].size();i++)
{
int v=G[u][i];
if(d[k][v]>d[k][u]+dis[u][i])
{
d[k][v]=d[k][u]+dis[u][i];
if(!vis[v]) vis[v]=,q.push(node(v,d[k][v]));
}
}
}
}
int p[maxn],s[maxn],c[maxn],tot; ll B[maxn*],fcy[maxn];
int main()
{
int n,m,k,u,v,w;
scanf("%d%d%d",&n,&m,&k);
int S=,T=n+k+;
for(int i=;i<=n;i++)
{
scanf("%d",&p[i]);
sum+=p[i];
}
for(int i=;i<=m;i++)
{
scanf("%d%d%d",&u,&v,&w);
G[u].push_back(v);
dis[u].push_back(w);
G[v].push_back(u);
dis[v].push_back(w);
}
for(int i=;i<=k;i++)
{
scanf("%d%d",&s[i],&c[i]);
dij(i,s[i],n);
}
for(int i=;i<=k;i++)
for(int j=;j<=n;j++) B[++tot]=d[i][j];
sort(B+,B+tot+);
tot=unique(B+,B+tot+)-(B+);
ll l=,r=tot,mid,ans;
while(l<=r)
{
mid=(l+r)/;
solve.init(T+);
for(int i=;i<=n;i++) fcy[i]=;
for(int i=;i<=n;i++)
for(int j=;j<=k;j++)
if(d[j][i]<=B[mid])
solve.AddEdge(i,n+j,inf),fcy[i]+=c[j];
bool F=true;
for(int i=;i<=n;i++)
if(fcy[i]<p[i]){ F=false; break;}
if(!F) {l=mid+;continue;}
for(int i=;i<=n;i++)
solve.AddEdge(S,i,p[i]);
for(int i=;i<=k;i++)
solve.AddEdge(n+i,T,c[i]);
ll res=solve.Maxflow(S,T);
if(res==sum) r=mid-,ans=B[mid];
else l=mid+;
}
printf("%lld\n",ans);
}

J .Janitor Troubles

题意：给定ABCD四条边，让你组成一个面积最大的四边形，保证有解。

思路：比赛时三分做的，1A。但是赛后拿去做hduwa掉了，还是没弄清楚。

三分：枚举相邻边的组合形式，然后对对角线进行三分。 结论是当三分到对角线的对面两个角都是直角时，答案最大，此时正好有外接圆。

定理：根据上面结论，相当于四边形是圆内接四边形，我们可以至直接用海伦公式的变形；

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
double X[],ans;
double S(double x,double y,double x1,double y1,double z)
{
double p=(x+y+z)/,p1=(x1+y1+z)/;
return sqrt(p*(p-x)*(p-y)*(p-z))+sqrt(p1*(p1-x1)*(p1-y1)*(p1-z));;
}
double get(int i,int j,int k,int p)
{
double Mn=max(X[j]-X[i],X[p]-X[k]),Mx=min(X[j]+X[i],X[p]+X[k]);
if(Mn>Mx) return 0.0;
int T=; double L=Mn,R=Mx,res=0.0;
while(T--){
double Mid1=L+(R-L)/,Mid2=R-(R-L)/;
double F1=S(X[i],X[j],X[k],X[p],Mid1);
double F2=S(X[i],X[j],X[k],X[p],Mid2);
if(F1>=F2) res=max(F1,res),R=Mid2;
else res=max(F2,res),L=Mid1;
}
return res;
}
int main()
{
rep(i,,) scanf("%lf",&X[i]);
sort(X+,X++);
rep(i,,)
rep(j,i+,)
rep(k,,)
rep(p,k+,)
if(i!=k&&i!=p&&j!=k&&j!=p){
ans=max(ans,get(i,j,k,p));
}
printf("%.10lf\n",ans);
return ;
}

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
double X[],ans;
double get()
{
double res=0.0,p=(X[]+X[]+X[]+X[])/;
res=sqrt((p-X[])*(p-X[])*(p-X[])*(p-X[]));
return res;
}
int main()
{
rep(i,,) scanf("%lf",&X[i]);
ans=get();
printf("%.10lf\n",ans);
return ;
}

K .Kingpin Escape

题意：给定一棵有根树，现在让你加最少的边，使得无论原图上哪条边被砍掉，从任意点出发都可以回到根节点。

思路：首先把保证至少两个点与根相邻，所以如果根只有一条边与它相邻，它需要加边。 其次，所有叶子节点需要加边，因为他和父亲被砍断后就GG了。

我们按照一定顺序把这些点连边即可。 假设X个点需要连边，那么最少需要加(X+1)/2条边。然后两两连边即可，但是注意至少一条新加的边连通了根的两个子树。

所以我们按照DFS序，然后q[i]+q[X/2+i]连边，这样可以保证，并不是所有新加的边都在同一个子树（这样可能合法，因为砍掉根与儿子后，子树无法到达根）里。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int Laxt[maxn],Next[maxn],To[maxn];
int sz[maxn],cnt,q[maxn],tot;
void add(int u,int v){
Next[++cnt]=Laxt[u]; Laxt[u]=cnt; To[cnt]=v;
}
void dfs(int u,int f)
{
if(sz[u]==) q[++tot]=u-;
for(int i=Laxt[u];i;i=Next[i])
if(To[i]!=f) dfs(To[i],u);
}
int main()
{
int N,M,u,v;
scanf("%d%d",&N,&M); M++;
rep(i,,N-){
scanf("%d%d",&u,&v);
u++; v++; sz[u]++; sz[v]++;
add(u,v); add(v,u);
}
dfs(M,);
int ans=(tot+)/;
printf("%d\n",ans);
rep(i,,tot/) printf("%d %d\n",q[i],q[i+tot/]);
if(tot&) printf("%d %d\n",q[],q[tot]);
return ;
}