Inna and Binary Logic

Codeforces Round #234 (Div. 2) E:http://codeforces.com/problemset/problem/400/E

题意：给你n个数，然后每相邻的两个数可以通过and运算生成一个新的数，然后这些新生成的n-1个数每相邻的两个数又通过and运算成n-2个数，最后只会剩下一个数，然后让你求这n(n+1)/2个数的和，然后每一次会更新最底层的某个数，然后操作之后，输出刚才的总和。

题解：这一题，虽然是看题解，然后自己敲出来的，但是还是有点成就感和收获。首先，这一题的思路很巧妙。如果所有的数都是1或者0，加入说序列是1110001，通过计算，发现其实和就是(3+1)*3/2+(1+1)*1/2==7,与连续1的个数有关，加入连续1的个数是x，那么这连续的x个1，形成的和就是(x+1)*x/2;总和就是把所有连续的1和相加。想到这里，就可以知道，数的范围是1e5，最多是2^17，所以可以把每个数拆成17位，每一位要么是0或者是1，这样只要统计底层的连续1有多少就可以了。num[i][j]表示第j个数的第i位，一开始我们可以计算出总和，然后更新时，如果更新数的这一位和原来相同则这一位不用变化，否则，如果是1要0，可以把ans先减去原来连续个一所形成的和，然后加上这个数左边连续1的和以及右边连续1的和，然后把这一位变成0，如果是0变1，则相反。这一要注意数据范围，在过程中有可能爆int，所以要用long long，并且在求和过程中使用的局部变量也要用long long，由于自己没有注意到这样的问题结果wa 2发，int和long long之间转换出了问题。也许，有人会问这样会不会t，首先事实上没有t，而且跑的很快。从理论上讲，也不会，因为要出现很长的连续的1是很难的，必须保证这个连续的数在某些位上都是1，并且连续，很难，这样的数据很难。所以很快。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e5+20;
long long num[18][N];
long long ans;
long long a[N],v;
int n,m,p;
long long tmp[18];
int main(){
  while(~scanf("%d%d",&n,&m)){
      memset(num,0,sizeof(num));
      memset(a,0,sizeof(a));
      memset(tmp,0,sizeof(tmp));
      for(int i=1;i<=n;i++){
        scanf("%I64d",&a[i]);
        for(int j=0;j<=17;j++){
            num[j][i]=(a[i]&1);
            a[i]/=2;
        }
      }
       ans=0;
      for(int j=0;j<=17;j++){
            long long temp=0,tp=0;
         for(int i=1;i<=n;i++){
            if(num[j][i]==0){
                tp+=temp*(temp+1)/2;
                temp=0;
            }
            else{
                temp++;
            }
            if(i==n){
                tp+=temp*(temp+1)/2;
                temp=0;
            }
         }
        ans+=(tp<<j);
      }
      for(int i=1;i<=m;i++){
        scanf("%d%I64d",&p,&v);
         for(int j=0;j<=17;j++){
             tmp[j]=(v&1);
              v/=2;
        }
        for(int j=0;j<=17;j++){
            long long sum=0;
            int tt=p;
            if(tmp[j]==num[j][p])continue;
            long long lnum=0,rnum=0;
            while(num[j][--tt])
                lnum++;
                tt=p;
            while(num[j][++tt])
                rnum++;
            if(tmp[j]==0&&num[j][p]==1){
                sum-=(lnum+rnum+1)*(lnum+rnum+2)/2;
                sum+=(lnum+1)*lnum/2;
                sum+=(rnum+1)*rnum/2;
                num[j][p]=0;
            }
            else{
                sum-=(lnum+1)*lnum/2;
                sum-=(rnum+1)*rnum/2;
                sum+=(lnum+rnum+1)*(lnum+rnum+2)/2;
                num[j][p]=1;
             }
           ans+=(sum<<j);
        }
         printf("%I64d
",ans);
      }
  }

}