BZOJ3652 大新闻

首先我们打一个i^j的表

我们发现，第一，这是一个关于对角线对称的（废话）

第二，对于任意一个从左上角开始，边长为2的幂的正方形（比如图中用白框框起来的），其下边和右边的边长为2的幂的正方形就是把这个正方形的每个元素加上这个2的幂

第三，对于任意一个从左上角开始，边长为2的幂的正方形，其右下的边长为2的幂的正方形和这个正方形是一样的

第四，对于从顶部开始，向下长度为2的幂的一列，0到2的幂-1的所有数都恰好出现一次

发现了这些规律之后我们看题，当p=1的时候，相当于求从第一列开始n列每列的最大值的和/n，p=0的时候，相当于求从左上角开始n*n的正方形内所有元素的和/(n*n)

当0<p<1时，只要把上边两个值加权加起来即可

先说如何求n*n的正方形的和

找到小于等于n的最大的2的幂的正方形，左上角的正方形可由规律4直接求出，右上的长方形可由规律2和规律4直接求出，由规律1左下的长方形和右上的长方形和是一样的，然后由规律3，我们可以递归求右下的正方形，这样复杂度是log的

在说如何求n列里每列最大值的和

找到小于n的最大的2的幂的正方形，由规律2和3，可知最大值一定在左下和右上，由规律4右上部分的可以直接求，然后我们递归左下的长方形

这样递归就有了3个变量x,y,u，x代表有多少行，y代表有多少列，u代表又规律2现在这个递归部分被加了多少值

在进行一次递归之后，y一定是2的幂

找到小于y的最大的2的幂xx，这里分两种情况讨论，若x>xx，即有左下部分，则算出右上然后递归左下，若x<=xx，即没有左下部分，由规律2我们递归左半边，然后右半边可以根据左半边的直接算出来

然后就好了

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<set>
#include<bitset>
#include<map>
using namespace std;
#define MAXN 1010
#define MAXM 1010
#define ll long long
#define eps 1e-8
#define INF 1000000000
#define MOD 1000000007
ll n;
double p;
ll low(ll x){
    ll t=1;
    for(;t<=x;t<<=1){
          
    }
    return t>>1;
}
double sum(ll x,ll y){
    return 1.0*(x+y)*(y-x+1)/2;
}
double cal2(ll x){
    if(x<=1){
        return 0;
    }
    ll xx=low(x);
    return sum(0,xx-1)*xx+sum(xx,xx+xx-1)*(x-xx)*2+cal2(x-xx);
}
double cal1(ll x,ll y,ll u){
    if(!x||!y){
        return 0;
    }
    if(y==1){
        return u;
    }
    ll xx=low(y-1);
    if(x>xx){
        return (1.0*u+xx*2-1)*(1.0*y-xx)+cal1(x-xx,xx,u+xx);
    }
    return cal1(x,xx,u)*2+1.0*xx*xx;
}
int main(){
    scanf("%lld%lf",&n,&p);
    printf("%.6lf\n",p*cal1(n,n,0)/n+(1.0-p)*cal2(n)/n/n);
    return 0;
}
  
/*
  
*/