zjoi2016 day1【bzoj4455】【bzoj4456】

首先做了T2的旅行者，看到bz上面过的人数比较多。。

考试的时候完全没有想太多。一闪而过了分块思想，然后就没有然后了。。

大视野上面有题解，竟然是一个初中生写的。。？

正解其实是“分治”，每次选择中轴线，不会势能分析，感觉考场上想出来也肯定不敢打。。

#include<cstdio>
#include<algorithm>
#define inf 1000000000
#define N 100010
using namespace std;
int dis[N],ans[N],q[N],id[N],vet[N],head[N],next[N],pri[N];
int edgenum,cas,cnt,n,m;
struct node{int x1,y1,x2,y2,id;}a[N],b[N];
void add(int u,int v,int w)
{
  edgenum++;vet[edgenum]=v;next[edgenum]=head[u];
  head[u]=edgenum;pri[edgenum]=w;
  //printf("%d %d %d\n",u,v,w);
}
void push_down(int k)
{
  <=cnt)
  {
     ],k3=q[k*+];
     +<=cnt&&dis[k1]>dis[k3]&&dis[k3]<dis[k2])
     {
       id[k1]=k*+;id[k3]=k;swap(q[k],q[k*+]);
       k=k*+;
     }else if(dis[k1]>dis[k2])
     {
       id[k1]=k*;id[k2]=k;swap(q[k],q[k*]);k=k*;
     }else break;
  }
}
void push_up(int k)
{
  >&&dis[q[k]]<dis[q[k/]])
  {
    id[q[k]]=k/;id[q[k/]]=k;swap(q[k],q[k/]);
    k=k/;
  }
}
void dijk(int sta,int x1,int x2,int y1,int y2)
{
  cnt=;q[]=sta;id[sta]=;dis[sta]=;
  for(int i=x1;i<=x2;i++)
    for(int j=y1;j<=y2;j++)
     {
       )*m+j;
       if(u!=sta){dis[u]=inf;q[++cnt]=u;id[u]=cnt;}
     }
  while(cnt)
  {
     ];q[]=q[cnt--];id[q[]]=;push_down();int e=head[u];//printf("reng=%d\n",u);
     )
     {
       ,yy=v%m;)yy=m;)xx--;
       //printf("ok=%d %d %d\n",v,xx,yy);
       if(xx>=x1&&xx<=x2&&yy>=y1&&yy<=y2)
       if(dis[v]>dis[u]+pri[e])
       {

         dis[v]=dis[u]+pri[e];
         //printf("ok=%d %d %d\n",u,v,dis[v]);
         push_up(id[v]);
       }
       e=next[e];
     }
  }
  //printf("sta=%d\n",sta);
  for(int i=x1;i<=x2;i++)
    for(int j=y1;j<=y2;j++)
     {
       )*m+j;
       //printf("dis=%d %d\n",u,dis[u]);
     }
}
void solve(int x1,int x2,int y1,int y2,int l,int r)
{
  if(l>r)return;int mid;
  if(x2-x1>y2-y1)
  {
    mid=(x1+x2)>>;
    for(int i=y1;i<=y2;i++)
    {
      )*m+i;dijk(u,x1,x2,y1,y2);
      for(int j=l;j<=r;j++)
      {
        )*m+a[j].y1,u2=(a[j].x2-)*m+a[j].y2;
        int t=dis[u1]+dis[u2];
        ans[a[j].id]=min(ans[a[j].id],t);
      }
    }
      ,k=r+;
      for(int i=l;i<=r;i++)if(a[i].x1<mid&&a[i].x2<mid)b[++j]=a[i];else
                           if(a[i].x1>mid&&a[i].x2>mid)b[--k]=a[i];
      ,y1,y2,l,j);
      ,x2,y1,y2,k,r);
  }else
  {
    mid=(y1+y2)>>;
    for(int i=x1;i<=x2;i++)
    {
      )*m+mid;dijk(u,x1,x2,y1,y2);
      for(int j=l;j<=r;j++)
      {
        )*m+a[j].y1,u2=(a[j].x2-)*m+a[j].y2;
        int t=dis[u1]+dis[u2];
        ans[a[j].id]=min(ans[a[j].id],t);
      }
    }
      ,k=r+;
      for(int i=l;i<=r;i++)if(a[i].y1<mid&&a[i].y2<mid)b[++j]=a[i];else
                           if(a[i].y1>mid&&a[i].y2>mid)b[--k]=a[i];
      ,l,j);
      ,y2,k,r);
  }
}
int main()
{
freopen("4456.in","r",stdin);
freopen("4456.out","w",stdout);
  scanf("%d %d",&n,&m);
  ;i<=n;i++)
   ;j<m;j++)
    {
      )*m+j,v=(i-)*m+j+;scanf("%d",&x);
      add(u,v,x);add(v,u,x);
    }

  ;i<n;i++)
   ;j<=m;j++)
    {
      )*m+j,v=i*m+j;scanf("%d",&x);
      add(u,v,x);add(v,u,x);
    }
  scanf("%d",&cas);
  ;i<=cas;i++)
  {
    scanf("%d%d%d%d",&a[i].x1,&a[i].y1,&a[i].x2,&a[i].y2);
    a[i].id=i;ans[i]=inf;
  }
  solve(,n,,m,,cas);
  ;i<=cas;i++)printf("%d\n",ans[i]);
} 

T3 star小星星

考场的时候竟然没想到状压。。。。好想撞墙呀。。。。

40感觉很好拿。。。。（雾

一开始并不是很能理解题解里“映射”的概念，后来是一边写代码一边想的，感觉这种习惯不太好

个人理解：所有点映射到全集的方案数-所有点映射到某个点集的方案数+所有点映射到某个点集的方案数++--……

这里关于“所有点映射到***的方案数”很好实现，一个树形DP就可以了

比较难想的是这个容斥，感觉容斥功底弱爆了。。。。O（2^n n^3)算了一下感觉不太对（要跑6s？）

但是A得很轻易啊。。。。。最后别忘了数组要开大，否则。。。边表爆炸

#include<cstdio>
#include<algorithm>
#define ll long long
#define N 45
using namespace std;
,cnt;
int vet[N],next[N],n,m,map[N][N],head[N],a[N];
ll dp[N][N],bin[N],sum,ans;
void add(int u,int v)
{
  edgenum++;vet[edgenum]=v;next[edgenum]=head[u];head[u]=edgenum;
}
void dfs(int u,int fa)
{
   int e=head[u];
   )
   {
     int v=vet[e];
     if(v!=fa)dfs(v,u);
     e=next[e];
   }
   ;i<=cnt;i++)
   {
     dp[u][i]=;e=head[u];
     )
     {
       ;if(v==fa){e=next[e];continue;}
       ;j<=cnt;j++)
         if(map[a[i]][a[j]])tmp+=dp[v][j];
       dp[u][i]*=tmp;)break;
       e=next[e];
     }
   }
}
int main()
{
freopen("4455.in","r",stdin);
int u,v;
  scanf("%d%d",&n,&m);
  ;i<=m;i++)
  {
    scanf(;
  }bin[]=;
  ;i<=n;i++)bin[i]=bin[i-]<<;
  ;i<n;i++)
  {
    scanf("%d%d",&u,&v);add(u,v);add(v,u);
  }
  ans=;
  ;i<bin[n];i++)
  {
    sum=cnt=;;k<=n;k++)])a[++cnt]=k;
    dfs(,);;k<=cnt;k++)sum+=dp[][k];
    )ans-=sum;else ans+=sum;
    //printf("%lld\n",ans);
  }
  printf("%lld",ans);
}