codeforces #198 div.2 解题报告
A 题 签到题。不过我挺水的中途被hack成功了。。。后来一直在测数据。
这题用gcd求出lcm。然后根据公式就可以直接计算获得
#include<stdio.h>
int Gcd(int a,int b)
{
if(b>0)
{
a=a%b;
return Gcd(b,a);
}
else return a;
}
int main()
{
int x,y,a,b,i,j,c,k;
while(scanf("%d %d %d %d",&x,&y,&a,&b)!=EOF)
{
c=0;
if(x>y) k=Gcd(x,y);
else k=Gcd(y,x);
k=x*y/k;
if(a%k) printf("%d\n",b/k-a/k);
else printf("%d\n",b/k-a/k+1);
}
}
B题，这道题目出题者明显出难了。作者也在事后给了道歉原话：I want to apologize for not estimating the real difficulty of this task. It turns out that it was more complicated than we thought it might be.
我就是我这么机智的弱菜怎么不会。其实这个题目用到的知识我是见过的。但没想到。这里首先你必须知道任意多边形的面积公式。
|x1 y1| |x2 y2| |x3 y3| |x5 y5| |x6 y6| |x7 y7| |xn-1 yn-1| |xn yn | |x2 y2|+|x3 y3|+|x4 y4|+|x6 y6|+|x7 y7|+|x8 y8|...+|xn yn |+ |x1 y1 |*0.5=S
这里用到的知识点就是你要枚举平面内所有的直线，再根据直线枚举所有的点与这条直线构成的三角形的面积。
网上的一般代码都是这个公式的简化，所以希望读者能过牢记这个万能的计算任意多边形（凹、凸）的公式。
而这个叉积的面积当点为逆时针排列时S为负数，当点为正时针排列时S为正数。也就是说如果从与这条直线所构成的三角形中选两个组成四边形，那么这两个点必须分布在直线两侧，使得他们分别与直线构成的三角形通过公式计算出的值为一正一负。设直线为 ab那么枚举的点c与之组成的序列acb可能是正也可能是负数。我们要用贪婪的方式保留这两个最大子面积。
#include<stdio.h>
#include<string.h>
#include<math.h>
typedef struct Point
{
int x,y;
}Point;
Point p[305];
double area(Point a,Point b,Point c)
{
return a.x*b.y-b.x*a.y+b.x*c.y-c.x*b.y+c.x*a.y-a.x*c.y; /*公式直接计算*/
}
double max(double s1,double s2)
{
return s1>s2?s1:s2;
}
int main()
{
int n;
while(scanf("%d",&n)!=EOF)
{
int i,j,k;
double s1,s2,s3,s;
for(i=0;i<n;i++) scanf("%d %d",&p[i].x,&p[i].y);
s3=0;
for(i=0;i<n;i++)
{
for(j=i+1;j<n;j++)
{
s1=s2=0;
for(k=0;k<n;k++)
{
if(k!=i&&k!=j)
{
s=area(p[i],p[k],p[j]);
if(s>0) s1=max(s1,fabs(s)/2);
else if(s<0) s2=max(s2,fabs(s)/2);
}
}
if(s1!=0&&s2!=0) s3=max(s1+s2,s3); //子面积必须分布在直线两侧，若某个子面积为0.则说明所有点在这条直线的同一侧。见测试#23
}
}
printf("%.6lf\n",s3);
}
return 0;
}
C 题
Consider 6 possible routes:
- [2, 3, 5]: total distance traveled: |2 – 0| + |3 – 2| + |5 – 3| = 5;
- [2, 5, 3]: |2 – 0| + |5 – 2| + |3 – 5| = 7;
- [3, 2, 5]: |3 – 0| + |2 – 3| + |5 – 2| = 7;
- [3, 5, 2]: |3 – 0| + |5 – 3| + |2 – 5| = 8;
- [5, 2, 3]: |5 – 0| + |2 – 5| + |3 – 2| = 9;
- [5, 3, 2]: |5 – 0| + |3 – 5| + |2 – 3| = 8.
The average travel distance is 
=
=
.
这周似乎都是数学题，这道题是组合问题。首先你要推一下公式。
我们看到最终结果中分母是6。这个6就是A(3,3);
而分子是不同元素相减的结果。可以发现除了0以外其他元素都是数组中的。
可以将分子C分解成c1=|2-0|+|2-0|+|3-0|+|3-0|+.....+|5-0|+|5-0|;
c2=|3-2|+|5-3|+.....|2-3|;
可以发现每个组合会出现2次即A(2,2)次；
那么推广一下c1=sum(a[n])*A(n-1,n-1);
c2=C(2,n)*A(n-1,n-1)*sum(abs(a[i]-a[j]));
c2+c1=A(n-1,n-1)*(sum(a[n])+2*sum(a[i]-a[j])); a[i]>a[j]
与分母A(n,n)进行化简得；
结果为(sum(a[n])+2*sum(a[i]-a[j]))/n; a[i]>a[j]
可是要枚举所有的组合是会超时的，但根据先用冒泡排序，再根据c2的组成结构我们也能很好得发现规律。
如例子中c2=2*((3-1)+(3-2)+(2-1));
轻松得发现3出现了2次，2出现1次，-2出现1次，-1出现2次，1出现0次；
这说明当排序好后c2=c2+(i-j)*a[i]; 即a[i]在c2中的贡献是比a[i]小的数字个数减去比a[i]数字大的个数乘a[i];
例如例子中比2大的数只有1个，比2小的数也只有1个那么2的贡献为0；例如数字1；比1大的有两个那么1的贡献为-2*1;
#include<stdio.h>
#include<stdlib.h>
#include <algorithm>
using namespace std;
__int64 gcd(__int64 a,__int64 b)
{
if(b==0) return a;
return gcd(b,a%b);
}
__int64 a[100005];
int main()
{
int n;
while(scanf("%d",&n)!=EOF)
{
int i,j;
__int64 c1,c2,cnt,g;
c1=c2=0;
for(i=0;i<n;i++)
{
scanf("%I64d",&a[i]);
c2+=a[i];
}
/*for(i=0;i<n;i++)
{
for(j=i+1;j<n;j++)
{
if(a[i]>a[j]) c1+=a[i]-a[j];
else c1+=a[j]-a[i];
}
}*/
sort(a,a+n); //
for(i=0,j=n-1;i<n;i++,j--) //
{ //
c1+=(i-j)*a[i]; //
} //
c1*=2;
cnt=c1+c2;
g=gcd(cnt,n);
printf("%I64d %I64d\n",cnt/g,n/g);
}
}
