【loli的胡策】NOIP训练7.17（模拟+乱搞二进制+乱搞BFS）

T2题目：
在美丽的比特镇一共有n 个景区，编号依次为1 到n，它们之间通过若干条双向道路连接。 Byteasar 慕名来到了比特镇旅游，不过由于昂贵的门票费，他只能负担起4 个景区的门票费。他可以在任意景区开始游览，然后结束在任意景区。Byteasar 的旅游习惯比较特殊，一旦他路过了一个景区，他就一定会进去参观，并且他永远不会参观同一个景区两次。所以他想知道，有多少种可行的旅游路线，使得他可以恰好参观4 个景区呢？即，有多少条简单路径恰好经过了4 个点。（n<=1500)
样例输入：
4
0101
1010
0101
1010
输出：8
题解：
假设路径是a-b-c-d，考虑枚举中间这条边b-c ，计算有多少可行的a-d.设dx表示点x 的度数，那么边b-c对答案的贡献为（db-1）*（dc-1）-经过bc 这条边的三元环个数。
三元环个数怎么求呢？
我们可以用bitset，这个东西在我理解就是一个强转二进制的数组，所以读入这种010101啥的简直小菜
考虑三元环的特征：i，j都和一个点相连，某一位应该都是1，我们可以用&，同1为1，不然为0，然后统计1的个数就是形成三元环的个数，唔，这么大的数字要用longlong啊
一些奇妙的操作还是要多多了解啊
代码：

#include <cstdio>
#include <bitset>
#include <cstring>
#define LL long long
using namespace std;
int n,i,j;LL ans=0;
int dian[1505];
bitset<1505> a[1505];
int main()
{
    freopen("tour.in","r",stdin);
    freopen("tour.out","w",stdout);
    scanf("%d",&n);
    char c;
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++)
      {
        c=getchar();
        while (c!='0'&& c!='1') c=getchar();
        a[i][j]=c-'0';
        if (a[i][j]) dian[i]++;
      }

    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++)
        if (i!=j && a[i][j])
        {
            ans+=(dian[i]-1)*(dian[j]-1);
            ans-=(a[i]&a[j]).count();
          }
    printf("%lld",ans);
}

T3题目：
给一张n个点的无向图，点i的权值是vali，如果vali & valj=valj 的话，从i向j连边，i，j之间还连了m条边，每条边权值为1，求1到所有点的最短路
n<=20000 m<=300000 val[i]<2^20
题解：
40分就是送分
特殊边的处理是非常恶心的，因为我们不能枚举判定，也不能直接相连，所以需要考虑转化。
考虑新增 2^20个点，这些点中 i 向它所有位为1的位置去掉这个1的影响的点连一条边权为0的有向边.
例如1111->1110（只相差一位1），我们可以通过走多条这种有向边的方式到达他的所有子集点，即i&j=i的点。对于原来的 n 个点，先把 m条边连好，连权值为1的有向边，然后对于 i号点，由它向新增的第 vali 个点（就是他权值所对应的点）连一条权值为 1 的有向边，再由新增的第 vali 个点向它连一条权值为 0的有向边。
BFS 的时候，每次要把用 0 权值的边连接的所有点都加入队尾，以保证距离不降。

代码：

#include <cstdio>
#include <queue>
#include <cstring>
#define INF 1e9
#define N 1248578
using namespace std;
queue <int> q;
int val[N+5],dis[N+5],n,qw;
int tot,nxt[N*2+5],point[N*2+5],v[N*2+5],point0[N*2+5];
void addline(int x,int y)
{
    tot++; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;
}
void addline0(int x,int y)
{
    tot++; nxt[tot]=point0[x]; point0[x]=tot; v[tot]=y;
}
void dfs(int x,int len)
{
    int i;
    if (dis[x]>=0) return;
    q.push(x);
    dis[x]=len;
    for (i=point0[x];i;i=nxt[i])
      dfs(v[i],len);
    if (x>=qw) return;
    for (i=0;i<20;i++)
      if (x>>i&1)//如果x的这一位是1 
        dfs(x^(1<<i),len); 
}
int main()
{
    freopen("walk.in","r",stdin);
    freopen("walk.out","w",stdout);
    int i,j,m,x,y;
    scanf("%d%d",&n,&m);
    qw=1<<20; 
    for (i=1;i<=n;i++) 
    {
        scanf("%d",&x);
        addline(qw+i,x);
        addline0(x,qw+i);
    } 
    for (i=1;i<=m;i++) 
    {
        scanf("%d%d",&x,&y); addline(qw+x,qw+y);
    }
    for (i=0;i<=qw+n;i++) dis[i]=-1; 
    dfs(qw+1,0);
    while (!q.empty())
    {
        int now=q.front(); q.pop();
        for (i=point[now];i;i=nxt[i])
          dfs(v[i],dis[now]+1);
    }
    for (i=1;i<=n;i++)
      printf("%d\n",dis[qw+i]);
}