问题 E: 张三 (fwkt)

时间限制: 1 Sec 内存限制: 128 MB

题目描述
听说罗翔老师最近很火。
张三今天迫不及待想犯罪。就决定是爆炸罪了，正好罗老师讲了这个案例。
他一眼就看上了 CaPeF_Yyx 的家 （或许是他太 duliu 了） 。
作为张三，她蛮不讲理，所以他的爆炸会波及部分街道。每次爆炸会造成一定的损坏。
作为张三，她出人意料，所以他的爆炸波及的范围每次都不一样。
作为张三，她十恶不赦，所以他的爆炸会进行m 次 （CaPeF_Yyx 的家也太坚固了吧）
作为张三，她井井有条，所以他的爆炸波及的范围是一个闭区间。
作为张三，她精益求精，所以他的爆炸以等差数列的方式造成伤害。
CaPeF_Yyx 家在偏远角落，街道编号从1~n 。
CaPeF_Yyx 街道刚翻新，每间房屋破坏程度都为0。
罗翔老师可以预测每次爆炸的范围。会给到 CaPeF_Yyx 形如[l,r] 的闭区间。
罗翔老师可以预测每次爆炸的伤害。会给到 CaPeF_Yyx 形如l,r,bg,ed 的 4 个数。
表示形如{ql,ql+1,…,qr-1,qr},ql=bg,qr=ed的等差数列，对于街道ai(l≤i≤r)，造成qi的伤害。
CaPeF_Yyx 找到了 Rainy7 修复街道。可是她必须先知道街道到底损坏了多少。
她还忙着拯救世界呢，CaPeF_Yyx 被吓得不轻，于是只能让你帮助她了。
出题人不想太 duliu ，她让你只要输出所有房屋损害程度的最大值和异或和。
Q：话说为什么不直接用魔法阻止张三？
A：因为 唯我法外狂徒张三永世长存 !!!
输入
第一行，两个整数n,m表示房屋的数量，和爆炸次数。

接下来m行，每行四个整数l,r,bg,ed，分别表示等差数列的首项标号，末项标号，首项值，末项值。
输出
一行，两个数，分别表示所有房屋损害程度的最大值和异或和。
样例输入
【样例1】
5 2
1 5 2 10
2 4 1 1
【样例2】
6 2
1 5 2 10
2 4 1 1
样例输出
【样例1】
3 10
【样例2】
3 10
提示
样例1解释:
第一次爆炸产生的伤害：[2,4,6,8,10]。
第二次爆炸产生的伤害：[0,1,1,1,0]。
所有爆炸结束后每个房屋的损伤程度：[2,5,7,9,10]。
输出异或和：

输出最大值：


解题思路：
$可以进行二次差分。设原数组是a[N]，a的差分是d1[N]，d1[N]的差分是d2[N]。$
$每次询问输入l,r,bg,ed，公差d=（ed-bg）/（l-r）。$

先分析d1:
$d1[l]=a[l]-a[l-1],a[l]+=bg,则d1[l]+=bg;$
$d1[l+1]=a[l+1]-a[l],a[l+1]+=bg+d，a[l]+=bg,d1[l+1]+=d;$
$d1[r+1]=a[r+1]-a[r],d1[r+1]+=-((r-l)\ast d+bg)=-ed$
$而且，d1[i]+=d,l<i<=r.$

再分析d2:
$d2[l]+=d1[l]-d1[l-1],d2[l]+=bg;$
$d2[i]不变，l<i<=r;$
$d2[r+1]=d1[r+1]-d1[r],d2[r+1]+=(l-r-1)\ast d-bg;$
$d2[r+2]=d1[r+2]-d1[r+1],d2[r+2]+=(r-l)\ast d+bg=ed。$

我们在写的时候直接只设d2数组就行了，两次求前缀和，就得出a来了。

#include <bits/stdc++.h>
#define ll long long
#define local
using namespace std;
const int N = 1e7+5;
const int inf = 0x3f3f3f3f;//1061109567
ll n,m,l,r,bg,ed,maxx=-1e18,sum;
ll d2[N];
void add(ll l,ll r,ll bg,ll d)
{
    d2[l]+=bg;
    d2[l+1]+=(d-bg);
    d2[r+1]-=(r-l+1)*d+bg;
    d2[r+2]-=(l-r)*d-bg;
}
void preSum()
{
    for(int i=1;i<=n;i++)
        d2[i]+=d2[i-1];
}
int main()
{
#ifdef local
    freopen("input.txt","r",stdin);
#endif // local
     scanf("%lld%lld",&n,&m);
     while(m--)
     {
         scanf("%lld%lld%lld%lld",&l,&r,&bg,&ed);
         ll q = (ed-bg)/(r-l);
         ll cnt=0;
         add(l,r,bg,q);
     }
     preSum();
     preSum();
     for(int i=1;i<=n;i++)
     {
         sum^=d2[i];
         maxx=max(maxx,d2[i]);
     }
     printf("%lld %lld\n",sum,maxx);
     return 0;
}