入门赛4

前言

蒟蒻的第四场赛事，被E题完虐。

本题涉及题目:
A: [Arpa and a research in Mexican wave]

B: [Arpa and an exam about geometry]

C: [Five Dimensional Points]

D: [PolycarpusDice]

E: [Little Elephant and LCM]

---

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！！！！！