注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

时光机TimeMachine

——一个退役OIer

 
 
 

日志

 
 

1923: [Sdoi2010]外星千足虫  

2014-11-11 17:01:09|  分类: Problems |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

1923: [Sdoi2010]外星千足虫

NOIP被虐爆了各种烦躁……
刷道题冷静一下
题解在代码下方

#include<cstdio>
#include<cstdlib>
#include<bitset>
#include<algorithm>
using namespace std;
const int N = 1100;
typedef bitset<N> Matrix[N*2];
int n,m;
Matrix mat;
int scan(){int i=0;scanf("%d",&i);return i;}
int gauss(Matrix a){
int re=0,i,j,k;
for(i=1;i<=n;i++){
int t = i;
for(j=i;j<=m;j++)
if(a[j][i]){t = j;break;}
re = max(re,t);//update ans
if(!a[j][i])return -i;

if(t!=i)swap(a[i],a[t]);
for(j=i+1;j<=m;j++)
if(a[j][i])a[j] = a[j] ^ a[i];
}
//back
for(i=n;i;i--){
for(j=i+1;j<=n;j++)
if(a[i][j])a[i][n+1] = a[j][n+1] ^ a[i][n+1];
}
return re;
}
void print(){
int i,j;
for(i=1;i<=m;i++){
for(j=1;j<=n+1;j++)printf("%d",(int)mat[i][j]);
printf("\n");
}printf("\n");
}
int main(){
int i,j,a;char cc;
//freopen("in.txt","r",stdin);
n = scan();m = scan();
for(i=1;i<=m;i++){
cc=' ';
while(cc==' '||cc=='\r'||cc=='\n')cc=getchar();
for(j=1;j<=n;j++){
if(cc=='0')mat[i][j] = 0;
else mat[i][j] = 1;
cc=getchar();
}
scanf("%d",&a);
mat[i][n+1] = a;
}
a = gauss(mat);
if(a<0){

printf("Cannot Determine\n");
}else{
printf("%d\n",a);
for(i=1;i<=n;i++){
if(mat[i][n+1])printf("?y7M#\n");
else printf("Earth\n");
}
}
return 0;
}

题解:
高斯消元解XOR方程组
比较裸……但是第一问要怎么搞?我一开始以为要二分答案什么的= =、
但是实际上在高斯消元交换行的时候每次用最下面用到的行的行数更新答案就行了
然后据说要压位,实际上直接用bitset秒就行了(STL大法好)
  评论这张
 
阅读(283)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017