HDU5738(平面上的点,统计共线的点的子集)

题目链接
题目的意思很容易就能够转化成，对于给出的点集，共线的点的集合有哪些，每个大小是n的集合对答案的贡献是C(n,2)+C(n,3)+…+C(n,n)=(2^n )- n -1，比较麻烦的是重点的处理，参考题解，以及卿学姐的实现，写了一下 。

暴力统计每个点所在重点的集合大小，假设是num，那么贡献就是2^num-num-1，然后再计算与这个点共线的点集，点共线那么斜率是一样的，但是用浮点数精度比较烦，所以，化简y/x,存这两个值就好了，
然后排序，由于暴力两个for循环计算，A到B，会被计算A,B与B,A两次，所以，在统计的时候，我们跳过x值是正数的(因为我们排序是从小到大的，因为AB 斜率是正的，反过来BA斜率就是负的了，只统计某一个符号就好了)。这样就不重不漏了。
= =….菜的抠脚

#include<bits/stdc++.h>using namespace std;#define cl(a,b) memset(a,b,sizeof(a))#define LL long long#define pb push_back#define gcd __gcdconst int maxn = 1e3+200;const int inf  = 1 << 23;const LL P = 1e9+7;struct node{    int x,y;}p[maxn];vector<pair<int,int> > v;LL mod[maxn];bool vis[maxn];int main(){    mod[0]=1;mod[1]=2;    for(int i=2;i<=1000;i++)mod[i]=mod[i-1]*2%P;    int T;scanf("%d",&T);    while(T--){        int n;scanf("%d",&n);        for(int i=0;i<n;i++){            scanf("%d%d",&p[i].x,&p[i].y);        }        int num = 0;        LL ans=0;        cl(vis,0);        for(int i=0;i<n;i++){            if(vis[i])continue;            v.clear();num=0;            for(int j=0;j<n;j++){                if(p[i].x==p[j].x&&p[i].y==p[j].y){                    vis[j]=1;num++;continue;                }                int x = p[i].x - p[j].x;                int y = p[i].y - p[j].y;                int t = abs(gcd(x,y));                v.push_back(make_pair(x/t,y/t));            }            ans += (mod[num]-num-1)%P;            ans = (ans+P)%P;            sort(v.begin(),v.end());            int l=0,r=1;            while(l<v.size()){                if(v[l].first>0||(v[l].first==0&&v[l].second>0))break;                for(r=l+1;r<v.size();r++)if(v[l]!=v[r])break;                ans+=(mod[r-l]-1)*(mod[num]-1)%P;                ans=(ans+P)%P;                l=r;            }        }        printf("%lld\n",ans%P);    }    return 0;}