BZOJ 1800: [Ahoi2009]fly 飞行棋 (尺取+数学)
1800: [Ahoi2009]fly 飞行棋
Time Limit: 10 Sec Memory Limit: 64 MB Submit: 1530 Solved: 1220 [Submit][Status][Discuss]
Description
给出圆周上的若干个点,已知点与点之间的弧长,其值均为正整数,并依圆周顺序排列。 请找出这些点中有没有可以围成矩形的,并希望在最短时间内找出所有不重复矩形。
Input
第一行为正整数N,表示点的个数,接下来N行分别为这N个点所分割的各个圆弧长度
Output
所构成不重复矩形的个数
Sample Input
8 1 2 2 3 1 1 3 3
Sample Output
3
HINT
N<= 20

思路:一开始的想法是枚举四条边,如果a==c&&b==d,就认为是找到了一个矩形。
重点在于判重,非常不好判断。一开始想要根绝四条边的长度hash一下,但是设想一个所有弧长都相等,且弧长较多的情况。
此时所有矩形的边长长度都相同,但是显然是不同的矩形,因此这种判断方法是错误的。
比较棒的思路是:矩形的对角线一定是外接圆的直径。
因此只需要找有多少条执行,假设为c,那么答案就是c*(c-1)/2(因为任意两条对角线就可以构成一个矩形)
这种做法的好处很明显,不需要判重。
至于如何找直径,直径是把圆周等分的弦,所以尺取找一下就好了。
复杂度O(n)
1/* ***********************************************
2Author :111qqz
3Created Time :2017年01月26日 星期四 18时43分21秒
4File Name :code/bzoj/r1800.cpp
5************************************************ */
6
7#include <cstdio>
8#include <cstring>
9#include <iostream>
10#include <algorithm>
11#include <vector>
12#include <queue>
13#include <set>
14#include <map>
15#include <string>
16#include <cmath>
17#include <cstdlib>
18#include <ctime>
19#define fst first
20#define sec second
21#define lson l,m,rt<<1
22#define rson m+1,r,rt<<1|1
23#define ms(a,x) memset(a,x,sizeof(a))
24typedef long long LL;
25#define pi pair < int ,int >
26#define MP make_pair
27
28using namespace std;
29const double eps = 1E-8;
30const int dx4[4]={1,0,0,-1};
31const int dy4[4]={0,-1,1,0};
32const int inf = 0x3f3f3f3f;
33const int N=25;
34int n,m;
35int a[N];
36int sum=0;
37LL ans=0;
38void ruler()
39{
40 int head;
41 int tail;
42 int cur = 0 ;
43 head = tail = 1;
44 while (head<n)
45 {
46 cur+=a[head];
47 while (cur>sum) cur-=a[tail++];
48 if (cur==sum)
49 {
50 ans++;
51 }
52 head++;
53 }
54
55}
56int main()
57{
58 #ifndef ONLINE_JUDGE
59// freopen("code/in.txt","r",stdin);
60 #endif
61
62 scanf("%d",&n);
63 for ( int i = 1 ;i <= n ; i++)
64 {
65 scanf("%d",&a[i]);
66 sum = sum + a[i];
67 }
68 //cout<<"sum:"<<sum<<endl;
69 if (sum%2==1)
70 {
71 puts("0");
72 return 0;
73 }
74 sum/=2;
75 ruler();
76 //cout<<"sum:"<<sum<<" ans:"<<ans<<endl;
77 ans=ans*(ans-1)/2;
78 printf("%lld\n",ans);
79 #ifndef ONLINE_JUDGE
80 fclose(stdin);
81 #endif
82 return 0;
83}