前言
蒟蒻的第四场赛事,被E题完虐。
本题涉及题目:
A: [Arpa and a research in Mexican wave]
T1 Arpa and a research in Mexican wave
题面
原题地址: [A]
题意
e~~~,自己看吧后面的样例解释吧.
思路
我们发现t1~t5,^的个数分别为1,2,3,4,5;t14~t10,^的个数分别为1,2,3,4,5;又因为数据范围1 ≤ n ≤ 10^9, 1 ≤ k ≤ n, 1 ≤ t < n + k,所以必有代码中的分段求值……
代码
#include<bits/stdc++.h>
using namespace std;
int n,k,t;
int main()
{
scanf("%d%d%d",&n,&k,&t);
if(t<=k) printf("%d",t);
else
if(t>=n) printf("%d",n+k-t);
else
printf("%d",k);
return 0;
}
小结
水题不解释.
T2 Arpa and an exam about geometry
题面
原题地址: [B]
题意
有三个点a,b,c,问能否找到一个点k与一个角度p,使a,b以k旋转p度与b,c重叠.
思路
a,b的中垂线与b,c的中垂线的交点就是k,再满足线段a-b与b-c的长度相等,就行了.
但此时还有一个特例:a,b,c连成一条线时是没有k的!!!(WA的那次).
代码
#include<bits/stdc++.h>
using namespace std;
long long ax,ay,bx,by,cx,cy,x,y;
int main()
{
scanf("%lld%lld%lld%lld%lld%lld",&ax,&ay,&bx,&by,&cx,&cy);
x=(ax-bx)*(ax-bx)+(ay-by)*(ay-by);
y=(bx-cx)*(bx-cx)+(by-cy)*(by-cy);
if(x!=y) printf("No");//线段不等;
else
{
if((bx-cx)*(ay-by)==(ax-bx)*(by-cy)) printf("No");//一条线特判;
else printf("Yes");
}
}
小结
水题不解释.
T3 Five Dimensional Points
题面
原题地址: [C]
题意
好点:与其他任意两点形成的夹角都不小于90。
求n个点中的好点个数(五维空间中)。
思路
由于i点与任意两点的角度都要大于90度,在二维空间上,最多有4个点,三维最多6个,五维10个,再加上原点的那个点,max(n)=11.当n>11时,直接输出0。
(自己证为什么)。
公式优化
然后去看一下向量点乘,反三角函数,大于90度与小于90度三角函数值的不同。
此时我们就可以发现,我们只用判断向量点乘后的正负性就行了。
代码
#include<bits/stdc++.h>
using namespace std;
struct Point
{
int a,b,c,d,e;
}point[12];
int n,i,j,e,cut,ans[12];
bool p;
int compare(int k,int x,int y)//向量点乘;
{
return (point[k].a-point[x].a)*(point[k].a-point[y].a)+(point[k].b-point[x].b)*(point[k].b-point[y].b)+(point[k].c-point[x].c)*(point[k].c-point[y].c)+(point[k].d-point[x].d)*(point[k].d-point[y].d)+(point[k].e-point[x].e)*(point[k].e-point[y].e);
}
int main()
{
scanf("%d",&n);
if(n>11)
{
printf("0");
return 0;
}
for(i=1;i<=n;i++) scanf("%d%d%d%d%d",&point[i].a,&point[i].b,&point[i].c,&point[i].d,&point[i].e);
for(e=1;e<=n;e++)//枚所有点的所有形成角的可能;
{
p=true;
for(i=1;i<n;i++)
{
if(i==e) continue;//点不能重复;
for(j=i+1;j<=n;j++)
{
if(j==e) continue;//点不能重复;
if(compare(e,i,j)>0)//判 断正负性;
{
p=false;
break;
}
}
if(!p) break;
}
if(p)
{
ans[0]++;
cut++;
ans[cut]=e;
}
}
for(i=0;i<=cut;i++) printf("%d\n",ans[i]);
return 0;
}
小结
向量点乘和反三角函数都没学过…
T4 Polycarpus’Dice
题面
原题地址: [D]
题意
不会讲QWQ…
思路
推公式…
推法:
每一个点所能拥有的正常偏差值max(0,d[i]-A-1+n)(负值是不允许的)和特殊偏差值max(0,A-sum+d[i]-1))的和。
(完全无法组织语言~~~,自己看吧~)
代码
#include<bits/stdc++.h>
#define LL long long
using namespace std;
LL d[200000],n,A,i,sum;
int main()
{
for(scanf("%lld%lld",&n,&A),i=0;i<n;i++) scanf("%lld",&d[i]),sum+=d[i];
for(i=0;i<n;i++) printf("%lld ",max((LL)0,d[i]-A-1+n)+max((LL)0,A-sum+d[i]-1));
return 0;
}
小结
推推公式就过了。
T5 Little Elephant and LCM
题面
原题地址: [E]
题意
有n个数,求序列b使bi不大于ai,且b的最小公倍数等于b中的最大值。
思路
由于最小公倍数等于最大值,而max(bi)=ai,所以只用找出a中的最大值,然后枚举b中的最大值就行了。
此时b中的所有数都是maxn的因数,此时要知道maxn的因数i有多少种选择,那比i小的因数同样也有这些选择(但这样会有重复)。
先将a排序,再在a中找到因数bi的位置,求出比bi大的数有k个,ans+=i^k,但这时,后面的数如b(i+1),就会将之前的数重新再算一遍,所以k要改成比bi大且b(i+1)小的数的个数。
最后一个数要特判,因为最后一个数必修要出现(前面的数可以不出现)。
代码
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 2e5 + 10, mod = 1e9 + 7;
int sum[N],n,mx,a[N],ans=0,s[N],top=0,res;
void pp(int n)
{
for(int i=1;i<=sqrt(n+0.5);i++)//扫因数;
{
if(n%i==0)
{
s[++top]=i;
if(i*i!=n) s[++top]=n/i;
}
}
}
int pow_mod(int a,int k,int mod)//快速幂;
{
int ans=1;
for(;k;k>>=1,a=(a*a)%mod) if(k&1) ans=(ans*a)%mod;
return ans;
}
void solve(int m)
{
sort(s+1,s+top+1),res=1;
for(int i=1;i<top;i++) res=(res*pow_mod(i,sum[s[i+1]-1]-sum[s[i]-1],mod))%mod;
ans=(ans+(res*(pow_mod(top,sum[N-100]-sum[s[top]-1],mod)-pow_mod(top-1,sum[N-100]-sum[s[top]-1],mod)))%mod)%mod;
}
main()
{
scanf("%lld",&n);
for(int i=1;i<=n;i++) scanf("%lld",&a[i]),mx=max(mx,a[i]),sum[a[i]]++;
for(int i=1;i<=N-10;i++) sum[i]+=sum[i-1];
for(int i=1;i<=mx;i++)
{
top=0,pp(i),solve(i);
}
printf("%lld\n",(ans+mod)%mod);
}
小结
数学优化我是怎么都没推出来(暴力不断WA)。
总结
怎么又是一堆没学过的东西aaaaa!!!!!