#415 Div2 Problem C Do you want a data? (math && 前后缀和 && 快速幂)

题意: 首先定义集合的F值为  这个集合里面最大值和最小值的差。 现给出一个拥有n个数的集合(没有相同的元素), 要求求出这个集合内所有子集的F的值的和。例如: {4、7}这个集合里面有子集{4}、{7}、{4, 7}, 则这些子集的F值分别为4-4=0、7-7=0、7-4=3, 所以最后的结果就是0+0+3 = 3!

以下分析引用至 : http://blog.csdn.net/dragon60066/article/details/72599167

分析: 不难想到要先使数组升序方便计算和思考, 然后观察发现如果选定任意一头一尾两个元素, 假如为ai、aj (ai<aj) 那这个两个元素之间围起来的区间可以看成ai....aj这样的数集, 且这集合里面所有子集贡献的F值都等于 aj - ai, 那这个区域包含了多少个这样的F值的集合呢?用笔简单列出表后不难发现规律, 这个和为个(需要用到快速幂), 至于具体为什么, 可以看自己多加思考一下且在参考博客有说明。当时想到这个后马上打出了枚举程序

#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int mod = 1e9+7;
LL arr[300005];
LL mi[300005];
LL qmod(LL a, LL b, LL c)
{
    LL ans = 1;
    while(b){
        if(b&1) ans = (ans*a)%c;
        a = (a*a)%c;
        b >>= 1;
    }
    return ans;
}
int main(void)
{
    int n;
    scanf("%d", &n);
    for(int i=0; i<n; i++){
        scanf("%lld", &arr[i]);
    }
    sort(arr, arr+n);
    LL sum = 0;
    LL temp;
    for(int i=0; i<n; i++){
        mi[i] = qmod(2, i, mod);
    }
    for(int i=0; i<n-1; i++){
        for(int j=n-1; j>i; j--){
            LL sub = arr[j] - arr[i];
            if(j-i-1!=0)temp = mi[j-i-1];
            else temp = 1;
            temp%=mod;
            sum += (((temp%mod)*(sub%mod))%mod)%mod;
            sum%=mod;
        }
    }
    printf("%lld\n", sum);
    return 0;
}

然后就可以顺利超时了, 因为像我这样做的复杂度是O(n^2)!需要优化!于是上网搜索看到了前缀和的解法, 甚是巧妙, 而且貌似榜上大佬大多也是这样做!

引用一下——>进一步思考,对于每个i, j,取值可能有,第一种前面系数是(a[2]-a[1]+a[3]-a[2]+...+a[n]-a[n-1])=a[n]-a[1],同理可以发现第二种是a[n]+a[n-1]-a[1]-a[2]......,那么就推出来了:

那这个直观一点的话到底是在干什么玩意呢?看下面根据描述打出来的表

看到规律了吧!发现系数是对称的, 那我们根据系数的规律枚举从0到n-2的2次方就行了!不过细想的话, 虽然是找到了系数的规律, 但是好像还是不好实现, 参考别人的程序看到别人使用了前后缀和的做法, 巧妙的完成了依据上面规律进行的枚举操作, 具体的话不好说, 看一下程序便知!

以下为拙略表达, 可以跳过: 不过我还是说一下, 按我的理解就是从前面到中间的系数(以n=6为例, 枚举从2^0到2^2的系数)倒是不难枚举出来, 但是后面的就略微麻烦了, 解决方法——>定义前后缀和(例:sumf[i]是从1到i的前缀和, sumb[i]是n-1到i的后缀和), 那还是以n=6为例, 2^0到2^2的系数分别可以用sumb[6] - sum[1] 、 sumb[5]-sumf[2]、 sumb[4]-sumf[3]来解决, 后面的2^3 可以 sumb[3] - sumf[4], 停!有没有发现 计算这个差的过程中实际就是 a6+a5+a4+a4-a4-a3-a1-a2! a3和a4巧妙的被约去了, 2^4的计算也是同样道理!类似于回文?好思想!

#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int mod = 1e9+7;
const int maxn = 300005;
LL arr[maxn], sumf[maxn], sumb[maxn], mi[maxn];
LL quick_mod(LL a, LL b, LL c)
{
    LL ans = 1;
    while(b){
        if(b&1) ans = (ans*a)%c;
        a = (a*a)%c;
        b >>= 1;
    }
    return ans;
}
int main(void)
{
    int n;
    scanf("%d", &n);
    for(int i=1; i<=n; i++)
        scanf("%d", &arr[i]);
    sort(arr+1, arr+1+n);///使序列单调
    sumf[0] = sumb[n+1] = 0;
    mi[0] = 1;
    for(LL i=1; i<n-1; i++)
        mi[i] = quick_mod(2, i, mod);///二的次方数组
    for(int i=1; i<=n; i++)
        sumf[i] = sumf[i-1] + arr[i];///前缀和
    for(int i=n; i>=1; i--)
        sumb[i] = sumb[i+1] + arr[i];///后缀和
    LL ans = 0;
    for(int i=1; i<n; i++){ //实际就是0 -- n-2 次循环
        LL temp = (-sumf[i] + sumb[n-i+1])%mod;
        ans = (ans%mod + (temp * mi[i-1])%mod)%mod;
    }
    printf("%lld\n", ans);
    return 0;
}

时间: 2024-11-03 21:58:26

#415 Div2 Problem C Do you want a data? (math && 前后缀和 && 快速幂)的相关文章

dutacm.club Water Problem(矩阵快速幂)

Water Problem Time Limit:3000/1000 MS (Java/Others)   Memory Limit:163840/131072 KB (Java/Others)Total Submissions:1228   Accepted:121 [Submit][Status][Discuss] Description 函数 f:Z+→Z .已知 f(1),f(2) 的值,且对于任意 x>1,有 f(x+1)=f(x)+f(x?1)+sin(πx2) . 求 f(n) 的

A Simple Math Problem(矩阵快速幂)(寒假闭关第一题,有点曲折啊)

A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 155 Accepted Submission(s): 110   Problem Description Lele now is thinking about a simple function f(x). If x < 10 f(x) = x.If

hdu 1757 A Simple Math Problem 矩阵快速幂

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1757 Lele now is thinking about a simple function f(x).If x < 10 f(x) = x.If x >= 10 f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + …… + a9 * f(x-10);And ai(0<=i<=9) can only be 0 or 1 .Now, I w

fzu Problem 2198 快来快来数一数 (快速幂+优化)

题目链接: Problem  2198  快来快来数一数 题目描述: 给出n个六边形排成一排,a[i]代表i个六边形能组成的生成树个数,设定s[i]等于a[1]+a[2]+a[3]+....+a[i-1]+a[i],问s[n]为多少? 解题思路: n取值范围[1, 1018],打表内存不够,然后就要考虑快速幂咯!纳尼!!!!快速幂写出来竟然超时,敢信?果然还是见题太少了.(GG) 对于a[n] = 6*a[n-1] - a[n-2],可以很明显看出. 然后求和的时候就要化简一番了,但是并不是很难

LightOJ 1070 Algebraic Problem (推导+矩阵快速幂)

题目链接:LightOJ 1070 Algebraic Problem 题意:已知a+b和ab的值求a^n+b^n.结果模2^64. 思路: 1.找递推式 得到递推式之后就是矩阵快速幂了 注意:模2^64,定义成unsigned long long 类型,因为无符号类型超过最大范围的数与该数%最大范围 的效果是一样的. AC代码: #include<stdio.h> #include<string.h> #define LL unsigned long long struct Ma

BestCoder Round #29——A--GTY&#39;s math problem(快速幂(对数法))、B--GTY&#39;s birthday gift(矩阵快速幂)

GTY's math problem Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 0    Accepted Submission(s): 0 Problem Description GTY is a GodBull who will get an Au in NOI . To have more time to learn alg

HDU 2256 Problem of Precision (矩阵快速幂)

HDU 2256 Problem of Precision (矩阵快速幂) ACM 题目地址:HDU 2256 Problem of Precision 题意: 给出一个式子,求值. 分析: 推起来最后那步会比较难想. 具体过程见: 表示共轭只听说过复数的和图的... 这构题痕迹好明显... 跟基友开玩笑说:如果遇到这种题,推到Xn+Yn*sqrt(6)这步时,打表最多只能打到10就爆int了,这是输出正解和Xn,说不定眼神好能发现ans = Xn * 2 - 1呢.= =... 代码: /*

HDU 5685 Problem A | 快速幂+逆元

Problem A Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 463    Accepted Submission(s): 162 Problem Description 度熊手上有一本字典存储了大量的单词,有一次,他把所有单词组成了一个很长很长的字符串.现在麻烦来了,他忘记了原来的字符串都是什么,神奇的是他竟然记得原来那些字符串的

hdu1757 A Simple Math Problem(矩阵快速幂)

题目: A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3522    Accepted Submission(s): 2130 Problem Description Lele now is thinking about a simple function f(x). If x < 10 f