100047. 【NOIP2017提高A组模拟7.14】基因变异

Description

21 世纪是生物学的世纪，以遗传与进化为代表的现代生物理论越来越多的 进入了我们的视野。 如同大家所熟知的，基因是遗传因子，它记录了生命的基本构造和性能。 因此生物进化与基因的变异息息相关，考察基因变异的途径对研究生物学有着 至关重要的作用。现在，让我们来看这样一个模型：
1、所有的基因都可以看作一个整数或该整数对应的二进制码；
2、在 1 单位时间内，基因 x 可能会在其某一个二进制位上发生反转；
3、在 1 单位时间内，基因 x 可能会遭到可感染基因库内任一基因 y 的影响 而突变为 x XOR y。
现在给出可感染基因库，Q 组询问，每组给出初始基因与终止基因，请你 分别计算出每种变异最少要花费多少个单位时间。

Input

第 1 行两个整数 N, Q； 第 2 行 N 个用空格隔开的整数分别表示可感染基因库内的基因； 接下来 Q 行每行两个整数 S、T，分别表示初始基因与终止基因。

Output

输出 Q 行，依次表示每组初始基因到终止基因间最少所花时间。

Sample Input

3 3
1 2 3
3 4
1 2
3 9

Sample Output

2
1
2

Data Constraint

对于 20%的数据，N = 0；
对于额外 40%的数据，1 ≤ Q ≤ 100，所有基因表示为不超过 10^4 的非负整 数；
对于 100%的数据，0 ≤ N ≤ 20，1 ≤ Q ≤ 10^5，所有基因表示为不超过 10^6 的 非负整数。

Solution

20%:直接比较起止基因二进制位不同数，即 st Xor ed 中 1 的个数。
40%:bfs。从起点开始bfs到终点统计答案。
100%：我们要求的是st ^ x1 ^ x2 ^ … ^ xk =ed。两边同时异或st，得：x1 ^ x2 ^ … ^ xk=st ^ ed。所以实际上我们只要求出0到st ^ ed 的最少步数就可以了。因而只需要从0开始做一遍bfs，每次询问在f数组里O(1)获取答案。时间复杂度O(Q)。

Code1

Uses math;

var
    bucket:array[-1..5100000] of boolean;
    queue:array[-1..5100000] of longint;
    dis:array[-1..2100000] of longint;
    num:array[-1..100] of longint;
    u,n,m,i,tmp,maxn,head,tail:longint;
begin
    read(n,m);
    for i:=1 to n do read(num[i]);
    for i:=0 to 19 do begin inc(n); num[n]:=1 << i; end;
    head:=0; tail:=1; queue[1]:=0; dis[0]:=0; bucket[0]:=True;
    while(head<tail) do
    begin
        inc(head);
        u:=queue[head];
        for i:=1 to n do
        begin
            tmp:=u xor num[i];
            if (bucket[tmp]=False) then
            begin
                inc(tail); queue[tail]:=tmp;
                dis[tmp]:=dis[u]+1;
      	        bucket[tmp]:=True;
            end;
        end;
    end;
    for i:=1 to m do
    begin
        read(head,tail);
        writeln(dis[head xor tail]);
    end;
end.

Code2

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 120011
using namespace std;
int n,q,s,t,a[30],bz[N*10],d[N*100],f[N*10];
int main(){
	scanf("%d%d",&n,&q);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=20;i++) a[++n]=1<<(i-1);
	int i=0,j=1;
	memset(f,127,sizeof(f));
	f[0]=0;
	while(i++<j){
		int x=d[i];
		for(int l=1;l<=n;l++){
			x^=a[l];
			if(f[d[i]]+1<f[x]){
				f[x]=f[d[i]]+1;
				d[++j]=x;
			}
			x^=a[l];
		}
	}
	while(q--){
		scanf("%d%d",&s,&t);
		printf("%d\n",f[s^t]);
	}
	return 0;
}

---

---

作者：zsjzliziyang
QQ:1634151125
转载及修改请注明
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/85015665