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));
}
}