bzoj 2197： [Usaco2011 Mar]Tree Decoration

　　树形dp。。f[i]表示处理完以i为根的子树的最小时间。

　　因为一个点上可以挂无数个，所以在点i上挂东西的单位花费就是i所在子树里的最小单位花费。。

　　所以每次求f[i]只要使子树里的数量都满足要求就好了。。i的祖先还要更多的话随时可以选某个节点多挂一些。。

　　f[i]=sum{f[j]}+mincost[i]*max(need[i]-sum{need[j]},0)。。（j是i的儿子，mincost[i]表示子树i里的最小花费，need[i]表示i这颗子树需要的最少数量）

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #define ll long long
  5.  using namespace std;
  6.  const int maxn=;
  7.  struct zs{
  8.      int too,pre;
  9.  }e[maxn];
  10.  int last[maxn],mn[maxn],tot;
  11.  int i,j,k,n,m,a;
  12.  ll f[maxn],need[maxn];
  13.  int ra;char rx;
  14.  inline int read(){
  15.      rx=getchar();ra=;
  16.      while(rx<''||rx>'')rx=getchar();
  17.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  18.  }
  19.  void dfs(int x){
  20.      int i,to;ll sum=;
  21.      for(i=last[x],to=e[i].too;i;mn[x]=mn[x]>mn[to]?mn[to]:mn[x],sum+=need[to],f[x]+=f[to],i=e[i].pre,to=e[i].too)dfs(to);
  22.      if(sum>=need[x])need[x]=sum;else f[x]+=(need[x]-sum)*(ll)mn[x];
  23.  }
  24.  int main(){
  25.      n=read();for(rx=getchar();rx!='-';rx=getchar());rx=getchar();
  26.      need[]=read();mn[]=read();
  27.      for(i=;i<=n;i++)a=read(),e[i].pre=last[a],last[a]=e[i].too=i,need[i]=read(),mn[i]=read();
  28.      dfs();
  29.      printf("%lld\n",f[]);
  30.      return ;
  31.  }

bzoj 2097： [Usaco2010 Dec]Exercise 奶牛健美操

　　没思路跑去看题解系列。。。二分答案。。

　　假设二分出来的答案是mid，那么对于每颗子树，设根节点为i，那么删边后子树的直径不能超过mid。

　　子树内的路径分两种情况，一种是经过i的，另一种是不经过i的。不经过i的显然可以递归i的子节点去删边，删完边后，如果i的任意两个儿子可以通过i凑成一条超过mid的路径的话，就贪心地把i和对长度贡献较大的孩子之间的连边删掉。

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  #include<cstdlib>
  6.  using namespace std;
  7.  const int maxn=;
  8.  struct child{
  9.      int x,len;
  10.  };
  11.  struct zs{
  12.      int too,pre;
  13.  }e1[maxn<<],e[maxn];
  14.  int tot,t1;
  15.  int l1[maxn],last[maxn],mxdep[maxn];
  16.  int i,j,mx,n,m,cnum,a,b,l,r,mid;
  17.  child ch[maxn];
  18.  int ra;char rx;
  19.  inline int read(){
  20.      rx=getchar();ra=;
  21.      while(rx<''||rx>'')rx=getchar();
  22.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  23.  }
  24.  inline void ins1(int a,int b){
  25.      e1[++t1].too=b;e1[t1].pre=l1[a];l1[a]=t1;
  26.      e1[++t1].too=a;e1[t1].pre=l1[b];l1[b]=t1;
  27.  }
  28.  inline void insert(int a,int b){e[++tot].too=b;e[tot].pre=last[a];last[a]=tot;}
  29.  void dfs1(int x,int fa){for(int i=l1[x];i;i=e1[i].pre)if(e1[i].too!=fa)insert(x,e1[i].too),dfs1(e1[i].too,x);}
  30.  bool cmp(child a,child b){return a.len<b.len;}
  31.  int dfs(int x,int val){
  32.      int i,ans=,to;
  33.      for(i=last[x];i&&ans<=mx;i=e[i].pre)
  34.          ans+=dfs(e[i].too,val);
  35.      if(ans>mx)return mx+;
  36.      cnum=;
  37.      for(i=last[x],to=e[i].too;i;i=e[i].pre,to=e[i].too)
  38.          ch[++cnum].x=to,ch[cnum].len=mxdep[to]+;
  39.      sort(ch+,ch++cnum,cmp);ch[].len=;
  40.      for(i=cnum;i&&ch[i].len+ch[i-].len>val;i--)ans++;
  41.      mxdep[x]=ch[i].len;
  42.      return ans;
  43.  }
  44.  int main(){
  45.      srand();
  46.      n=read();mx=read();if(mx==n-){puts("");return ;}if(mx==n-){puts("");return ;}
  47.      for(i=;i<n;i++)a=read(),b=read(),ins1(a,b);
  48.      int rt=rand()%n+;
  49.      dfs1(rt,);
  50.      l=;r=n-;
  51.      while(l<r){
  52.          mid=(l+r)>>;
  53.          if(dfs(rt,mid)<=mx)r=mid;else l=mid+;
  54.      }
  55.      printf("%d\n",l);
  56.      return ;
  57.  }

bzoj 1694： [Usaco2007 Demo]Grazing on the Run

　　又是英文题面+双倍经验= =同bzoj1742。。。。双倍经验==>双倍#1。。233

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<cstdlib>
  5.  #include<algorithm>
  6.  using namespace std;
  7.  const int maxn=;
  8.  int f[maxn][];
  9.  int len,s,mxj,n;
  10.  int a[maxn],b[maxn];
  11.  int ra;char rx;
  12.  inline int read(){
  13.      rx=getchar();ra=;
  14.      while(rx<''||rx>'')rx=getchar();
  15.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  16.  }
  17.  int main(){
  18.      int i,j,mxj,tmp;
  19.      n=read();s=read();
  20.      for(i=;i<=n;i++)a[i]=read();sort(a+,a++n);
  21.      for(i=;i<=n;b[i]=a[i+]-a[i],i++)f[i][]=f[i][]=abs(s-a[i])*n;
  22.      for(i=,mxj=n-;i<n;mxj--,i++)for(j=;j<=mxj;j++)
  23.      tmp=f[j][],f[j][]=min(f[j+][]+mxj*b[j],f[j+][]+mxj*(a[j+i]-a[j])),
  24.                  f[j][]=min(tmp+mxj*(a[j+i]-a[j]),f[j][]+mxj*b[j+i-]);
  25.      printf("%d\n",min(f[][],f[][]));
  26.      return ;
  27.  }

bzoj 2501： [usaco2010 Oct]Soda Machine

　　区间加后查询最大值。。写了离散化后二分。。其实离散化姿势正确的话是不用二分找出区间左右端点位置的TAT。

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  using namespace std;
  6.  const int maxn=;
  7.  int mp[maxn<<],l[maxn],r[maxn],sum[maxn<<];
  8.  int i,j,ans,n,size,nowsum;
  9.  int ra;char rx;
  10.  inline int read(){
  11.      rx=getchar();ra=;
  12.      while(rx<''||rx>'')rx=getchar();
  13.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  14.  }
  15.  inline int get(int x){
  16.      int l=,r=size,mid;
  17.      while(l<r){
  18.          mid=(l+r+)>>;
  19.          if(mp[mid]<=x)l=mid;else r=mid-;
  20.      }return l;
  21.  }
  22.  int main(){
  23.      n=read();
  24.      for(i=;i<=n;i++)l[i]=mp[i]=read(),r[i]=mp[i+n]=read();
  25.      sort(mp+,mp++n*);size=;n<<=;
  26.      for(i=;i<=n;i++)if(mp[i]!=mp[i-])mp[++size]=mp[i];n>>=;
  27.      for(i=;i<=n;i++)sum[get(l[i])]++,sum[get(r[i])+]--;
  28.      for(i=;i<=size;i++){
  29.          nowsum+=sum[i];
  30.          if(nowsum>ans)ans=nowsum;
  31.      }
  32.      printf("%d\n",ans);
  33.      return ;
  34.  }

bzoj 1915： [Usaco2010 Open]奶牛的跳格子游戏

　　dp+单调队列优化 http://www.cnblogs.com/czllgzmzl/p/5084120.html

bzoj 1729： [Usaco2005 dec]Cow Patterns 牛的模式匹配

　　正解似乎是kmp。。。具体见网上题解

　　然而我偷懒写了hash= =。。。结果卧槽又长又慢

　　如何判断网上题解挺详细的。。。然而事实告诉我们其实并不用那么严格。。选几个类似（并好维护些）的指标（比方说轮廓啊什么的）。。符合的话就再随机几个点确定一下是不是真的全都符合= =。。。。这样大概就可以了吧

　　真相是看网上题解时几个字看错TAT。。结果拍了半天。。最后直接再rand几个数判断= =。。写了4个hashTAT

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  #include<ctime>
  6.  #include<cstdlib>
  7.  #define ll long long
  8.  #define ull unsigned int
  9.  using namespace std;
  10.  const int maxn=;
  11.  const int maxk=;
  12.  int sum[],num[],sumb[],tmpsum[],bel[],belb[],pre1[],pre2[],pre0[];
  13.  int b[maxk],preb[maxk];
  14.  int a[maxn],ans[maxn],prea[maxn],presm[maxn],prebg[maxn];
  15.  int i,j,k,n,m,judlen,s;
  16.  int jud[];
  17.  ull jc1[maxk];
  18.  bool u[],u1[maxk];
  19.  int ra;char rx;
  20.  inline int read(){
  21.      rx=getchar();ra=;
  22.      while(rx<''||rx>'')rx=getchar();
  23.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  24.  }
  25.  struct zs{
  26.      int too,pre;
  27.  }e0[maxn],e1[maxn],e2[maxn];
  28.  int l0[maxn],l1[maxn],l2[maxn],t0,t1,t2;
  29.  ull nowh1=,nowh2=,nowh0=;
  30.  ull hash=,hnow=,jc,h1=,h2=,h0=;
  31.  int kindb,kind;
  32.  inline void ins0(int a,int b){e0[++t0].too=b;e0[t0].pre=l0[a];l0[a]=t0;}
  33.  inline void ins1(int a,int b){e1[++t1].too=b;e1[t1].pre=l1[a];l1[a]=t1;}
  34.  inline void ins2(int a,int b){e2[++t2].too=b;e2[t2].pre=l2[a];l2[a]=t2;}
  35.  inline void clr(int x){
  36.      int i,to;
  37.      for(i=l0[x];i;i=e0[i].pre){
  38.          to=e0[i].too;prea[to]=;if(to>x+k)continue;
  39.          nowh0-=(ull)(to-x)*jc1[x+k-to];
  40.      }
  41.      for(i=l1[x];i;i=e1[i].pre){
  42.          to=e1[i].too;presm[to]=;if(to>x+k)continue;
  43.          nowh1-=(ull)(to-x)*jc1[x+k-to];
  44.      }
  45.      for(i=l2[x];i;i=e2[i].pre){
  46.          to=e2[i].too;prebg[to]=;if(to>x+k)continue;
  47.          nowh2-=(ull)(to-x)*jc1[x+k-to];
  48.      }
  49.  }
  50.  ull mphash=,mphashb[],jjc;
  51.  int main(){
  52.      srand();jc=jc1[]=;
  53.      n=read();k=read();s=read();
  54.      judlen=min(k,min(,/n));
  55.      for(i=;i<=judlen;i++){
  56.          for(j=rand()%k;u1[j];j=rand()%k);
  57.          jud[i]=j;u1[j]=;
  58.      }
  59.      for(i=;i<=n;i++)a[i]=read();
  60.      for(i=;i<=k;i++){
  61.          b[i]=read();u[b[i]]=;
  62.          if(i>){hash*=;hash+=b[i]>=b[i-]?:;jc*=;}
  63.          presm[i]=pre1[b[i]];prebg[i]=pre2[b[i]];preb[i]=pre0[b[i]];
  64.  
  65.          for(j=;j<b[i];j++)pre2[j]=i;for(j=b[i]+;j<=s;j++)pre1[j]=i;pre0[b[i]]=i;
  66.  
  67.          h1*=;if(presm[i])h1+=i-presm[i];
  68.          h2*=;if(prebg[i])h2+=i-prebg[i];
  69.          h0*=;if(preb[i])h0+=i-preb[i];jc1[i]=jc1[i-]*;
  70.      }
  71.      for(i=;i<=s;i++)if(u[i])kindb++,belb[i]=kindb;
  72.      for(i=;i<=k;i++)b[i]=belb[b[i]],sumb[b[i]]++;
  73.      memset(pre1,,sizeof(pre1));memset(pre2,,sizeof(pre2));memset(pre0,,sizeof(pre0));
  74.      for(i=;i<=n;i++){
  75.          presm[i]=pre1[a[i]];prebg[i]=pre2[a[i]];prea[i]=pre0[a[i]];
  76.  
  77.          pre0[a[i]]=i;for(j=;j<a[i];j++)pre2[j]=i;for(j=a[i]+;j<=s;j++)pre1[j]=i;pre0[a[i]]=i;
  78.          ins0(prea[i],i);ins1(presm[i],i);ins2(prebg[i],i);
  79.      }
  80.      for(i=;i<=k;i++){
  81.          if(i>){hnow*=;hnow+=a[i]>=a[i-]?:;}
  82.          kind+=!sum[a[i]],sum[a[i]]++;
  83.          nowh1*=;if(presm[i])nowh1+=i-presm[i];
  84.          nowh2*=;if(prebg[i])nowh2+=i-prebg[i];
  85.          nowh0*=;if(prea[i])nowh0+=i-prea[i];
  86.      }
  87.      int tmp;jjc=;
  88.      for(i=;i<=k;i++){
  89.          for(j=;j<;j++)mphashb[j]*=,mphashb[j]+=b[i]+j-;
  90.          mphash*=,mphash+=a[i];jjc*=;
  91.      }
  92.      int val[]={,};
  93.      for(i=k;i<=n;
  94.          i++,hnow*=,hnow+=val[a[i]>=a[i-]],hnow-=jc*val[a[i-k+]>=a[i-k]]
  95.          ,sum[a[i-k]]--,kind-=(!sum[a[i-k]])-(!sum[a[i]]),sum[a[i]]++
  96.          ,nowh1*=,nowh1+=presm[i]?i-presm[i]:
  97.          ,nowh2*=,nowh2+=prebg[i]?i-prebg[i]:
  98.          ,nowh0*=,nowh0+=prea[i]?i-prea[i]:
  99.          ,clr(i-k)
  100.          ,mphash*=,mphash+=a[i],mphash-=jjc*a[i-k]
  101.          )
  102.      if(hnow==hash&&kind==kindb&&nowh1==h1&&nowh0==h0&&nowh2==h2){
  103.          bool f1=;
  104.          for(j=;j<;j++)if(mphash==mphashb[j]){ans[++ans[]]=i-k+;f1=;break;}
  105.          if(!f1)continue;
  106.          for(j=;j<;j++)if(mphash==mphashb[j]){ans[++ans[]]=i-k+;f1=;break;}
  107.          if(!f1)continue;
  108.          for(tmp=,j=;j<=s;j++)if(sum[j])if(sum[j]!=sumb[++tmp]){f1=;break;}
  109.              else bel[j]=tmp;
  110.          if(!f1)continue;
  111.          for(j=;j<=judlen;j++)if(bel[a[i-jud[j]]]!=b[k-jud[j]]){f1=;break;}
  112.          if(!f1)continue;
  113.          ans[++ans[]]=i-k+;
  114.      }
  115.      printf("%d\n",ans[]);
  116.      for(i=;i<=ans[];i++)printf("%d\n",ans[i]);
  117.      return ;
  118.  }

bzoj 1780： [Usaco2010 Feb]corral 覆盖牛棚

　　首先把被完全覆盖的围栏去掉，剩下的按左端点升序排序。那么右端点肯定也是升序的了。。然后计算出每段围栏，它接下去一段围栏可达到的最远距离。

　　枚举起点，贪心地一段一段接下去就可得到该起点的最优解。

　　直接这样做显然会T。。我们找出一圈围栏后，把这些围栏都指向刚好跨过一圈的那条围栏就可以避免同个解不同起点的重复计算了（看了kpm代码才会的QAQ）。

　　复杂度不会算...

大概是#1吧。。。

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  using namespace std;
  6.  const int maxn=;
  7.  struct fence{
  8.      int l,r;
  9.  }fen[maxn];
  10.  int next[maxn],pos[maxn],jump[maxn],tmpcos[maxn];
  11.  int i,j,L,n,m,ans,k;
  12.  int ra;char rx;
  13.  inline int read(){
  14.      rx=getchar();ra=;
  15.      while(rx<''||rx>'')rx=getchar();
  16.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  17.  }
  18.  bool cmp1(fence a,fence b){return (a.l<b.l)||(a.l==b.l&&a.r>b.r);}
  19.  inline void get(int x){
  20.      if(jump[x])return;
  21.      if(fen[x].r>=L)jump[x]=x,tmpcos[x]=;else get(next[x]),jump[x]=jump[next[x]],tmpcos[x]=tmpcos[next[x]]+;
  22.  }
  23.  int main(){
  24.      L=read();n=read();
  25.      for(i=;i<=n;i++)fen[i].l=read(),fen[i].r=read()+fen[i].l;
  26.      sort(fen+,fen++n,cmp1);
  27.      int mx=-;int tmp=;
  28.      for(i=;i<=n;i++)if(fen[i].r>mx)fen[++tmp]=fen[i],mx=fen[i].r;n=tmp;
  29.      tmp=;
  30.      for(i=;i<=n;i++){
  31.          next[i]=next[i-];
  32.          for(;tmp<n&&fen[tmp+].l<=fen[i].r;tmp++);
  33.          if(tmp==i)next[i]=;else next[i]=tmp,pos[i]=fen[tmp].r;
  34.      }
  35.      int tmpsum,now;
  36.      ans=;
  37.      for(i=;i<=n;i++){
  38.          if(!jump[i])get(i);tmpsum=tmpcos[i]+;now=jump[i];
  39.          while(next[now]&&fen[now].r<fen[i].l+L&&tmpsum<ans)now=next[now],tmpsum++;
  40.          if(fen[now].r>=fen[i].l+L&&tmpsum<ans)ans=tmpsum;
  41.      }
  42.      printf("%d\n",ans);
  43.      return ;
  44.  }

　　正常做法的话大概就是 每个围栏 向 它接下去最远的围栏连条有向边，然后就是基环内向树上的询问了?

bzoj 1698: [Usaco2007 Feb]Lilypad Pond 荷叶池塘

　　最短路计数。。题解说得挺详细的。。

　　一开始想到把连到空格子的边权设为1，结果发现这样会重复计算方案（兜了一圈边权0的再回来，放的荷叶并没有变化）

　　那么根据网上的题解TAT，处理出两两空格子（包括起点和终点）之间能否不经其他空格子到达。。可以到达的话就连边，权值为1。

　　在新图中跑最短路计数就行了。KPM实力#1。。。伏地膜

　　注意因为新图中的边权都为1所以最短路计数随便搞。。。如果是一般些的图的话得考虑重复计数的情况。

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #define ll long long
  5.  using namespace std;
  6.  const int maxn=;
  7.  const int xx[]={-,-,,,,,-,-},yy[]={-,-,-,-,,,,};
  8.  struct zs{
  9.      int too,pre;
  10.  }e[maxn*maxn*maxn*maxn];
  11.  int map[maxn][maxn];
  12.  int id[maxn][maxn],poi;
  13.  int last[maxn*maxn],tot;
  14.  int l,r,now,i,j,k,nowx,nowy,x,y,n,m,a,S,T,to;
  15.  int dlx[maxn*maxn],dly[maxn*maxn],dis[maxn*maxn];
  16.  ll num[maxn*maxn];
  17.  bool u[maxn*maxn];
  18.  inline void insert(int a,int b){
  19.      e[++tot].too=b;e[tot].pre=last[a];last[a]=tot;
  20.  //  printf("    %d-->%d\n",a,b);
  21.  }
  22.  
  23.  int ra;char rx;
  24.  inline int read(){
  25.      rx=getchar();ra=;
  26.      while(rx<''||rx>'')rx=getchar();
  27.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  28.  }
  29.  int main(){
  30.      n=read();m=read();
  31.      for(i=;i<=n;i++)for(j=;j<=m;j++){
  32.          a=read(),id[i][j]=++poi,map[i][j]=(a==||a==||a==)?:a;
  33.          if(a==)S=poi;if(a==)T=poi;
  34.      }
  35.      for(i=;i<=n;i++)for(j=;j<=m;j++)if(!map[i][j]){
  36.          memset(u,,poi+);
  37.          l=;r=;dlx[]=i;dly[]=j;u[id[i][j]]=;
  38.          while(l<r){
  39.              nowx=dlx[++l];nowy=dly[l];
  40.              for(k=;k<;k++){
  41.                  x=nowx+xx[k];y=nowy+yy[k];
  42.                  if(x<||y<||x>n||y>m||u[id[x][y]]||map[x][y]==)continue;u[id[x][y]]=;
  43.                  if(map[x][y])dlx[++r]=x,dly[r]=y;else insert(id[i][j],id[x][y]);
  44.              }
  45.          }
  46.      }
  47.      memset(u,,poi+);
  48.      l=;r=;dlx[]=S;num[S]=;dis[S]=;
  49.      while(l<r){
  50.          now=dlx[++l];u[now]=;
  51.          if(now==T)break;
  52.          for(i=last[now],to=e[i].too;i;i=e[i].pre,to=e[i].too)if(!u[to])
  53.              if(!num[to])dlx[++r]=to,num[to]=num[now],dis[to]=dis[now]+;
  54.              else if(dis[to]==dis[now]+)num[to]+=num[now];
  55.      }
  56.      if(num[T])printf("%d\n",dis[T]-),printf("%lld\n",num[T]);
  57.      else puts("-1");
  58.  //  printf("  %d\n",tot);
  59.      return ;
  60.  }

bzoj 1737： [Usaco2005 jan]Naptime 午睡时间

　　dp。。设f[i][j][0]表示前i个时间段里，睡了j个时间段的最大值，但时间i没睡。f[i][j][1]表示时间i睡了，其他一样。val[i]表示在时间i睡觉的总效用值。

　　f[i][j][0]=max{　f[i-1][j][0],f[i-1][j][1]　}

　　f[i][j][1]=max{　f[i-1][j-1][0],f[i-1][j-1][1]+val[i]　}。

　　第一次我们强制不跨环，初始化f[1][1][1]=f[1][0][0]=0，此时答案是f数组最大值；第二次我们强制跨环，初始化f[1][1][1]=0,f[1][0][0]=-inf，此时答案是f[n][m][1]+val[1]。

　　最终答案就是两种情况的较大值。

又是#1辣。。

  1.  #include<cstdio>
  2.  #include<cstring>
  3.  #include<iostream>
  4.  using namespace std;
  5.  const int maxn=;
  6.  int f[maxn][],pre;
  7.  int val[maxn];
  8.  int i,j,n,m,ans;
  9.  int ra;char rx;
  10.  inline int read(){
  11.      rx=getchar();ra=;
  12.      while(rx<''||rx>'')rx=getchar();
  13.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  14.  }
  15.  int main(){
  16.      n=read();m=read();
  17.      for(i=;i<=n;i++)val[i]=read();
  18.      memset(f,,(m+)<<);
  19.      f[][]=f[][]=;
  20.      for(i=;i<=n;i++)for(j=m;j;j--){
  21.          if(f[j][]>f[j][])f[j][]=f[j][];
  22.          f[j][]=f[j-][]+val[i];if(f[j-][]>f[j][])f[j][]=f[j-][];
  23.      }ans=max(f[m][],f[m][]);
  24.      memset(f,,(m+)<<);
  25.      f[][]=;
  26.      for(i=;i<=n;i++)for(j=m;j;j--){
  27.          if(f[j][]>f[j][])f[j][]=f[j][];
  28.          f[j][]=f[j-][]+val[i];if(f[j-][]>f[j][])f[j][]=f[j-][];
  29.      }ans=max(ans,f[m][]+val[]);
  30.      printf("%d\n",ans);
  31.      return ;
  32.  }

bzoj 3312: [Usaco2013 Nov]No Change

　　一开始以为是贪心。。结果看数据范围果断状压= =

　　f[i]表示已取硬币状态为i时，最多可支付的东西数。pre[i]表示前i个东西的价格和。val[i]表示第i个硬币的面值

　　f[i]=max{　f[i-(1<<k)]+x　}，（(1<<k)存在于状态i中，x是使pre[ x+f[i-(1<<k) ]-pre[ f[i-(1<<k)] ] <=val[k]的最大值，也就是最多能再买多少件东西  ）。

　　时间复杂度(2^k*k*logn)

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  using namespace std;
  5.  const int maxn=;
  6.  int pre[maxn],val[],two[];
  7.  int f[(<<)+],pos[(<<)+];
  8.  int i,j,K,n,l,r,mid,k,st,tmp,ans;
  9.  int ra;char rx;
  10.  inline int read(){
  11.      rx=getchar();ra=;
  12.      while(rx<''||rx>'')rx=getchar();
  13.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  14.  }
  15.  int main(){
  16.      K=read();n=read();two[]=;pos[]=;
  17.      for(i=;i<K;i++)val[i]=read();
  18.      for(i=;i<K;i++)two[i]=two[i-]<<,pos[two[i]]=i;
  19.      for(i=;i<=n;i++)pre[i]=read()+pre[i-];
  20.      f[]=;int mx=<<K;
  21.      ans=-;
  22.      for(i=;i<mx;i++){
  23.          for(tmp=,k=mx-i-,j=pos[k&-k];k;k-=k&-k,j=pos[k&-k])tmp+=val[j];
  24.          if(tmp<=ans)continue;
  25.          for(k=i,j=pos[k&-k];k;k-=k&-k,j=pos[k&-k]){
  26.              l=st=f[i^(k&-k)];r=n;
  27.              while(l<r&&r>f[i]){
  28.                  mid=(l+r+)>>;
  29.                  if(pre[mid]-pre[st]>val[j])r=mid-;else l=mid;
  30.              }
  31.              if(l>f[i])f[i]=l;
  32.          }
  33.          if(f[i]==n)ans=tmp;
  34.  //      printf("   %d %d\n",i,f[i]);
  35.      }
  36.      printf("%d\n",ans);
  37.      return ;
  38.  }

bzoj 1750: [Usaco2005 qua]Apple Catching

　　傻逼dp。。f[i][j][k]表示i秒后，牛已移动j次，在k处（k=0在左，k=1在右）

　　f[i][j][k]=max{　f[i-1][j][k],f[i-1][j-1][1-k]　}+map[i][k]，（map[i][k]表示k处在第i秒的时候有无苹果。）

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  using namespace std;
  5.  int f[][];
  6.  int i,j,x,n,m,ans;
  7.  int ra;char rx;
  8.  inline int read(){
  9.      rx=getchar();ra=;
  10.      while(rx<''||rx>'')rx=getchar();
  11.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  12.  }
  13.  int main(){
  14.      n=read();m=read();
  15.      f[][]=-;
  16.      for(i=;i<=n;i++){
  17.          x=read();x--;
  18.          if(x)for(j=m;j;j--){
  19.              if(f[j-][]>f[j][])f[j][]=f[j-][];
  20.              if(f[j-][]>f[j][])f[j][]=f[j-][];
  21.              f[j][]++;
  22.          }else for(j=m;j;j--){
  23.              if(f[j-][]>f[j][])f[j][]=f[j-][];
  24.              if(f[j-][]>f[j][])f[j][]=f[j-][];
  25.              f[j][]++;
  26.          }
  27.          f[j][x]++;
  28.      }
  29.      ans=max(f[m][],f[m][]);
  30.      for(i=m-;i>=;i--){
  31.          if(f[i][]>ans)ans=f[i][];
  32.          if(f[i][]>ans)ans=f[i][];
  33.      }
  34.      printf("%d\n",ans);
  35.      return ;
  36.  }

bzoj 1746： [Usaco2005 open]Lazy Cows

　　再次验证了我是傻逼= =

　　因为只有两行，所以帆布的状态只有四种：只覆盖上面一行，只覆盖下面一行，两行用一块布覆盖，两行用两块布分别覆盖。

　　用f[i][j][0..3]表示用j块布覆盖了前i个洞后，对应以上四种状态的最小面积。预处理pos[i]表示第i个破洞的列数。

　　f[i][0]=min(　　min{ f[i-1][j][0],f[i-1][j][3] }+pos[i]-pos[i-1]　　,　　min{ f[i-1][j-1][0..3] }+1　　)，（第i个洞在第一行）

　　f[i][1]基本同上。

　　f[i][2]=min{　　f[i-1][j][2]+2*(pos[i]-pos[i-1])　　,　　min{ f[i-1][j-1][0..3] }+2　　}

　　f[i][3]=min{　　min{ f[i-1][j-1][0..1] }+pos[i]-pos[i-1]+1　　,　　min{ f[i-1][j-2][0..3] }+2　　,　　f[i-1][j][3]+2*(pos[i]-pos[i-1])　　}

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  using namespace std;
  6.  const int maxn=;
  7.  const int inf=;
  8.  struct poi{
  9.      bool x;int y;
  10.  }a[maxn];
  11.  int num[maxn],dis[maxn];
  12.  int f[maxn][];
  13.  int n,m,flag,K,mn,ans;
  14.  int ra;char rx;
  15.  inline int read(){
  16.      rx=getchar();ra=;
  17.      while(rx<''||rx>'')rx=getchar();
  18.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  19.  }
  20.  bool cmp(poi a,poi b){return a.y<b.y;}
  21.  int main(){
  22.      int i,j,now;
  23.      n=read();K=read();m=read();
  24.      for(i=;i<=n;i++)a[i].x=read()==,a[i].y=read();
  25.      sort(a+,a++n,cmp);int tmp=;a[].y=a[].y;
  26.      for(i=;i<=n;){
  27.          flag=a[i].x+;
  28.          for(j=i+;j<=n&&a[j].y==a[i].y;j++)flag|=a[j].x+;
  29.          num[++tmp]=flag-;
  30.          dis[tmp]=a[i].y-a[i-].y;
  31.          i=j;
  32.      }n=tmp;
  33.      memset(f,,(n+)<<);
  34.      f[][num[]]=max(num[],);f[][]=f[][]=;
  35.      for(i=;i<=n;i++)for(now=dis[i],j=K-,mn=min(f[j][],f[j][]<f[j][]?f[j][]:f[j][]),f[++j][]+=now,f[j][]+=now,f[j][]+=now,f[j][]+=now;j;){
  36.          if(!num[i]){
  37.              if(mn<f[j][])f[j][]=mn+;
  38.              if(f[j][]<f[j][])f[j][]=f[j][];
  39.              f[j][]=inf;
  40.          }else if(num[i]==){
  41.              if(mn<f[j][])f[j][]=mn+;
  42.              if(f[j][]<f[j][])f[j][]=f[j][];
  43.              f[j][]=inf;
  44.          }else f[j][]=f[j][]=inf;
  45.  
  46.          f[j][]+=now;if(mn+<f[j][])f[j][]=mn+;
  47.          f[j][]+=now;
  48.          j-=;if(j>=){
  49.              mn=f[j][];if(f[j][]<mn)mn=f[j][];if(f[j][]<mn)mn=f[j][];if(f[j][]<mn)mn=f[j][];
  50.              if(mn+<f[j+][])f[j+][]=mn+;
  51.          }
  52.          f[++j][]+=now;f[j][]+=now;f[j][]+=now;f[j][]+=now;
  53.          if(f[j][]<f[j+][])f[j+][]=f[j][]+;
  54.          if(f[j][]<f[j+][])f[j+][]=f[j][]+;
  55.      }
  56.      ans=inf;
  57.      i=min(n,K);
  58.      for(j=;j<;j++)if(f[i][j]<ans)ans=f[i][j];
  59.      printf("%d\n",ans);
  60.      return ;
  61.  }

bzoj 1727: [Usaco2006 Open]The Milk Queue 挤奶队列

　　抱了ctl学长的大腿。。。

　　贪心。。如果牛A排在牛B前面比B排在A前面更优所以blabla。。具体列个式子后快排就好了。。

  1.  #include<iostream>
  2.  #include<cstring>
  3.  #include<cstdio>
  4.  #include<algorithm>
  5.  using namespace std;
  6.  const int maxn=;
  7.  int id[maxn],a[maxn],b[maxn];
  8.  int i,end1,end2,j,n;
  9.  int ra;char rx;
  10.  inline int read(){
  11.      rx=getchar();ra=;
  12.      while(rx<''||rx>'')rx=getchar();
  13.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  14.  }
  15.  bool cmp(int x,int y){return a[x]+b[y]+(b[x]>a[y]?b[x]:a[y])<a[y]+b[x]+(b[y]>a[x]?b[y]:a[x]);}
  16.  int main(){
  17.      n=read();
  18.      for(i=;i<=n;i++)a[i]=read(),b[i]=read(),id[i]=i;
  19.      sort(id+,id++n,cmp);
  20.      for(i=;i<=n;i++){
  21.          j=id[i];end1+=a[j];
  22.          if(end1<=end2)end2+=b[j];else end2=end1+b[j];
  23.      }
  24.      printf("%d\n",max(end1,end2));
  25.      return ;
  26.  }

bzoj 1736： [Usaco2005 jan]The Wedding Juicer 婚宴的榨汁机

　　从图的边界开始往内部跑类似最短路的。。dis[i]表示i点最大可接受高度。

　　dis[i]也就是从图的边界上的任意点出发，高度不递减地跑到i，途中经过的最大高度的最小值（就是水流出去的过程反过来）。

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  #include<queue>
  6.  #define ll long long
  7.  using namespace std;
  8.  const int maxn=;
  9.  const int xx[]={,,,-},yy[]={,-,,};
  10.  int map[maxn][],num[maxn],poinum,id[][];
  11.  int dis[maxn],mp[maxn];
  12.  bool used[maxn],full[maxn];
  13.  int i,j,k,n,m,nx,ny,x,y,mx,now,size,to;
  14.  ll ans;
  15.  
  16.  struct poi{int pos,dis;};
  17.  priority_queue<poi>q;
  18.  bool operator <(poi a,poi b){return a.dis>b.dis;}
  19.  
  20.  int ra;char rx;
  21.  inline int read(){
  22.      rx=getchar();ra=;
  23.      while(rx<''||rx>'')rx=getchar();
  24.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  25.  }
  26.  int main(){
  27.      m=read();n=read();
  28.      for(i=;i<=n;i++)for(j=;j<=m;j++)id[i][j]=++poinum,mp[poinum]=read(),mx=max(mx,mp[poinum]);
  29.      for(i=;i<=n;i++)for(j=;j<=m;j++)for(k=,x=id[i][j];k<;k++){
  30.          nx=i+xx[k];ny=j+yy[k];if(nx<||nx>n||ny<||ny>m)continue;
  31.          y=id[nx][ny];map[x][num[x]++]=y;
  32.      }
  33.      for(i=;i<=n;i++)for(j=;j<=m;j++)
  34.          if(i==||i==n||j==||j==m)dis[id[i][j]]=mp[id[i][j]],full[id[i][j]]=;
  35.          else dis[id[i][j]]=mx;
  36.      for(i=;i<m;i++)q.push((poi){i,mp[i]}),q.push((poi){id[n][i],mp[id[n][i]]}),size+=;
  37.      for(i=;i<=n;i++)q.push((poi){id[i][],mp[id[i][]]}),q.push((poi){id[i][m],mp[id[i][m]]}),size+=;
  38.      while(size){
  39.          while(size&&used[q.top().pos])q.pop(),size--;
  40.          if(!size||q.top().dis==mx)break;
  41.          now=q.top().pos;x=q.top().dis;used[now]=;q.pop();size--;
  42.          for(i=;i<num[now];i++)if(dis[to=map[now][i]]>x&&!full[to]){
  43.              if(mp[to]>x)dis[to]=mp[to],full[to]=;else dis[to]=x;
  44.              q.push((poi){to,dis[to]});size++;
  45.          }
  46.      }
  47.      for(i=;i<=n;i++)for(j=;j<=m;j++)if(!full[id[i][j]])ans+=(ll)dis[id[i][j]]-mp[id[i][j]];
  48.      printf("%lld\n",ans);
  49.      return ;
  50.  }

bzoj 3408： [Usaco2009 Oct]Heat Wave 热浪

　　裸最短路。。。为何这么少人写

　　又一次调优先队列结果和正常的spfa一样快。

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<queue>
  5.  using namespace std;
  6.  struct zs{int too,pre,dis;}e[];
  7.  struct poi{int pos,dis;};
  8.  priority_queue<poi>q;
  9.  bool operator <(poi a,poi b){return a.dis>b.dis;}
  10.  int last[],dis[];
  11.  bool used[];
  12.  int i,j,s,t,a,b,x,y,z,n,m,tot,size,to;
  13.  int ra;char rx;
  14.  inline int read(){
  15.      rx=getchar();ra=;
  16.      while(rx<''||rx>'')rx=getchar();
  17.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  18.  }
  19.  int main(){
  20.      n=read();m=read();s=read();t=read();
  21.      for(;m;m--){
  22.          a=read(),b=read();
  23.          e[++tot].dis=e[tot+].dis=read();
  24.          e[tot].too=b;e[tot].pre=last[a];last[a]=tot;
  25.          e[++tot].too=a;e[tot].pre=last[b];last[b]=tot;
  26.      }
  27.      memset(dis,,(n+)<<);dis[s]=;size=;q.push((poi){s,});
  28.      while(size&&!used[t]){
  29.          while(size&&used[q.top().pos])q.pop(),size--;
  30.          if(!size||dis[q.top().pos]>=dis[t]||used[t])break;
  31.          x=q.top().pos;y=q.top().dis;q.pop();--size;used[x]=;
  32.  
  33.          for(to=e[i=last[x]].too,z=y+e[i].dis;i;to=e[i=e[i].pre].too,z=y+e[i].dis)if(dis[to]>z&&z<dis[t])
  34.              q.push((poi){to,dis[to]=z}),size++;
  35.      }
  36.      printf("%d\n",dis[t]);
  37.      return ;
  38.  }

bzoj 2590： [Usaco2012 Feb]Cow Coupons

　　排序显然是错的= =懒得重新写了>_<

　　如果没有优惠券的话直接排序后贪心。有了优惠券后，把每头牛拆成有优惠券和没优惠券两种，把2n个数排序，注意如果买了一头优惠的牛，把它没优惠券的版本设为不可取。

　　无法理解为啥sort比kpm调优先队列慢TAT（虽然要排2n个数。。但常数应该是比优先队列小很多的吧。。）

  1.  #include<cstdio>
  2.  #include<iostream>
  3.  #include<cstring>
  4.  #include<algorithm>
  5.  #define ll long long
  6.  using namespace std;
  7.  bool u[];
  8.  int id[],a[];
  9.  int i,j,n,k,ans,n1;
  10.  ll m;
  11.  int ra;char rx;
  12.  inline int read(){
  13.      rx=getchar();ra=;
  14.      while(rx<''||rx>'')rx=getchar();
  15.      while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
  16.  }
  17.  bool cmp(int x,int y){return a[x]<a[y];}
  18.  int main(){
  19.      n=read();k=read();scanf("%lld",&m);
  20.      for(i=;i<=n;i++)id[i]=i,id[i+n]=i+n,a[i]=read(),a[i+n]=read();
  21.      n1=n<<;
  22.      sort(id+,id++n1,cmp);
  23.      for(i=,j=id[];i<=n1&&m>;j=id[++i])if(m>=a[j]&&!u[j>n?j-n:j]&&(j<=n||k))
  24.      m-=a[j],u[j>n?j-n:j]=,ans++,k-=j>n;
  25.      printf("%d\n",ans);
  26.      return ;
  27.  }

完结撒花[鼓掌熊]总共大概106道吧。。

bzoj usaco 金组水题题解（2.5）的更多相关文章

  1. bzoj usaco 金组水题题解（1）

    UPD:我真不是想骗访问量TAT..一开始没注意总长度写着写着网页崩了王仓(其实中午的时候就时常开始卡了= =)....损失了2h(幸好长一点的都单独开了一篇)....吓得赶紧分成两坨....TAT. ...

  2. bzoj usaco 金组水题题解（2）

    续.....TAT这回不到50题编辑器就崩了.. 这里塞40道吧= = bzoj 1585: [Usaco2009 Mar]Earthquake Damage 2 地震伤害 比较经典的最小割?..然而 ...

  3. BZOJ USACO 银组 水题集锦

    最近刷银组刷得好欢快，好像都是水题，在这里吧他们都记录一下吧（都是水题大家一定是道道都虐的把= =）几道比较神奇的题到时再列出来单独讲一下吧= =（其实我会说是BZOJ蹦了无聊再来写的么 = =） [ ...

  4. NOIP2018初赛普及组原题&题解

    NOIP2018初赛普及组原题&题解 目录 NOIP2018初赛普及组原题&题解 原题&答案 题解 单项选择题 第$1$题 第$2$题 第$3$题 第$4$题 第$5$题 第$ ...

  5. World Finals 2017 (水题题解)

    看大佬做2017-WF,我这种菜鸡,只能刷刷水题,勉强维持生活. 赛后补补水题. 题目pdf链接,中文的,tls翻译的,链接在这里 个人喜欢在vjudge上面刷题. E Need for Speed ...

  6. noip2008普及组3题题解-rLq

    (第一次写题解,随意喷) (只是前一天的作业哈) (先凑个数) 题目描述 上体育课的时候,小蛮的老师经常带着同学们一起做游戏.这次,老师带着同学们一起做传球游戏.游戏规则是这样的:n个同学站成一个圆圈 ...

  7. 2006-2007 ACM-ICPC | POJ3380 POJ3384 POJ3385 水题题解

    // CF比赛链接:http://codeforces.com/gym/101650 // POJ链接:http://poj.org/searchproblem?field=source&ke ...

  8. noip2008普及组4题题解-rLq

    (啊啊啊终于补到了今天的作业了) 本题地址:http://www.luogu.org/problem/show?pid=1058 题目描述 小渊是个聪明的孩子,他经常会给周围的小朋友们将写自己认为有趣 ...

  9. Cmd2001的毒瘤水题题解

    怕不是我再不写题解这题就该成没人做也没人会的千古谜题了...... T1: 仔细分析题面,发现相同就是广义SAM上节点相同,相似就是广义SAM上为从根到某个点路径的前缀..直接SAM上跑从根开始,每个 ...

随机推荐

  1. Vagrant安装完lnmp后，配置linux和windows共享文件并配置虚拟主机访问项目

    虚拟机目录下的Vagrantfile文件是vagrant的配置文件,如果想把虚拟机当作一台服务器,可以通过ip访问,需要修改配置文件进行配置. (1)第一步:打开虚拟机目录下的Vagrantfile文 ...

  2. Mac中Eclipse安装和使用svn

    Eclipse版本为Neon Release (4.6.0) 安装svn 安装HomeBrew 在终端中输入 ruby -e "$(curl -fsSL https://raw.github ...

  3. awk 字符拼接

    i-6gdkkemx 2 2 #cat cpu.txt | awk '{print $1,$2"核"$3"G"}' i-6gdkkemx 2核2G

  4. Kotlin——最详细的数据类、密封类详解

    在前面几个章节章节中,详细的讲解了Koltin中的接口类(Interface).枚举类(Enmu),还不甚了解的可以查看我的上一篇文章Kotlin--接口类.枚举类详解.当然,在Koltin中,除了接 ...

  5. 496. Next Greater Element I

    You are given two arrays(without duplicates)nums1andnums2wherenums1's elements are subset ofnums2. F ...

  6. Hawk原理：通过IEnumerable实现通用的ETL管道

    针对IEnumerable已经有多篇文章,本篇介绍如何使用IEnumerable实现ETL. ETL,是英文 Extract-Transform-Load 的缩写,用来描述将数据从来源端经过萃取(ex ...

  7. PHP （超文本预处理器）

    PHP(外文名:PHP: Hypertext Preprocessor,中文名:"超为本预处理器")是一种通用开源脚本语言.语法吸收了C语言.java和Rerl的特点,利于学习,使 ...

  8. python基础知识——字符串详解

    大多数人学习的第一门编程语言是C/C++,个人觉得C/C++也许是小白入门的最合适的语言,但是必须承认C/C++确实有的地方难以理解,初学者如果没有正确理解,就可能会在使用指针等变量时候变得越来越困惑 ...

  9. BCL和CoreFx的区别

    bcl是.netframework clr 的基础库corefx是.net core clr的基础库

  10. Ajax异步提交登录(2)--登录使用

    http://cjp1989.iteye.com/blog/1740964 1.Ajax的原理: Ajax的原理就是:通过javascript的方式,将前台数据通过xmlhttp对象传递到后台,后台在 ...