Comet OJ - Contest #12


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的更多相关文章
- Comet OJ - Contest #12 D
题目描述 求(x,y)的对数满足x∈[0,a],y∈[0,b],x⊕y=0且|x-y|<=m 题解 一种比较sb的做法是考虑x-y的借位,根据借位以及差值进行转移 还有一种比较正常的做法,假设一 ...
- Comet OJ - Contest #2 简要题解
Comet OJ - Contest #2 简要题解 cometoj A 模拟,复杂度是对数级的. code B 易知\(p\in[l,r]\),且最终的利润关于\(p\)的表达式为\(\frac{( ...
- Comet OJ - Contest #2简要题解
Comet OJ - Contest #2简要题解 前言: 我没有小裙子,我太菜了. A 因自过去而至的残响起舞 https://www.cometoj.com/contest/37/problem/ ...
- Comet OJ - Contest #4--前缀和
原题:Comet OJ - Contest #4-B https://www.cometoj.com/contest/39/problem/B?problem_id=1577传送门 一开始就想着暴力打 ...
- Comet OJ - Contest #11 题解&赛后总结
Solution of Comet OJ - Contest #11 A.eon -Problem designed by Starria- 在模 10 意义下,答案变为最大数的最低位(即原数数位的最 ...
- Comet OJ - Contest #8
Comet OJ - Contest #8 传送门 A.杀手皇后 签到. Code #include <bits/stdc++.h> using namespace std; typede ...
- Comet OJ - Contest #13-C2
Comet OJ - Contest #13-C2 C2-佛御石之钵 -不碎的意志-」(困难版) 又是一道并查集.最近做过的并查集的题貌似蛮多的. 思路 首先考虑,每次处理矩形只考虑从0变成1的点.这 ...
- Comet OJ - Contest #13 「火鼠的皮衣 -不焦躁的内心-」
来源:Comet OJ - Contest #13 芝士相关: 复平面在信息学奥赛中的应用[雾 其实是道 sb 题??? 发现原式貌似十分可二项式定理,然后发现确实如此 我们把 \(a^i\) 替换成 ...
- Comet OJ - Contest #13 「佛御石之钵 -不碎的意志-」(hard)
来源:Comet OJ - Contest #13 一眼并查集,然后发现这题 tmd 要卡常数的说卧槽... 发现这里又要用并查集跳过访问点,又要用并查集维护联通块,于是开俩并查集分别维护就好了 一开 ...
随机推荐
- 203-基于ARM和双TI DSP TMS320C6678的6UCPCI高清编解码处理平台
基于ARM和双TI DSP TMS320C6678的6UCPCI高清编解码处理平台 1.产品简介 该板卡由我公司自主研发,以TI Cortex-A8.TI 双DSP TMS320C6678为设计核心, ...
- Move Over and Click Link
Move Over and Click Link [Documentation] 等待悬浮菜单中的元素出现并单击元素 [Arguments] ${hover_locator} ${opt_locato ...
- loj2542 「PKUWC2018」随机游走 MinMax 容斥+树上高斯消元+状压 DP
题目传送门 https://loj.ac/problem/2542 题解 肯定一眼 MinMax 容斥吧. 然后问题就转化为,给定一个集合 \(S\),问期望情况下多少步可以走到 \(S\) 中的点. ...
- BZOJ2756 [SCOI2012]奇怪的游戏 最大流
好久没有写博客了.不过这个博客也没有多少人看 最近在写网络流,为了加深理解,来写一两篇题解. 对整个棋盘进行黑白染色以后可以发现,一次操作就是让二分图的两个点的值分别 \(+1\). 这样,我们就可以 ...
- SVN查看所有人的日志提交记录
1. svn默认显示最近一周的文件提交和修改记录,怎么查看更长时间的日志记录呢? 2. TortoiseSVN 3. 点击show all 或者NEXT 100,就可显示更长时间的文件提交记录.
- asp.net core web api 版本控控制
通过微软的一个库Microsoft.AspNetCore.Mvc.Versioning实现asp.net core web api的版本控制. 以两种形式组织了Controller: 文件夹分开 命名 ...
- JAVA学习笔记--方法中的参数调用是引用调用or值调用
文献来源:<JAVA核心技术卷Ⅰ>,第4章第5节 (没有相关书籍的可看传送门) ps:测试demo因为偷懒,用的是String对象 结论:Java使用的是对象的值引用.即将任何对象所在内存 ...
- 如何用命令行启动android的模拟器(简要描述)
参考来源 http://blog.sina.com.cn/s/blog_5033827f0101cxhz.html 环境前置条件: 已安装JDK,已下载并安装了SDK 1.查看已安装的SDK平台:an ...
- 51nod 1122:机器人走方格 V4 (矩阵快速幂)
题目链接 昨天上随机信号分析讲马氏链的时候突然想到这题的解法,今天写一下 定义矩阵A,Ans=A^n,令A[i][j]表示,经过1次变换后,第i个位置上的机器人位于第j个位置的情况数,则Ans[i][ ...
- boost IOStreams
Boost.IOStreams provides numerous implementations of the two concepts. Devices which describes data ...