题意:二维平面上右一点集(S),共(n)个元素,开始位于平面上任意点(P),(P)不一定属于(S),每次操作为选一条至少包含(S)中两个元素和当前位置(P)的直线,每条直线选取概率相同,同一直线上每个点(Q in S) 选取概率相同,(Q)次询问 包含两个元素(t,m) 即点(P)到(t)共操作(m)次的最大概率
打了场(CF) 结果(D)题死活调不出来 只能一大早来补题了
可以想到记录(f[i][j][k])表示从点(i)到点(j)走(k)步的概率 这个过程我们可以通过记录(2^x)的矩阵来存储
之后可以发现对于一个询问(t,m) 我们可以通过矩阵的转移得到走(m-1)步的答案 之所以不能直接走(m)步是因为第一步的点P不一定在(S)内 分析一下就可以发现 一条线(l)上的点的概率就(frac {sum_{(i in S)}probility[i]}{sum_{(i in S)}1}) 对所有直线取个(max)就是答案了
复杂度 (O((n + q) cdot n^2 cdot log m))
#include<bits/stdc++.h>
using namespace std;
#define FO(x) {freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);}
#define pa pair<int,int>
#define mod 1000000007
#define ll long long
#define mk make_pair
#define pb push_back
#define fi first
#define se second
#define cl(x) memset(x,0,sizeof x)
#ifdef Devil_Gary
#define bug(x) cout<<(#x)<<" "<<(x)<<endl
#define debug(...) fprintf(stderr, __VA_ARGS__)
#else
#define bug(x)
#define debug(...)
#endif
const int INF = 0x7fffffff;
const int N=2e2+5;
/*
char *TT,*mo,but[(1<<15)+2];
#define getchar() ((TT==mo&&(mo=(TT=but)+fread(but,1,1<<15,stdin),TT==mo))?-1:*TT++)//*/
inline int read(){
int x=0,rev=0,ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')rev=1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return rev?-x:x;
}
int n,flg,x[N],y[N];
double ans,pro[N],tmp[N];
vector<vector<int> > lines;
struct matrix{
double g[N][N];
matrix operator * (const matrix&a){
matrix c;cl(c.g);
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
for(int k=0;k<n;k++)
c.g[i][j]+=g[i][k]*a.g[k][j];
return c;
}
}f[15];
pa fix(pa a){
if(a.fi==0) return mk(0,1);
if(a.se==0) return mk(1,0);
int d=__gcd(a.fi,a.se);
a.fi/=d,a.se/=d;
if(a.fi<0) a.fi*=-1,a.se*=-1;
return a;
}
void add(int p){
map<pa ,vector<int>>cnt;
for(int i=0;i<n;i++) if(i!=p){
int dx=x[i]-x[p],dy=y[i]-y[p];
cnt[fix(mk(dx,dy))].pb(i);
}
int sz=cnt.size();
for(auto u:cnt){
u.se.pb(p),flg=1;
for(auto v:u.se){
f[0].g[p][v]+=1.0/u.se.size()/sz;
if(v<p) flg=0;
}
if(flg) lines.pb(u.se);
}
}
int main(){
#ifdef Devil_Gary
freopen("in.txt","r",stdin);
#endif
n=read();
for(int i=0;i<n;i++) x[i]=read(),y[i]=read();
for(int i=0;i<n;i++) add(i);
for(int i=1;i<=14;i++) f[i]=f[i-1]*f[i-1];
for(int Q=read();Q;Q--){
int t=read()-1,m=read()-1;
for(int i=0;i<n;i++) pro[i]=0;pro[t]=1.0;
for(int i=0;i<=14;i++){
if(m&(1<<i)){
for(int j=0;j<n;j++) tmp[j]=pro[j],pro[j]=0;
for(int j=0;j<n;j++) for(int k=0;k<n;k++){
pro[j]+=tmp[k]*f[i].g[j][k];
}
}
}
ans=0;
for(auto line:lines){
double temp=0;
for(auto u:line) temp+=pro[u];
ans=max(ans,temp/line.size());
}
printf("%.12lf
",ans);
}
}