洛谷p2119魔法阵

一蓑烟雨任生平 一蓑烟雨任生平     2022-09-29     177

关键词:

题目描述

六十年一次的魔法战争就要开始了,大魔法师准备从附近的魔法场中汲取魔法能量。

大魔法师有m个魔法物品,编号分别为1,2,...,m。每个物品具有一个魔法值,我们用Xi表示编号为i的物品的魔法值。每个魔法值Xi是不超过n的正整数,可能有多个物品的魔法值相同。

大魔法师认为,当且仅当四个编号为a,b,c,d的魔法物品满足xa<xb<xc<xd,Xb-Xa=2(Xd-Xc),并且xb-xa<(xc-xb)/3时,这四个魔法物品形成了一个魔法阵,他称这四个魔法物品分别为这个魔法阵的A物品,B物品,C物品,D物品。

现在,大魔法师想要知道,对于每个魔法物品,作为某个魔法阵的A物品出现的次数,作为B物品的次数,作为C物品的次数,和作为D物品的次数。

输入输出格式

输入格式:

 

输入文件的第一行包含两个空格隔开的正整数n和m。

接下来m行,每行一个正整数,第i+1行的正整数表示Xi,即编号为i的物品的魔法值。

保证1 le n le 150001n15000,1 le m le 400001m40000,1 le Xi le n1Xin。每个Xi是分别在合法范围内等概率随机生成的。

 

输出格式:

 

共输出m行,每行四个整数。第i行的四个整数依次表示编号为i的物品作 为A,B,C,D物品分别出现的次数。

保证标准输出中的每个数都不会超过10^9。

每行相邻的两个数之间用恰好一个空格隔开。

 

输入输出样例

输入样例#1: 复制
30 8
1
24
7
28
5
29
26
24
输出样例#1: 复制
4 0 0 0
0 0 1 0
0 2 0 0
0 0 1 1
1 3 0 0
0 0 0 2
0 0 2 2
0 0 1 0
输入样例#2: 复制
15 15
1 
2 
3 
4 
5
6 
7 
8 
9
10
11
12
13
14
15
输出样例#2: 复制
5 0 0 0
4 0 0 0
3 5 0 0
2 4 0 0
1 3 0 0
0 2 0 0
0 1 0 0
0 0 0 0
0 0 0 0
0 0 1 0
0 0 2 1
0 0 3 2
0 0 4 3
0 0 5 4
0 0 0 5

说明

【样例解释1】

共有5个魔法阵,分别为:

物品1,3,7,6,其魔法值分别为1,7,26,29;

物品1,5,2,7,其魔法值分别为1,5,24,26;

物品1,5,7,4,其魔法值分别为1,5,26,28;

物品1,5,8,7,其魔法值分别为1,5,24,26;

物品5,3,4,6,其魔法值分别为5,7,28,29。

以物品5为例,它作为A物品出现了1次,作为B物品出现了3次,没有作为C物品或者D物品出现,所以这一行输出的四个数依次为1,3,0,0。

此外,如果我们将输出看作一个m行4列的矩阵,那么每一列上的m个数之和都应等于魔法阵的总数。所以,如果你的输出不满足这个性质,那么这个输出一定不正确。你可以通过这个性质在一定程度上检查你的输出的正确性。

【数据规模】

技术分享

技术分享
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 40010
using namespace std;
int n,m;
int x[MAXN],a[MAXN],b[MAXN],c[MAXN],d[MAXN];
bool judge(int i,int j,int k,int l){
    if(x[i]<x[j]&&x[j]<x[k]&&x[k]<x[l]&&x[j]-x[i]==2*(x[l]-x[k])&&double(x[j]-x[i])<double(x[k]-x[j])/3.0){
        a[i]++;b[j]++;
        c[k]++;d[l]++;
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)    scanf("%d",&x[i]);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=m;j++)
            for(int k=1;k<=m;k++)
                for(int l=1;l<=m;l++)
                    judge(i,j,k,l); 
    for(int i=1;i<=m;i++)
        cout<<a[i]<<" "<<b[i]<<" "<<c[i]<<" "<<d[i]<<endl;
} 
m^4的暴力,只有35分
技术分享
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 40010
using namespace std;
int n,m,tot;
int num[MAXN],sum[MAXN];
int x[MAXN],a[MAXN],b[MAXN],c[MAXN],d[MAXN];
bool judge(int i,int j,int k,int l){
    if(x[i]<x[j]&&x[j]<x[k]&&x[k]<x[l]&&x[j]-x[i]==2*(x[l]-x[k])&&double(x[j]-x[i])<double(x[k]-x[j])/3.0){
        a[i]++;b[j]++;
        c[k]++;d[l]++;
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        scanf("%d",&x[i]);
        sum[x[i]]++;
    }    
    for(int i=1;i<=15000;i++)
        if(sum[i])    num[++tot]=i;
    for(int i=1;i<=tot;i++)
        for(int j=i+1;j<=tot;j++)
            for(int k=j+1;k<=tot;k++){
                if(4*num[j]-3*num[i]-num[k]>=0)    continue;
                for(int l=k+1;l<=tot;l++)
                    if(num[j]-num[i]==2*(num[l]-num[k])){
                        a[num[i]]+=sum[num[j]]*sum[num[k]]*sum[num[l]];
                        b[num[j]]+=sum[num[i]]*sum[num[k]]*sum[num[l]];
                        c[num[k]]+=sum[num[i]]*sum[num[j]]*sum[num[l]];
                        d[num[l]]+=sum[num[i]]*sum[num[j]]*sum[num[k]];
                    }
            }
    for(int i=1;i<=m;i++)
        cout<<a[x[i]]<<" "<<b[x[i]]<<" "<<c[x[i]]<<" "<<d[x[i]]<<endl;
} 
/*
30 8
1
24
7
28
5
29
26
24
*/
70分的暴力
技术分享
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int a[40050];
int ans[15050][10];
int t[15050];
int main() {
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        cin>>a[i];
        t[a[i]]++; 
    }
    for(int i=1;i<=(n-2)/9;i++){
        int sum=0,x=9*i+1;
        for(int xd=9*i+2;xd<=n;xd++){
            int xc=xd-i;
            sum+=t[xd-x]*t[xd-x+2*i];
            ans[xd][4]+=sum*t[xc];
            ans[xc][3]+=sum*t[xd]; 
        }
        sum=0;
        x=n-9*i-1;
        for(int xa=n-9*i-1;xa>=1;xa--){
            int xb=xa+2*i;
            sum+=t[n-(x-xa)]*t[n-(x-xa)-i];
            ans[xa][1]+=sum*t[xb];
            ans[xb][2]+=sum*t[xa];
        }
    }
    for(int i=1;i<=m;i++){
        for(int j=1;j<=4;j++){
            cout<<ans[a[i]][j]<<" ";
        }
        cout<<endl;
    }
    return 0;
}
数学推论 AC

 

[p2119]魔法阵(模拟?搜索?)(代码片段)

很玄学我暴力都没做出来#include<cstdio>intn,m,x[40001],vis[15001],a[15001],b[15001],c[15001],d[15001];intmain()//freopen("magic.in","r",stdin);//freopen("magic.out","w",stdout);scanf("%d%d",&n,&m); 查看详情

[noip2016普及组]魔法阵

题目:洛谷P2119、VijosP2012、codevs5624。题目大意:有n件物品,每件物品有个魔法值。要求组成魔法阵(Xa,Xb,Xc,Xd),该魔法阵要满足Xa<Xb<Xc<Xd,Xb-Xa=2(Xd-Xc),并且Xb-Xa<(Xc-Xb)/3。求每件物品作为a、b、c、d的次数。解题思路:这... 查看详情

2016魔法阵

题解:枚举+超级无敌大优化。看了半天,以为是动归,仔细一想,感觉不太行,只能枚举。所以就想到了枚举+优化,想了半天,不知道怎么优化(我好渣啊),代码不提供,算法:枚举+超级无敌大优化。 查看详情

魔法阵(noip2016普及组第四题)

描述六十年一次的魔法战争就要开始了,大魔法师准备从附近的魔法场中汲取魔法能量。大魔法师有m个魔法物品,编号分别为1,2,...,m。每个物品具有一个魔法值,我们用x_ixi?表示编号为i的物品的魔法值。每个魔法值x_ixi?是不超... 查看详情

vijos1285&1283&1282&1284佳佳的魔法药水/魔杖/魔法照片/魔法阵

题目链接:https://vijos.org/p/1285https://vijos.org/p/1283https://vijos.org/p/1282https://vijos.org/p/1284  查看详情

noip2016普及组复赛魔法阵(代码片段)

...分别枚举xa和xd,根据之间的关系,用前缀和求出每一种魔法阵的每一种物品的次数。#include<cmath>#include<iostream>#include<cstdio>#include<cstdlib>#in 查看详情

c.奇奇怪怪的魔法阵(未搞懂)(代码片段)

C.奇奇怪怪的魔法阵题意:n个点m条边,定义集合S为独立集,当且仅当任意x,y∈S,x与y之间没有边。空集也是独立集现在对于每一个点的集合T,有多少子集为独立集设N=0,1,…,n-1,AT=∑S⊂T[S是独立集]A_T=... 查看详情

洛谷——p2387魔法森林

题目描述为了得到书法大家的真传,小E同学下定决心去拜访住在魔法森林中的隐士。魔法森林可以被看成一个包含n个节点m条边的无向图,节点标号为1,2,3,…,n,边标号为1,2,3,…,m。初始时小E同学在1号节点,隐士则住在n... 查看详情

bzoi——3343:教主的魔法||洛谷——p2801教主的魔法

http://www.lydsy.com/JudgeOnline/problem.php?id=3343  || https://www.luogu.org/problem/show?pid=2801题目描述教主最近学会了一种神奇的魔法,能够使人长高。于是他准备演示给XMYZ信息组每个英雄看。于是N个英雄们又一次聚集在了一起,这... 查看详情

牛客挑战赛53c.奇奇怪怪的魔法阵(状压)(代码片段)

Link我们要求fT=∑S∈T[S为独立集]f_T=\\sum\\limits_S\\inT[S为独立集]fT​=S∈T∑​[S为独立集]考虑集合TTT的最低位二进制位代表的集合是yyy(这个点设为uuu),其他二进制位代表的是xxx那么T=x+yT=x+yT=x+y显然,若SSS集... 查看详情

牛客挑战赛53c.奇奇怪怪的魔法阵(状压)(代码片段)

Link我们要求fT=∑S∈T[S为独立集]f_T=\\sum\\limits_S\\inT[S为独立集]fT​=S∈T∑​[S为独立集]考虑集合TTT的最低位二进制位代表的集合是yyy(这个点设为uuu),其他二进制位代表的是xxx那么T=x+yT=x+yT=x+y显然,若SSS集... 查看详情

洛谷——p1875佳佳的魔法药水

 P1875佳佳的魔法药水题目描述发完了k张照片,佳佳却得到了一个坏消息:他的MM得病了!佳佳和大家一样焦急万分!治好MM的病只有一种办法,那就是传说中的0号药水……怎么样才能得到0号药水呢?你要知道佳佳的... 查看详情

洛谷——p2483[sdoi2010]魔法猪学院

https://www.luogu.org/problem/show?pid=2483题目描述iPig在假期来到了传说中的魔法猪学院,开始为期两个月的魔法猪训练。经过了一周理论知识和一周基本魔法的学习之后,iPig对猪世界的世界本原有了很多的了解:众所周知,世界是由元... 查看详情

洛谷p2387魔法森林

这个题就是先按照a从小到大排序,然后lct尽可能维护b值最小,在这个过程中寻找答案#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<cmath>usingnamespacestd;intx,y,a,b;structin{intf,l 查看详情

洛谷3833shoi2012魔法树

【题解】  树链剖分模板题。。#include<cstdio>#include<algorithm>#include<queue>#defineN500010#definergregister#definels(u<<1)#definers(u<<1|1)#definemid((a[u].l+a[u].r)>>1)#defin 查看详情

洛谷p2801教主的魔法

题目描述教主最近学会了一种神奇的魔法,能够使人长高。于是他准备演示给XMYZ信息组每个英雄看。于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1、2、……、N。每个人的身高一开始都是不超过100... 查看详情

洛谷p2483[sdoi2010]魔法猪学院

P2483[SDOI2010]魔法猪学院 k短路模板1#include<queue>2#include<cstdio>3#include<cstring>4#include<algorithm>5usingnamespacestd;6constintmaxn=10005;7intn,m,s,t,k,num,head1[maxn*100],h 查看详情

洛谷——p1875佳佳的魔法药水

https://www.luogu.org/problemnew/show/1875题目背景发完了k张照片,佳佳却得到了一个坏消息:他的MM得病了!佳佳和大家一样焦急万分!治好MM的病只有一种办法,那就是传说中的0号药水……怎么样才能得到0号药水呢?你要知道... 查看详情