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 __gcd

const 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;
}