BZOJ 刷题记录 PART 4

【BZOJ1143】CTSC的题目。。。先用floyed传递闭包，然后直接上匈牙利算法。

【BZOJ1452】从未写过的二维树状数组。好像很简单。。

struct two_bit
{
  int f[305][305];
  inline void add(int x,int z,int A)
  {
    for (;x<=n;x+=L(x))
      for (int y=z;y<=m;y+=L(y))
        f[x][y]+=A;
  }
  inline int ask(int x,int z)
  {
    int ans=0;
    for (;x;x-=L(x))
      for (int y=z;y;y-=L(y))
        ans+=f[x][y];
    return ans;
  }
}a[105];

【BZOJ1876】高精GCD。但是实现卡的很紧。如果是裸的GCD，是否超时不论，光是代码长度就令人毛骨悚然。（详见我的VIJOS1047题解）于是可以用二进制GCD来优化。

若A,B是偶数，GCD(A,B)=GCD(A/2,B/2)

若A是奇数，B是偶数，GCD(A,B)=GCD(A,B/2)

若A是奇数，B是奇数，GCD(A,B)=GCD(A-B,B)

直到A,B中的一个是0，或者A=B。

【BZOJ3038&3211】开始还以为是神题，后来发现了一个性质：一个数即使暴力根号操作，几次后就会变成1。（主要还要特判0的情况！！！）于是我们可以用线段树暴力修改。

【BZOJ1695&1743】这道权限题萎了我好长时间。开始我倒着走，在字母树中有子节点更新父节点——但是狂WA不止。于是开始更新子节点，总算A了。

for (i=1;i<=n;i++)
    for (j=m;j;j--)
      for (k=1;k<=num[map[i][j]];k++)
      {
        x=record[map[i][j]][k];
        if (flag[x]) for (p=1;p<=j;p++) c[x][p]++;
        for (son=0;son<26;son++)
          if ((t=tree[x][son])>-1)
            for (p=1;p<=j;p++)
              c[x][p]+=c[t][j];
      }
  for (son=0;son<26;son++) 
    if (tree[0][son]>-1)
      ans+=c[tree[0][son]][1];

【BZOJ3170】经典的坐标重构。

for (i=1;i<=n;i++) 
  {
    scanf("%lld%lld",&x,&y);a[i].x=x+y;a[i].y=x-y;
    a[i].id=i;endx+=a[i].x;endy+=a[i].y;
  }

【BZOJ3319】我淡淡地说：暴力RANK1，你能信？

【BZOJ1765&1356】找最大的面积？卡时间！先把所有对角线全部都找出来，排一遍序。从大到小枚举每一个相等的对角线的集合，一但匹配成功就机智的BREAK。

【BZOJ1879】看题目这么神，就是一个简单的状压DP。

f[0][0]=1;
    for (i=0;i<len;i++)
      for (status=0;status<(1<<n);status++)
        if (f[i][status])
          for (j=0;j<26;j++)
            up(f[i+1][status|c[i+1][j]],f[i][status]);

【BZOJ1226】神状压！注意到一个性质：BI<=7,也就是说，对于某一个点，影响他的只有前7个和后7个点。

PS:a or b-a and b不就是a xor b吗？

    f[1][0][-1+A]=0;
    for (i=1;i<=n;i++)
      for (status=0;status<256;status++)
        for (k=-8;k<=7;k++)
          if (f[i][status][k+A]<INF)
          {
            if (status&1) {up(f[i+1][status>>1][k-1+A],f[i][status][k+A]);continue;}
            r=INF;
            for (j=0;j<=B[i];j++)
              if (!(status&(1<<j)))
              {
                if (i+j>r) break;
                up(r,B[i+j]+i+j);
                up(f[i][status|(1<<j)][j+A],f[i][status][k+A]+cost(i+k,i+j));
              }
          }