奋斗群群赛11总结与心得

总体情况

https://vjudge.net/contest/186188
本次题目非常难，我只能A第一道。下午晚上又搞了两个小时才搞定B。

T1

题目

给出a,b,n,k,输出在n，k之间a,b最小公倍数的个数。

思路

一开始我脑残了一下WA了好几次，后来又从n-k暴枚，我都不知道我干了什么。代码如下。
使用了？以及：之后好用多了。B题也是，不然太长了。

#include<bits/stdc++.h>
using namespace std;
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
int lcm(int a,int b){return a*b/gcd(a,b);}
int main()
{
int x,y,a,b,s=0,i,j;
cin>>x>>y>>a>>b;
int t=lcm(x,y);
i=a/t*t;j=b/t*t;
if (a==i) s++;
s+=(j-i)/t;
cout<<s;
}

T2

题目

给出平面直角坐标系中n（n<=300）个点，求能组成的面积最大的四边形的面积。

思路

一开始不会，后来大佬告诉我，枚举i,j两个点，再枚举k，判断k在ij直线的上方还是下方，把上下面积最大的两个三角形的面积和加起来。交上去之后T了。我把函数加个inline，在23个点wa了。后来我去调数据，这个测试点只有4个点，组成的是一个凹四边形，这个代码扫出了最外面大的三角形的面积。最后我不得不对n=4的情况特判，以998ms危险AC.(要求1秒之内)我果然是个蒟蒻。

#include<bits/stdc++.h>
using namespace std;
double x[310],y[310],d1,d2,d3,point[310][310],mianji=0;

inline double helen(double a,double b,double c){double p=(a+b+c)/2;return sqrt(p*(p-a)*(p-b)*(p-c));}//海伦公式
inline double distance(double x1,double y1,double x2,double y2){return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));}//两点间距离公式
inline bool rk(int i,int j){return x[i]==x[j]?0:1;}//两点的横坐标是否相等（两点连线是否和y轴平行）

inline bool upordown(int i,int j,int p)//p点在i,j连线上方还是下方
{
double k,b;
if (y[i]==y[j]) {k=0;b=y[i];}
else 
  {
  k=(y[j]-y[i])/(x[j]-x[i]);
  b=(y[j]*x[i]-y[i]*x[j])/(x[i]-x[j]);
  }
double t=x[p]*k+b;
return t<y[p]?1:0;
}

inline bool in_san(int i,int j,int k,int p)//p点是否在i,j,k三点组成的三角形之内，从p点以竖直方向作一条射线，判断它与三角形是否刚好有一个交点（交点是顶点不算）
{
int jiaodian=0;double k1,b1,k2,b2,k3,b3;
if (y[i]==y[j]) {k1=0;b1=y[i];}
else 
  {
  k1=(y[j]-y[i])/(x[j]-x[i]);
  b1=(y[j]*x[i]-y[i]*x[j])/(x[i]-x[j]); 
  } 
if (y[i]==y[k]) {k2=0;b2=y[i];}
else
  {
  k2=(y[k]-y[i])/(x[k]-x[i]);
  b2=(y[k]*x[i]-y[i]*x[k])/(x[i]-x[k]);
  }
if (y[j]==y[k]) {k3=0;b3=y[j];}
else 
  {
  k3=(y[k]-y[j])/(x[k]-x[j]);
  b3=(y[k]*x[j]-y[j]*x[k])/(x[j]-x[k]);
  }
if (x[p]*k1+b1>min(y[i],y[j])&&x[p]*k1+b1<max(y[i],y[j])&&x[p]!=x[i]&&x[p]!=x[j]) jiaodian++;
if (x[p]*k2+b2>min(y[i],y[k])&&x[p]*k2+b2<max(y[i],y[k])&&x[p]!=x[i]&&x[p]!=x[k]) jiaodian++; 
if (x[p]*k3+b3>min(y[j],y[k])&&x[p]*k2+b2<max(y[j],y[k])&&x[p]!=x[i]&&x[p]!=x[k]) jiaodian++; 
return jiaodian==1?1:0;
}

inline void tepan(int n)//n=4时n^4暴力特判求四边形面积
{
int i,j,k,l;double d4,d5;
for (i=1;i<=n;i++) for (j=1;j<=n;j++) if (i!=j&&rk(i,j))
  {
  d1=point[i][j];
  for (k=1;k<=n;k++) for (l=1;l<=n;l++) if (k!=i&&k!=j&&l!=i&&l!=j&&k!=l&&!in_san(i,j,k,l)&&upordown(i,j,k)!=upordown(i,j,l))//都不一样，l不在三角形之内，K和l对于i，j直线的方向不一样
    {
    d2=point[i][k];
    d3=point[j][k];
    d4=point[i][l];
    d5=point[j][l];
    double t=helen(d1,d2,d3)+helen(d1,d4,d5);
    mianji=max(mianji,t);
    }
  }
}

int main()
{
int n,i,j,k;double up=0,down=0;
cin>>n;
for (i=1;i<=n;i++) scanf("%lf%lf",&x[i],&y[i]);
for (i=1;i<=n-1;i++) for (j=i+1;j<=n;j++) point[i][j]=point[j][i]=distance(x[i],y[i],x[j],y[j]);
if (n==4) {tepan(n);printf("%lf",mianji);return 0;}
for (i=1;i<=n;i++) for (j=1;j<=n;j++) if (i!=j&&rk(i,j))
  {
  up=down=0;
  d1=point[i][j];
  for (k=1;k<=n;k++) if (k!=i&&k!=j)
    {
    d2=point[i][k];
    d3=point[j][k];
    double t=helen(d1,d2,d3);
    if (upordown(i,j,k)) up=max(up,t);//上下求最大
    else down=max(down,t);
    }
  mianji=max(mianji,up+down); 
  }
printf("%lf",mianji);
}

T3

题目

给出一数字串，对于该串的每一个全排列，计算a1+|a2-a1|+|a3-a2|+…+|a[n]-a[n-1]|之和的平均值，输出两个互质整数作为分数的分子和分母。

思路

在打的时候推错了公式一直都做不出。首先把a1单独拿出来讨论，在所有情况中a数组中的每一个元素都出现了（n-1）!次。然后固定第2个元素，在这（n-1）!种可能性中该元素放在第2个的可能性有(n-2)!次。交换这两个元素，结果乘以2；接下来不断推公式就可以了。

#include<bits/stdc++.h>
using namespace std;
const int boss=1e5;
int a[boss+10];
long long gcd(long long a,long long b)
{
return b==0?a:gcd(b,a%b);
}
void yuefen(long long &s,long long &n)
{
long long p=gcd(s,n);
s/=p;n/=p;
}
int main()
{
long long s=0,fenzi,fenmu,sum,n,i;
scanf("%I64d",&n);
for (i=1;i<=n;i++)
  {
  scanf("%d",&a[i]);
  s+=(long long)a[i];
  }//求和
sort(a+1,a+n+1);
sum=s;
for (i=1;i<n;i++)
  {
  sum-=a[i];
  s+=2*(sum-(n-i)*a[i]);
  }
yuefen(s,n);
printf("%I64d %I64d",s,n);
}

T4

题目

给出一个1-n数列,对其进行冒泡排序,每排序一次在图里加一条边,求排序完毕后最大不直接相邻点的数目.(我应该没有翻译错吧.)

思路

这题可以转化成最长不下降子序列,就是必须用nlogn算法来做.

#include<bits/stdc++.h>
using namespace std;
const int boss=1e5;
int a[boss+10],d[boss+10];
int main()
{
int n,i,j,len=1;
cin>>n;
for (i=1;i<=n;i++) scanf("%d",&a[i]);
if (n==0) return printf("0"),0;
d[1]=a[1];
for (i=2;i<=n;i++)
  {
  if (a[i]>=d[len]) d[++len]=a[i];
  else
    {
    int t=upper_bound(d+1,d+len+1,a[i])-d;
    d[t]=a[i];
    }
  }
cout<<len;
}

T5

题目

给出一个1-n数列,其中有个别数字一定在其位置上,其他的位置是-1,对剩下的数字进行错排,求方法的总数mod1000000007.

思路

可以利用错排公式.

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int boss=2000,mod=1e9+7;
ll a[boss+10],b[boss+10],dp[boss+10];
int main()
{
ll x=0,y=0,n,i;
cin>>n;
for (i=1;i<=n;i++) 
  {
  scanf("%I64d",&a[i]);
  if (a[i]!=-1) b[a[i]]=1;
  }
for (i=1;i<=n;i++) if (a[i]==-1) if (b[i]) x++;else y++;
dp[0]=1;
for (i=1;i<=x;i++) dp[0]=dp[0]*i%mod;
for (i=1;i<=y;i++)
  {
  dp[i]=(x+i-1)*dp[i-1]%mod;
  if (i>1) dp[i]=(dp[i]+(i-1)*dp[i-2]%mod)%mod;//错排公式f(i)=(i-1)*(f(i-1)+f(i-2)),f(2)=1,f(3)=2.
  }
printf("%I64d",dp[y]);
}