USACO 2016 February Contest, Gold解题报告

1.Circular Barn 

 http://www.usaco.org/index.php?page=viewproblem2&cpid=621

贪心

#include <cstdio>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
long long sum(long long v)
{
  return v*(v+)*(*v+)/;
}
int main()
{
      int N,c=;
      scanf("%d",&N);
      vector<long long> A(N);
      for(int i=;i<N;i++)
    {
        scanf("%d",&A[i]);
        c=max(0ll,c+A[i]-);
      }
      for(int i=;;i++)
    {
        if(c==)
        {
              rotate(A.begin(),A.begin()+i,A.begin()+N);
              break;
        }
        c=max(0ll,c+A[i]-);
      }
      long long res=;
      for (int i=;i<N;i++)
    {
        res+=sum(A[i]+c-)-sum(c-);
        c=max(0ll,c+A[i]-);
      }
      printf("%d",res);
}

2.Circular Barn Revisited

http://www.usaco.org/index.php?page=viewproblem2&cpid=622

首先，记录从i开始的j个房间的奶牛都从i从i号门进入所走的总距离。

然后进行dp，找出最小路程

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 100
#define MAXK 7
int n,r[MAXN+],k;
long long f[MAXN*+][MAXN+][MAXK+],ans=0x7fffffffffffffffll,dist[MAXN+][MAXN+];
void Read(int &x){
    char c;
    while(c=getchar(),c!=EOF)
        if(c>=''&&c<=''){
            x=c-'';
            while(c=getchar(),c>=''&&c<='')
                x=x*+c-'';
            ungetc(c,stdin);
            return;
        }
}
void read(){
    Read(n),Read(k);
    int i;
    for(i=;i<=n;i++)
        Read(r[i]);
}
void prepare(){
    int i,j;
    for(i=;i<=n;i++){
        dist[i][]=;
        for(j=;j<=n;j++)
            dist[i][j]=r[(i+j-)%n+]*(j-)+dist[i][j-];
    }
}
void dp(){
    int i,j,l,t=*n,p;
    memset(f,0x3f,sizeof f);
    for(i=;i<=t;i++)
        f[i][][]=;
    for(i=;i<=n;i++){
        for(j=;j<=n&&j<=i;j++)
            for(l=;l<=k;l++)
                for(p=;p<=j;p++)
                    f[i][j][l]=min(f[i][j][l],f[i-p][j-p][l-]+dist[i-p+][p]);
        ans=min(ans,f[i][n][k]);
    }
    for(i=n+;i<=t;i++){
        for(j=;j<=n;j++)
            for(l=;l<=k;l++)
                for(p=;p<=j;p++)
                    f[i][j][l]=min(f[i][j][l],f[i-p][j-p][l-]+dist[(i-p)%n+][p]);
        ans=min(ans,f[i][n][k]);
    }
}
int main()
{
    read();
    prepare();
    dp();
    printf("%lld\n",ans);
}

3.Fenced In

最小生成树

处理出连通块之间栅栏的长度然后做最小生成树即可。 但是这样会T，尽管算法时间复杂度为 O(n2log2n) 

优化：由于每一行或每一列中的相邻两个连通块之间的栅栏长度相等，我们可以考虑使用kruskal算法，仅仅对每两个栅栏之间的距离排序，然后然后再对这一行（列）的所有相邻连通块连边即可。 

#include<cstdio>
#include<algorithm>
#define MAXN 2000
using namespace std;
long long ans;
int n,a[MAXN+],b[MAXN+],A,B,m,tot,fa[(MAXN+)*(MAXN+)+];
void Read(int &x){
    char c;
    while(c=getchar(),c!=EOF)
        if(c>=''&&c<=''){
            x=c-'';
            while(c=getchar(),c>=''&&c<='')
                x=x*+c-'';
            ungetc(c,stdin);
            return;
        }
}
void read(){
    int i,t;
    Read(A),Read(B),Read(n),Read(m);
    for(i=;i<=n;i++)
        Read(a[i]);
    for(i=;i<=m;i++)
        Read(b[i]);
    sort(a+,a+n+);
    sort(b+,b+m+);
    a[++n]=A;
    b[++m]=B;
    for(i=n;i;i--)
        a[i]-=a[i-];
    for(i=m;i;i--)
        b[i]-=b[i-];
    sort(a+,a+n+);
    sort(b+,b+m+);
    t=n*m;
    for(i=;i<=t;i++)
        fa[i]=i;
}
int find(int x){
    return x==fa[x]?x:fa[x]=find(fa[x]);
}
void kruskal(){
    int i,j,k,l;
    i=j=;
    while(i<=n||j<=m){
        if(j>m||(i<=n&&a[i]<=b[j])){
            k=(i-)*m+;
            for(l=;l<m;l++,k++)
                if(find(k)!=find(k+))
                    fa[fa[k]]=fa[k+],ans+=a[i];
            ++i;
        }
        else{
            k=j;
            for(l=;l<n;l++,k+=m)
                if(find(k)!=find(k+m))
                    fa[fa[k]]=fa[k+m],ans+=b[j];
            ++j;
        }
    }
}
int main()
{
    read();
    kruskal();
    printf("%lld\n",ans);
}