ACM 第十六天

计算几何

练习题:

F - Beauty Contest

POJ - 2187

Bessie, Farmer John's prize cow, has just won first place in a bovine beauty contest, earning the title 'Miss Cow World'. As a result, Bessie will make a tour of N (2 <= N <= 50,000) farms around the world in order to spread goodwill between farmers and their cows. For simplicity, the world will be represented as a two-dimensional plane, where each farm is located at a pair of integer coordinates (x,y), each having a value in the range -10,000 ... 10,000. No two farms share the same pair of coordinates.

Even though Bessie travels directly in a straight line
between pairs of farms, the distance between some farms can be quite
large, so she wants to bring a suitcase full of hay with her so she has
enough food to eat on each leg of her journey. Since Bessie refills her
suitcase at every farm she visits, she wants to determine the maximum
possible distance she might need to travel so she knows the size of
suitcase she must bring.Help Bessie by computing the maximum distance
among all pairs of farms.

Input

* Line 1: A single integer, N

* Lines 2..N+1: Two space-separated integers x and y specifying coordinate of each farm

Output

* Line 1: A single integer that is the squared distance between
the pair of farms that are farthest apart from each other.

Sample Input

4
0 0
0 1
1 1
1 0

Sample Output

2

Hint

Farm 1 (0, 0) and farm 3 (1, 1) have the longest distance (square root of 2)

 #include<stdio.h>
 #include<iostream>
 #include<algorithm>
 #include<string.h>

 using namespace std;

 struct Point
 {
     int x,y;
     Point(int _x=,int _y=)
     {
         x=_x;
         y=_y;
     }
     Point operator -(const Point &b) const
     {
          return Point(x-b.x,y-b.y);
     }
     int operator ^(const Point &b) const
     {
         return x*b.y-y*b.x;
     }
     int operator *(const Point &b) const
     {
         return x*b.x+y*b.y;
     }
     void input()
     {
         scanf("%d%d",&x,&y);
     }
 };

 int dist2(Point a,Point b)
 {
     return (a-b)*(a-b);
 }
 const int maxn=;
 Point list[maxn];
 int Stack[maxn],top;
 bool _cmp(Point p1,Point p2)
 {
     int tmp=(p1-list[])^(p2-list[]);
     if(tmp>) return true;
     else if(tmp== && dist2(p1,list[])<= dist2(p2,list[]))
         return true;
     else return false;
 }
 void Graham(int n)
 {
     Point p0;
     int k=;
     p0=list[];
     for(int i=;i<n;i++)
     if(p0.y>list[i].y || (p0.y==list[i].y && p0.x>list[i].x))
     {
         p0=list[i];
         k=i;
     }
     swap(list[],list[k]);
     sort(list+,list+n,_cmp);
     if(n==)
     {
         top=;
         Stack[]=;
         return ;
     }
     if(n==)
     {
         top=;
         Stack[]=;
         Stack[]=;
         return ;
     }
     Stack[]=;
     Stack[]=;
     top=;
     for(int i=;i<n;i++)
     {
         while(top> && ((list[Stack[top-]]-list[Stack[top-]])^(list[i]-list[Stack[top-]]))<=)
             top--;
         Stack[top++]=i;
     }
 }
 //旋转卡壳，求两点间距离平方的最大值
 int rotating_calipers(Point p[],int n)
 {
     int ans=;
     Point v;
     int cur=;
     for(int i=;i<n;i++)
     {
         v=p[i]-p[(i+)%n];
         while((v^(p[(cur+)%n]-p[cur]))<)
             cur=(cur+)%n;
         ans=max(ans,max(dist2(p[i],p[cur]),dist2(p[(i+)%n],p[(cur+)%n])));
     }
     return ans;
 }
 Point p[maxn];
 int main()
 {
     int n;
     while(scanf("%d",&n)==)
     {
         for(int i=;i<n;i++)
         list[i].input();
         Graham(n);
         for(int i=;i<top;i++)
          p[i]=list[Stack[i]];
          printf("%d\n",rotating_calipers(p,top));
     }
     return ;
 }

A - Arcade Game

Gym - 100814A

Arcade mall is a new modern mall. It has a new hammer game called "Arcade Game". In this game you're presented with a number n which is hanged on a wall on top of a long vertical tube, at the bottom of the tube there is a button that you should hit with your hammer.

When you hit the button with all your force (as you always do), a ball is pushed all over the tube and hit the number n. The number n flies in the air and it's digits fall back in any random permutation with uniform probability.

If the new number formed is less than or equal to the previous number, the game ends and you lose what ever the new number is. Otherwise (if the number is greater than the previous number), you are still in the game and you should hit the button again.

You win if the new number formed is greater than the previous number and it is equal to the greatest possible permutation number.

Can you compute the probability of winning?

Input

The first line of the input contains the number of test cases T. Following that there are T lines represents T test cases. In each line, there is a single integer (1 ≤ n ≤ 10⁹) the target number. The digits of n are all unique, which means that any 2 digits of n are different.

Output

For each test case, print one line containing the answer. Print the answer rounded to exactly 9 decimal digits.

Examples

Input

3
952
925
592

Output

0.000000000
0.166666667
0.194444444

Note

In the first test case, the answer is 0 because 952 is greater than all 2,5 and 9 permutations so you can't win, whatever you do.

In the second test case, the answer is 0.166666667 because you may win by getting number 952 with probability 1/6.

In the third test case the answer is 0.194444444 because you may win by getting number 952 in round1 with probability 1/6 or you can win by getting number 925 in round 1 and then 952 in round 2 with probability 1/6 * 1/6.

原博客：https://blog.csdn.net/why850901938/article/details/51204620
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
using namespace std;
int n[]={,};
int a[];
double p,q;
int main()
{
    for(int i=;i<;i++)
    {
        n[i]=n[i-]*i;
    }
    int T;
    cin>>T;
    while(T--)
    {
        char s[];
        int cnt=;
        scanf("%s",s);
        int l=strlen(s);
        for(int i=;i<l;i++)
        {
            a[i]=s[i]-'';
        }
        while(next_permutation(a,a+l))
        {
            cnt++;
        }
        if(cnt==)
        {
            printf("0.000000000\n");
        }
        else
        {
            double temp=1.0/n[l];
            q=p=temp;
            for(int i=;i<cnt;i++)
            {
                q=p+p*temp;
                p=q;
            }
            printf("%.9lf\n",p);
        }
    }
    return ;
}

F - Geometry（简单题）

Gym - 100814F

Geometry is a very important field in mathematics, Squares and rectangles are essential shapes in geometry, both of them have 4 right angles, but a square is a special case of a rectangle where width and height are the same.

The figure below shows a square on the left and a rectangle on the right:

If you have the width and the height of a 4 right angled shape, can you figure out if it is a square or a rectangle?

Input

The first line contains T, the number of test cases, for each test case there is a line with two integers (1 ≤ w, h ≤ 1, 000, 000) representing width and height, respectively.

Output

For each test case, print one line consists of 'Square'
if the shape is a square, otherwise print 'Rectangle' if it is a
rectangle.

Examples

Input

3
10 10
13 200
300 300

Output

Square
Rectangle
Square

 #include<stdio.h>
 int main()
 {
     int n,a,b;
     scanf("%d",&n);
     while(n--)
     {
         scanf("%d %d",&a,&b);
         if(a==b)
             printf("Square\n");
         else
             printf("Rectangle\n");
     }
     return ;
 }

I - Salem

Gym - 100814I

Salem is known to be one of the best competitive programmers in the region. However, he always finds a hard time understanding the concept of the hamming distance. The hamming distance of two numbers is defined as the number of different digits in their binary representations (leading zeros are used if necessary to make the binary representations have the same length). For example the hamming distance between 12 and 5 is 2 since their binary representations are 1100 and 0101 respectively and they differ in the first and fourth positions.

Recently, Salem got a problem that he needs your help to solve. He is given N integers and asked to get the maximum among the hamming distances of all possible pairs of the given integers.

Input

The first line of the input will be a single integer T representing the number of test cases. Followed by T test cases. Each test case will start with a line with single integer (2 ≤ N ≤ 100) representing the number of the integers. Each of the following N lines contains one of the integers (1 ≤ A_i ≤ 10, 000) from the list of the integers.

Output

For each test case print one line consists of one integer
representing the maximum hamming distance between all possible pairs
from the given integers.

Examples

Input

2
2
12
5
3
1
2
3

Output

2
2

 #include<iostream>
 #include<cstdio>
 using namespace std;
 int check(int n)
 {
     int c=;
     while(n)
     {
         n=n&(n-);
         c++;
     }
     return c;
 }
 int main()
 {
     int n;
     int a[];
     cin>>n;
     while(n--)
     {
         int n1;
         cin>>n1;
         int maxx=;
         for(int i=;i<=n1;i++)
         {
             cin>>a[i];
         }
         for(int i=;i<=n1;i++)
         {
             for(int  j=i+;j<=n1;j++)
             {
                 int m=check(a[i]^a[j]);
 //                cout<<"m: "<<m<<endl;
                 if(m>maxx)
                     maxx=m;
             }
         }
         cout<<maxx<<endl;
     }
 }