Codeforces Round #242 (Div. 2) C. Magic Formulas

思路:


Q=q1^q2.......^qn = p1^p2......^pn^((1%1)^....(1%n))^((2%1)^......(2%n))^....


Q的求解过程分成两部分


第一部分是求p1^p2......^pn


第二部分是求((1%1)^....(1%n))^((2%1)^......(2%n))^....


将其化成矩形的形式


1%1  1%2  ........... 1%n


2%1  2%2  ............2%n


.....................................


n%1   n%2 ............. n%n


i小于除数时,即 i%j = i (i<j)


故该矩形的上对角线变成


0    1    1    1    1    1    .....     1              n-1)个1


     0    2    2    2    2    .....      2             n-2)个2


          0    3    3    3    .....      3              (n-3)个3


               0    4    4    .....      4              (n-4)个4


                   ........................

    

                      .....................


                          .................


                             ..............

 

                                ...........


                                     ......



                                        n-1                  1 n-1


                                          0                  n




注意偶数个相同的元素异或结果为0,奇数个相同元素异或结果为本身,0与其他元素异或结果为该元素


故对于矩形上三角 只需要考虑(n-i)的奇偶性故这部分代码简化为


if((n-i)%2 == 0) ans^=i;


现在考虑矩形下三角,将下三角取余数的


0


0    0


0    1     0


0     0    1    0


0     1    2    1    0


0     0     0    2    1    0


0     1     1    3    2    1    0


0     0     2    0    3    2    1     0


0     1     0    1    4    3    2     1    0


...............................................


...................................................




根据规律


第一列为0循环变化


第二列为0~1循环变化


i列为0 ~i-1 循环变化


...........


i列元素个数为    n-i+1


#include <stdio.h>
#include <iostream>
#include <string>
#include <string.h>
#include <algorithm>
#include <stdlib.h>
#include <math.h>
#include <vector>
#include <map>

using namespace std;
const int maxn = 1000005;
int p[maxn];
int yh[maxn];
int main(){
    int ss = 0;
    for(int i = 1;i < maxn;i++){//对1~n的异或和预处理
        ss ^= i;
        yh[i] = ss;
    }
    //for(int i = 1;i < 100;i++)
      //  printf("%d   ",yh[i]);
    int n;
    while(~scanf("%d",&n)){
        int ans = 0;
        for(int i = 1;i <= n;i++){
            scanf("%d",p+i);
            ans ^= p[i];
        }
        
        for(int i = 1;i <= n;i++)//上三角
            if((n-i)%2)
                ans ^= i;
        
        //printf("%d\n",ans);
        int tt = 1;
        for(int i = 2;i <= n;i++){//下三角
            if(n-i <= 0)
                continue;
            if(((n-i)/i)%2)
                ans ^= tt;
            int temp = (n-i)%i;
            ans ^= yh[temp];

            tt ^= i;
        }
        printf("%d\n",ans);
    }
    return 0;
}




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值