很好的一道数学推导题

45分做法

$O(N^4)$暴力枚举四个材料

55分做法

从第一个约束条件可得到所有可行答案都是单调递增的,所以可以排序一遍,减少枚举量,可以拿到55分

100分做法

首先可以发现每个$x$都小于n,而$n$最大值只是$15000$,所以可以开一个桶来存每个魔法值出现的次数

回忆一下3个约束条件

$xa<xb<xc<xd$ ①

$xb-xa=2(xd-xc)$ ②

$xb-xa<(xc-xb)/3$ ③

现在魔改一下这三个式子

设$t=xd-xc$

所以②可化为$xb-xa=2t$ ④

将④代入③

$2t<(xc-xb)/3$

移项一下,就变成

$6t<xc-xb$ ⑤

再魔改一下

设$6t+k=xc-xb$(就是把差的部分补上去)

于是可以画出来一个图

显然,$A$的最小值为$1$,$D$的最大值为$n$

由图可得$AD=9t+k$

所以我们可以尝试着枚举t,用t来表示各个魔法值的值

由上易得t的范围为$1<=t<=(n-1)/9$

在代码中为了避免除法写成$t*9<n$

再枚举D,因为我们已经枚举出了t,所以C的值是可以直接算出来的

$C=D-t$

又因为使$A,B,C,D$满足条件的k的最小值为1,所以对于当前的C和D,最大的A和B为$A=D-9t-1,B=D-7t-1$

那么如果A和B更小怎么办?

观察到在其他条件不变的情况下,只要$C$和$B$满足$Xc-Xb>6t$,那么这个魔法阵就一定成立,所以当$(a1<a2,b1<b2)$时,只要$a2$和$b2$能够和$C,D$组成魔法阵,$a1,b1$也一定能和$C,D$组成魔法阵,所以可以使用前缀和优化

然后又由乘法原理可得,当前魔法值作为$D$物品的个数为$SumD=SumA*SumB*SumC$

所以我们利用前缀和优化$SumA*SumB$

C的情况可以顺便在算D的时候算出来

那么还有一个问题是,我们枚举的D的范围是多少?

因为要统计前缀和,所以一定是要顺推下去的,由上面那张图我们可以知道,D的最大值为n,最小值则为当k=1且A=1的时候,所以D的最小值为$9*t+2$,再小是无法组成魔法阵的

同理可以枚举A

但是这个的情况又和枚举D的情况有一点不同

在其他条件不变的情况下,只要$C$和$B$满足$Xc-Xb>6t$,那么这个魔法阵就一定成立,所以当$(c1<c2,d1<d2)$时,只要$c1$和$d1$能够和$A,B$组成魔法阵,$c2,d2$也一定能和$A,B$组成魔法阵,所以可以使用后缀和优化

因为需要统计后缀和,所以需要逆推

枚举的范围:A的最大值为$(n-t*9-1)$(因为当k=1,D=n的时候A才最大),A的最小值则为1

所以就可以算出每个魔法值作为$A,B,C,D$物品的次数了,输出时直接输出当前魔法物品的魔法值的次数就可以了

#include <cstdio>
#include <algorithm>
using namespace std;
inline int read(){
int x=,f=;char c=getchar();
while(c<''||c>''){if(c=='-')f=-f;c=getchar();}
while(c>=''&&c<=''){x=(x<<)+(x<<)+c-'';c=getchar();}
return x*f;
}
#define N 50010
int n,m;
int a[N],b[N],c[N],d[N];
int x[N],vis[N];
int main(){
n=read();m=read();
for(int i=;i<=m;i++)
x[i]=read(),vis[x[i]]++;
for(int t=;t*<n;t++){
int sum=;
for(int D=*t+;D<=n;D++){
int A=D-*t-;
int B=A+*t;
int C=D-t;
sum+=vis[A]*vis[B];
c[C]+=vis[D]*sum;
d[D]+=vis[C]*sum;
}
sum=;
for(int A=n-*t-;A;A--){
int B=A+*t;
int C=B+*t+;
int D=A+*t+;
sum+=vis[C]*vis[D];
a[A]+=vis[B]*sum;
b[B]+=vis[A]*sum;
}
}
for(int i=;i<=m;i++){
printf("%d %d %d %d\n",a[x[i]],b[x[i]],c[x[i]],d[x[i]]);
}
return ;
}

参考资料:

学长的博客

(如果看不懂我的题解可以去看一下学长的qwq,我是看他的题解才懂这道题怎么写的,然后自己再归纳总结一下才写出来这篇博客)

[luogu2119]魔法阵 NOIP2016T4的更多相关文章

  1. 洛谷 P2119 魔法阵

    题目描述 六十年一次的魔法战争就要开始了,大魔法师准备从附近的魔法场中汲取魔法能量. 大魔法师有mm个魔法物品,编号分别为1,2,...,m1,2,...,m.每个物品具有一个魔法值,我们用X_iXi ...

  2. 【DFS】佳佳的魔法阵

    [vijos1284]佳佳的魔法阵 背景 也许是为了捕捉猎物(捕捉MM?),也许是因为其它原因,总之,佳佳准备设计一个魔法阵.而设计魔法阵涉及到的最关键问题,似乎就是那些带有魔力的宝石的摆放…… 描述 ...

  3. 洛谷P2119 魔法阵

    P2119 魔法阵 题目描述 六十年一次的魔法战争就要开始了,大魔法师准备从附近的魔法场中汲取魔法能量. 大魔法师有m个魔法物品,编号分别为1,2,...,m.每个物品具有一个魔法值,我们用Xi表示编 ...

  4. [NOIP2016普及组]魔法阵

    题目:洛谷P2119.Vijos P2012.codevs5624. 题目大意:有n件物品,每件物品有个魔法值.要求组成魔法阵(Xa,Xb,Xc,Xd),该魔法阵要满足Xa<Xb<Xc&l ...

  5. P2119 魔法阵

    原题链接  https://www.luogu.org/problemnew/show/P2119 YY同学今天上午给我们讲了这个题目,我觉得她的思路很好,特此写这篇博客整理一下. 50分:暴力枚举 ...

  6. 洛谷 题解 P2119【魔法阵】

    很好的一道数学推导题 45分做法 \(O(N^4)\) 暴力枚举四个材料 55分做法 从第一个约束条件可得到所有可行答案都是单调递增的,所以可以排序一遍,减少枚举量,可以拿到55分 100分做法 首先 ...

  7. 题解 【NOIP2016】魔法阵

    [NOIP2016]魔法阵 Description 六十年一次的魔法战争就要开始了,大魔法师准备从附近的魔法场中汲取魔法量. 大魔法师有m个魔法物品,编号分别为1,2,...,m.每个物品具有一个魔法 ...

  8. [NOIP2016PJ]魔法阵

    今天模拟赛的题,,,唯一没有Giao出来的题(不然我就AKIOI了~) 最开始没想到数学题,把所有部分分都说一遍吧: 35分:纯暴力O(M^4)枚举,对于每一组a,b,c,d验证其是否合法. 60分: ...

  9. 「Vijos 1284」「OIBH杯NOIP2006第二次模拟赛」佳佳的魔法阵

    佳佳的魔法阵 背景 也许是为了捕捉猎物(捕捉MM?),也许是因为其它原因,总之,佳佳准备设计一个魔法阵.而设计魔法阵涉及到的最关键问题,似乎就是那些带有魔力的宝石的摆放-- 描述 魔法阵是一个\(n ...

随机推荐

  1. django源码笔记-【1】(转)

    add by zhj:第二段代码有修改 原文:http://www.cnblogs.com/gaott/archive/2012/02/28/2371181.html 前言 Django是一个开放源代 ...

  2. "portrait"(竖屏) "landscape"(横屏) css设置

    取决于浏览器或者设备的方向,HTML元素总是会有"portrait"(竖屏) "landscape"(横屏) class. 你可以在css中如下使用它们: .p ...

  3. <span> 标签

    <span> 标签被用来组合文档中的行内元素. 如果不对 span 应用样式,那么 span 元素中的文本 与 其他文本不会有任何视觉上的差异.尽管如此,上例中的 span 元素仍然为 p ...

  4. C#调用VP 包含素材

    VS2012 +VP9.0 ***************** 自己运行的时只要修改VP里面素材的路径即可 链接: https://pan.baidu.com/s/1J6Bc5FcBYLZLgqe30 ...

  5. java反射之根据全类名创建对象

    现在的需求是根据类的全名.来创建对象 package 中介者设计模式; import java.util.Date; public class CreateObject { public static ...

  6. 【环境变量】删掉centos原有的openjdk并安装sun jdk

    一.卸载原有openjdk rpm -qa | grep java 之后,将展示出来的全部卸载掉,我这里是5个 rpm -e --nodeps java-1.7.0-openjdk-1.7.0.111 ...

  7. 摘自(http://www.ruanyifeng.com/blog/2011/07/linux_load_average_explained.html)

    理解Linux系统负荷   作者: 阮一峰 一.查看系统负荷 如果你的电脑很慢,你或许想查看一下,它的工作量是否太大了. 在Linux系统中,我们一般使用uptime命令查看(w命令和top命令也行) ...

  8. 对于session,request,cookie的理解

    session和request的生命周期 首先是session,比如我们在实现一个购物车功能时,在某一页面(这里称为页面A)选择了一些购物的商品,添加到购物车.那么当我们选择完成后点击我的购物车时会跳 ...

  9. Java: war包的作用及使用方法,如何解压后缀名为war的文件

    1.什么是war文件? 如果一个Web应用程序的目录和文件非常多,那么将这个Web应用程序部署到另一台机器上,就不是很方便了,我们可以将Web应用程序打包成Web归档(WAR)文件.这个过程和把Jav ...

  10. python 线程、多线程

    复习进程知识: python:主进程,至少有一个主线程 启动一个新的子进程:Process,pool 给每一个进程设定一下执行的任务:传一个函数+函数的参数 如果是进程池:map函数:传入一个任务函数 ...