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+;
 long long num[][N];
 long long ans;
 long long a[N],v;
 int n,m,p;
 long long tmp[];
 int main(){
   while(~scanf("%d%d",&n,&m)){
       memset(num,,sizeof(num));
       memset(a,,sizeof(a));
       memset(tmp,,sizeof(tmp));
       for(int i=;i<=n;i++){
         scanf("%I64d",&a[i]);
         for(int j=;j<=;j++){
             num[j][i]=(a[i]&);
             a[i]/=;
         }
       }
        ans=;
       for(int j=;j<=;j++){
             long long temp=,tp=;
          for(int i=;i<=n;i++){
             if(num[j][i]==){
                 tp+=temp*(temp+)/;
                 temp=;
             }
             else{
                 temp++;
             }
             if(i==n){
                 tp+=temp*(temp+)/;
                 temp=;
             }
          }
         ans+=(tp<<j);
       }
       for(int i=;i<=m;i++){
         scanf("%d%I64d",&p,&v);
          for(int j=;j<=;j++){
              tmp[j]=(v&);
               v/=;
         }
         for(int j=;j<=;j++){
             long long sum=;
             int tt=p;
             if(tmp[j]==num[j][p])continue;
             long long lnum=,rnum=;
             while(num[j][--tt])
                 lnum++;
                 tt=p;
             while(num[j][++tt])
                 rnum++;
             if(tmp[j]==&&num[j][p]==){
                 sum-=(lnum+rnum+)*(lnum+rnum+)/;
                 sum+=(lnum+)*lnum/;
                 sum+=(rnum+)*rnum/;
                 num[j][p]=;
             }
             else{
                 sum-=(lnum+)*lnum/;
                 sum-=(rnum+)*rnum/;
                 sum+=(lnum+rnum+)*(lnum+rnum+)/;
                 num[j][p]=;
              }
            ans+=(sum<<j);
         }
          printf("%I64d\n",ans);
       }
   }

 }