B

整个表格其实是一些联通块,取反操作不能跨连通块。所以直接统计一下每个连通块内数字不对的个数是不是偶数即可


#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define gt(x,y) ((x)*m+y)
using namespace std; const int M = 1100000;
int n,m,k,a[M],b[M],d[M],T,r,c,cnt=0,s[M];
char C[M];
bool check(int x,int y)
{
if(x<0 || y<0 || x>=n || y>=m) return 0;
return 1;
} void dfs(int x,int y)
{
if(!check(x,y)) return ;
if(d[gt(x,y)]) return ;
d[gt(x,y)]=cnt;
dfs(x-r,y); dfs(x+r,y); dfs(x,y-c); dfs(x,y+c);
} int main()
{
scanf("%d",&T);
for(;T;T--)
{
scanf("%d%d%d%d",&n,&m,&r,&c);
int B=1; cnt=0;
for(int i=0;i<n;i++)
for(int j=0;j<m;j++) if(!d[gt(i,j)])
{
cnt++; dfs(i,j);
}
for(int i=0;i<n;i++)
{
scanf("\n%s",C);
for(int j=0;j<m;j++) a[gt(i,j)]=C[j]-'0';
}
for(int i=0;i<n;i++)
{
scanf("\n%s",C);
for(int j=0;j<m;j++)
{
k=C[j]-'0';
if(k!=a[gt(i,j)]) s[d[gt(i,j)]]^=1;
}
}
for(int i=1;i<=cnt;i++) if(s[i]) B=0;
if(B) printf("Yes\n");
else printf("No\n");
for(int i=0;i<=n*m;i++) a[i]=b[i]=d[i]=s[i]=0;
}
}

C

直接建图跑最短路

注意注意注意 :如果在某一车站有两辆公交车同事进站,他们之间要连双向边!


#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#define gt(x,y) ((x)*m+y)
using namespace std; const int M = 2100000;
int t,m,k,b[M],d[M],T,cnt=0,n1,n2,s,t1,t2,head[M],nex[M*3],edge[M*3],ver[M*3],S;
priority_queue<pair<int,int> >q;
struct vv
{
int a,b;
} v[M]; bool cmp(vv a,vv b){return a.a<b.a;} void add(int x,int y,int z)
{
ver[++cnt]=y, nex[cnt]=head[x], head[x]=cnt; edge[cnt]=z;
} void dj()
{
for(int i=0;i<=t;i++) d[i]=0x3f3f3f3f, b[i]=0;
while(q.size()) q.pop(); d[S]=0;
q.push(make_pair(0,S));
while(q.size())
{
while(q.size() &&(b[q.top().second])) q.pop();
if(!q.size()) break;
int x=q.top().second; q.pop(); b[x]=1;
for(int i=head[x];i;i=nex[i])
{
int t=ver[i];
if(d[x]+edge[i]<d[t])
{
d[t]=d[x]+edge[i];
q.push(make_pair(-d[t],t));
}
}
}
} int main()
{
scanf("%d",&T);
for(;T;T--)
{
scanf("%d%d%d%d%d%d",&m,&n1,&n2,&s,&t1,&t2); t=m*(n1+n2)+2; S=t-1; cnt=0;
for(int i=0;i<m;i++)
{
for(int j=0;j<n1+n2;j++) scanf("%d",&v[j].a), v[j].b=j;
sort(v,v+n1+n2,cmp);
for(int j=0;j<n1+n2-1;j++)
{
add(gt(v[j].b,i),gt(v[j+1].b,i),v[j+1].a-v[j].a);
if(v[j+1].a==v[j].a) add(gt(v[j+1].b,i),gt(v[j].b,i),v[j+1].a-v[j].a);
}
if(i)
{
for(int j=0;j<n1;j++) add(gt(j,i-1),gt(j,i),0);
for(int j=0;j<n2;j++) add(gt(n1+j,i),gt(n1+j,i-1),0);
}
if(i!=s-1) continue;
for(int j=0;j<n1+n2;j++) if(v[j].a>=t1 && v[j].a<=t2) add(S,gt(v[j].b,i),v[j].a-t1);
for(int j=0;j<n1+n2;j++) if(v[j].a>=t1 && v[j].a<=t2) add(gt(v[j].b,i),t,t2-v[j].a);
}
dj();
printf("%d\n",min(d[t],t2-t1));
for(int i=0;i<=cnt;i++) edge[i]=ver[i]=nex[i]=0;
for(int i=0;i<=t;i++) head[i]=0;
}
}

D

异或和为n说明这两个数只有在n为1的位上值不同

也就是说如果n有第i位为1,那么两个数有且仅有一个第i位为1; 否则两个数这一位的取值相同

然后两数之差又不能超过m,就可以求出最大的n的子集记为g(窝就是这里求错了超时了15s来着

然后就可以数位dp辣!

$f[i][A][B][G]$表示二进制下的第i位,x是不是a的上限,y是不是b的上限,xy中的哪个或者都没有达到g的上限

由于没想到什么妙妙的转移方法,我选择手写转移


#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define LL long long
using namespace std; const int M = 100001;
int T;
LL n,m,k,f[2][2][2][4],a,b,w,g;
int main()
{
scanf("%d",&T);
for(;T;T--)
{
scanf("%lld%lld%lld%lld",&a,&b,&n,&m); g=0;
LL e=0,S=0;
for(int i=60;i>=0;i--) if((1ll<<i)&n)
{
S+=(1ll<<i);
if(e+(1ll<<i)-(n-S)<=m) g+=(1ll<<i), e+=(1ll<<i);
else e-=(1ll<<i);
}
if(n-g>g || e<0)
{
printf("0\n");
continue;
}
for(int i=60;i>=0;i--) if(g&(1ll<<i))
{
w=(1ll<<i);
break;
} memset(f,0,sizeof(f));
f[1][1][1][0]=1ll;
for(int i=60;i>=0;i--)
{
int t=i&1;
memset(f[t],0,sizeof(f[t]));
LL r=1ll<<i;
for(int A=0;A<=1;A++)
for(int B=0;B<=1;B++)
for(int k=0;k<=2;k++) if(f[!t][A][B][k])
{
if(n&r)
{
if(w==r)
{
int A1=0,B1=0,k1=0;
if((A && (a&r))||(!A))
{
int B1=0;
if(B && (b&r)==0) B1=1;
f[t][A][B1][1]=f[t][A][B1][1]+f[!t][A][B][k];
}
if((B && (b&r))||(!B))
{
int A1=0;
if(A && (a&r)==0) A1=1;
f[t][A1][B][2]=f[t][A1][B][2]+f[!t][A][B][k];
}
}
else
{
if(((k==1 && (g&r))||k!=1) &&( (A &&(a&r))|| !A))
{
int k1=0;
if(k==2 && ((g&r)==0)) k1=2;
if(k==1 && (g&r)) k1=1;
int B1=0;
if(B && (b&r)==0) B1=1;
f[t][A][B1][k1]=f[t][A][B1][k1]+f[!t][A][B][k];
}
if(((k==2 && (g&r))||k!=2) &&( (B &&(b&r))|| !B))
{
int k1=0;
if(k==1 && ((g&r)==0)) k1=1;
if(k==2 && (g&r)) k1=2;
int A1=0;
if(A && (a&r)==0) A1=1;
f[t][A1][B][k1]=f[t][A1][B][k1]+f[!t][A][B][k];
}
}
}
else
{
int A1=0, B1=0;
if(A && (a&r)==0) A1=1;
if(B && (b&r)==0) B1=1;
f[t][A1][B1][k]=f[t][A1][B1][k]+f[!t][A][B][k]; if((A && (a&r)==0) ||(B && (b&r)==0)) continue; A1=0, B1=0;
if(A && (a&r)) A1=1;
if(B && (b&r)) B1=1;
f[t][A1][B1][k]=f[t][A1][B1][k]+f[!t][A][B][k];
}
}
}
LL res=0;
for(int i=0;i<=1;i++)
for(int j=0;j<=1;j++)
for(int k=0;k<=2;k++) res=res+f[0][i][j][k];
printf("%lld\n",res);
}
}

Comet OJ - Contest #12的更多相关文章

  1. Comet OJ - Contest #12 D

    题目描述 求(x,y)的对数满足x∈[0,a],y∈[0,b],x⊕y=0且|x-y|<=m 题解 一种比较sb的做法是考虑x-y的借位,根据借位以及差值进行转移 还有一种比较正常的做法,假设一 ...

  2. Comet OJ - Contest #2 简要题解

    Comet OJ - Contest #2 简要题解 cometoj A 模拟,复杂度是对数级的. code B 易知\(p\in[l,r]\),且最终的利润关于\(p\)的表达式为\(\frac{( ...

  3. Comet OJ - Contest #2简要题解

    Comet OJ - Contest #2简要题解 前言: 我没有小裙子,我太菜了. A 因自过去而至的残响起舞 https://www.cometoj.com/contest/37/problem/ ...

  4. Comet OJ - Contest #4--前缀和

    原题:Comet OJ - Contest #4-B https://www.cometoj.com/contest/39/problem/B?problem_id=1577传送门 一开始就想着暴力打 ...

  5. Comet OJ - Contest #11 题解&赛后总结

    Solution of Comet OJ - Contest #11 A.eon -Problem designed by Starria- 在模 10 意义下,答案变为最大数的最低位(即原数数位的最 ...

  6. Comet OJ - Contest #8

    Comet OJ - Contest #8 传送门 A.杀手皇后 签到. Code #include <bits/stdc++.h> using namespace std; typede ...

  7. Comet OJ - Contest #13-C2

    Comet OJ - Contest #13-C2 C2-佛御石之钵 -不碎的意志-」(困难版) 又是一道并查集.最近做过的并查集的题貌似蛮多的. 思路 首先考虑,每次处理矩形只考虑从0变成1的点.这 ...

  8. Comet OJ - Contest #13 「火鼠的皮衣 -不焦躁的内心-」

    来源:Comet OJ - Contest #13 芝士相关: 复平面在信息学奥赛中的应用[雾 其实是道 sb 题??? 发现原式貌似十分可二项式定理,然后发现确实如此 我们把 \(a^i\) 替换成 ...

  9. Comet OJ - Contest #13 「佛御石之钵 -不碎的意志-」(hard)

    来源:Comet OJ - Contest #13 一眼并查集,然后发现这题 tmd 要卡常数的说卧槽... 发现这里又要用并查集跳过访问点,又要用并查集维护联通块,于是开俩并查集分别维护就好了 一开 ...

随机推荐

  1. ROT13加密和解密

    问题 ROT13(回转13位)是一种简易的替换式密码算法.它是一种在英文网络论坛用作隐藏八卦.妙句.谜题解答以及某些脏话的工具,目的是逃过版主或管理员的匆匆一瞥.ROT13 也是过去在古罗马开发的凯撒 ...

  2. ES6——函数-箭头函数

    箭头函数: 1.普通函数 function 函数名(){...} 2.箭头函数 注意:  1)如果只有一个返回值,{}return可以省略: let arr = [12,5,8,99,33,14,26 ...

  3. smbd - 向客户提供SMB/CIFS服务的服务器

    总览 SYNOPSIS smbd [-D] [-F] [-S] [-i] [-h] [-V] [-b] [-d <debug level>] [-l <log directory&g ...

  4. Linux学习笔记之文件与文件系统的压缩与打包

    四.文件与文件系统的压缩与打包 用途: 传输时减少带宽 Linux系统常见的压缩命令 Compass(旧版压缩,并且在centOS中默认为不安装) gzip, zcat bzip2, bzcat zi ...

  5. 二、SQL Server 分页

    一.SQL Server 分页 --top not in方式 select top 条数 * from tablename where Id not in (select top 条数*页数 Id f ...

  6. Codeforces Round #554 (Div. 2) C.Neko does Maths (gcd的运用)

    题目链接:https://codeforces.com/contest/1152/problem/C 题目大意:给定两个正整数a,b,其中(1<=a,b<=1e9),求一个正整数k(0&l ...

  7. Sass函数:数字函数-percentage()

    1.percentage() percentage()函数主要是将一个不带单位的数字转换成百分比形式: >> percentage(.2) 20% >> percentage( ...

  8. [AtCoder] NIKKEI Programming Contest 2019 (暂缺F)

    [AtCoder] NIKKEI Programming Contest 2019   本来看见这一场的排名的画风比较正常就来补一下题,但是完全没有发现后两题的AC人数远少于我补的上一份AtCoder ...

  9. BZOJ3331 BZOJ2013 压力

    考前挣扎 圆方树这么早就出现了嘛... 要求每个点必须被经过的次数 所以就是路径上的割点/端点++ 由于圆方树上所有非叶子圆点都是割点 所以就是树上差分就可以辣. 实现的时候出了一点小问题. 就是这里 ...

  10. 人生苦短_我用Python_OS对目录/文件操作_005

    # coding=utf-8 import os # 操作文件和目录 ", os.getcwd()) # 获取当前文件的目录 ", os.path.realpath(__file_ ...