写了300分得了120 被众多低年级大佬的暴力踩成了弱智
T1 avogadro
题目大意：
3行N列 第一行为N的排列 其余两行的数属于$[1,n]$
求最少删除多少列使剩下的列中 三行排序后一样
思路：
找到第二行或第三行里没有的 然后像拓扑一样搞就行了
（写的贼丑）
1 #include<iostream>
2 #include<cstdio>
3 #include<cstring>
4 #include<cstdlib>
5 #include<cmath>
6 #include<algorithm>
7 #include<queue>
8 #include<vector>
9 #include<map>
10 #define inf 2139062143
11 #define ll long long
12 #define MAXN 200100
13 using namespace std;
14 inline int read()
15 {
16 int x=0,f=1;char ch=getchar();
17 while(!isdigit(ch)) {if(ch==‘-‘)f=-1;ch=getchar();}
18 while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();}
19 return x*f;
20 }
21 //yyc score=0
22 int n,a[MAXN][3],hsh[MAXN],cnt[MAXN][3];
23 int vis[MAXN],q[MAXN],l=1,r,ans;
24 int main()
25 {
26 freopen("avogadro.in","r",stdin);
27 freopen("avogadro.out","w",stdout);
28 n=read();int b,c;
29 for(int i=1;i<=n;i++) a[i][0]=read(),hsh[a[i][0]]=i;
30 for(int i=1;i<=n;i++) a[i][1]=read(),cnt[a[i][1]][1]++;
31 for(int i=1;i<=n;i++) a[i][2]=read(),cnt[a[i][2]][2]++;
32 for(int i=1;i<=n;i++)
33 if((!cnt[i][1]||!cnt[i][2])&&!vis[hsh[i]]) vis[hsh[i]]=1,q[++r]=hsh[i];
34 while(l<=r)
35 {
36 c=q[l++],cnt[a[c][1]][1]--,cnt[a[c][2]][2]--;
37 if(!vis[hsh[a[c][1]]]&&!cnt[a[c][1]][1]) vis[hsh[a[c][1]]]=1,q[++r]=hsh[a[c][1]];
38 if(!vis[hsh[a[c][2]]]&&!cnt[a[c][2]][2]) vis[hsh[a[c][2]]]=1,q[++r]=hsh[a[c][2]];
39 }
40 for(int i=1;i<=n;i++) if(vis[i]) ans++;
41 printf("%d",ans);
42 }
T2 connect
题目大意：
N个点($n\le9$) M条边 每条边有$p_i$的概率消失
求所有点联通的概率
思路：
阶乘状压表示集合裸题
（但是忘记这是个二维dp导致算重了，而且转换也非常菜的算错了）
1 #include<iostream>
2 #include<cstdio>
3 #include<cstring>
4 #include<cstdlib>
5 #include<cmath>
6 #include<algorithm>
7 #include<queue>
8 #include<vector>
9 #include<map>
10 #define inf 2139062143
11 #define ll long long
12 #define MAXN 380100
13 #define MOD 1000000007
14 using namespace std;
15 inline int read()
16 {
17 int x=0,f=1;char ch=getchar();
18 while(!isdigit(ch)) {if(ch==‘-‘)f=-1;ch=getchar();}
19 while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();}
20 return x*f;
21 }
22 //yyc score=0
23 int n,m,lm[12],u[100],v[100],no[MAXN],g[MAXN][12],t[12];
24 double w[100],dp[MAXN],gr[MAXN];
25 inline int calc()
26 {
27 int res=0;
28 for(int i=1;i<=n;i++) res+=(t[i]-1)*(lm[n]/lm[i]);
29 return res+1;
30 }
31 int main()
32 {
33 freopen("connect.in","r",stdin);
34 freopen("connect.out","w",stdout);
35 n=read(),m=read();
36 for(int i=1;i<=m;i++)
37 {
38 u[i]=read(),v[i]=read();
39 scanf("%lf",&w[i]);
40 }
41 for(int i=lm[0]=1;i<=n;i++) lm[i]=lm[i-1]*i;
42 for(int i=1;i<=n;i++) g[1][i]=1;
43 for(int i=2,x;i<=lm[n];i++)
44 {
45 x=n;
46 while(1)
47 {
48 g[i][x]=g[i-1][x]+1;
49 if(g[i][x]>x) g[i][x]=1,g[i][x-1]++;
50 else break;
51 x--;
52 }
53 for(int j=1;j<x;j++) g[i][j]=g[i-1][j];
54 }
55 dp[lm[n]]=1.0;
56 for(int i=1;i<=m;i++)
57 {
58 for(int j=1;j<=lm[n];j++)
59 if(dp[j])
60 {
61 for(int k=1;k<=n;k++) t[k]=g[j][k];
62 for(int k=1;k<=n;k++) if(g[j][v[i]]==t[k]||t[k]==g[j][u[i]]) t[k]=min(g[j][u[i]],g[j][v[i]]);
63 gr[calc()]+=dp[j]*(1-w[i]),gr[j]+=dp[j]*w[i];
64 }
65 for(int j=1;j<=lm[n];j++) dp[j]=gr[j],gr[j]=0;
66 }
67 printf("%.3lf\n",dp[1]);
68 }
T3 dist
题目大意：
网格图上有一些点为黑点 M次询问 每次询问一个坐标离它最近的黑点的欧几里得距离
思路：
考场上没想那么多 KD tree莽就完事了
（结果估价函数写挂了）
正解是对每行处理斜率优化推一推式子
1 #include<iostream>
2 #include<cstdio>
3 #include<cstring>
4 #include<cstdlib>
5 #include<cmath>
6 #include<algorithm>
7 #include<queue>
8 #include<vector>
9 #include<map>
10 #define inf 2139062143
11 #define ll long long
12 #define MAXN 2500100
13 using namespace std;
14 inline int read()
15 {
16 int x=0,f=1;char ch=getchar();
17 while(!isdigit(ch)) {if(ch==‘-‘)f=-1;ch=getchar();}
18 while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();}
19 return x*f;
20 }
21 //yyc score=0
22 int n,tot,dim,m;
23 char ch[2010][2010];
24 struct node {int d[2],mn[2],mx[2],l,r;};
25 bool cmp (const node &a,const node &b) {return a.d[dim]<b.d[dim];}
26 inline int sqr(int x) {return x*x;}
27 struct kd_tree
28 {
29 node tr[MAXN],l,r,g;
30 int rt,ans;
31 inline void upd(int k)
32 {
33 l=tr[tr[k].l],r=tr[tr[k].r];
34 for(int i=0;i<2;i++)
35 {
36 if(tr[k].l) tr[k].mn[i]=min(tr[k].mn[i],l.mn[i]),tr[k].mx[i]=max(tr[k].mx[i],l.mx[i]);
37 if(tr[k].r) tr[k].mn[i]=min(tr[k].mn[i],r.mn[i]),tr[k].mx[i]=max(tr[k].mx[i],r.mx[i]);
38 }
39 }
40 inline int get(int k)
41 {
42 int res=0;
43 for(int i=0;i<2;i++)
44 res+=sqr(max(max(0,tr[k].mn[i]-g.d[i]),max(0,g.d[i]-tr[k].mx[i])));
45 return res;
46 }
47 inline int dis(int k)
48 {
49 int res=0;
50 for(int i=0;i<2;i++)
51 res+=sqr(tr[k].d[i]-g.d[i]);
52 return res;
53 }
54 int build(int l,int r,int D)
55 {
56 dim=D;int mid=l+r>>1;
57 nth_element(tr+l,tr+mid,tr+r+1,cmp);
58 for(int i=0;i<2;i++) tr[mid].mn[i]=tr[mid].mx[i]=tr[mid].d[i];
59 if(l<mid) tr[mid].l=build(l,mid-1,D^1);
60 if(mid<r) tr[mid].r=build(mid+1,r,D^1);
61 upd(mid);return mid;
62 }
63 int query(int k)
64 {
65 int dl=inf,dr=inf;
66 ans=min(ans,dis(k));
67 if(tr[k].l) dl=get(tr[k].l);
68 if(tr[k].r) dr=get(tr[k].r);
69 if(dl>dr)
70 {
71 if(dr<ans) query(tr[k].r);
72 if(dl<ans) query(tr[k].l);
73 }
74 else
75 {
76 if(dl<ans) query(tr[k].l);
77 if(dr<ans) query(tr[k].r);
78 }
79 }
80 }K;
81 int main()
82 {
83 freopen("dist.in","r",stdin);
84 freopen("dist.out","w",stdout);
85 n=read();
86 for(int i=0;i<n;i++)
87 {
88 scanf("%s",ch[i]);
89 for(int j=0;j<n;j++) if(ch[i][j]==‘1‘) K.tr[++tot].d[0]=i,K.tr[tot].d[1]=j;
90 }
91 n=tot,m=read();
92 K.rt=K.build(1,n,0);
93 while(m--)
94 {
95 K.g.d[0]=read()-1,K.g.d[1]=read()-1;
96 if(ch[K.g.d[0]][K.g.d[1]]==‘1‘) {puts("0");continue;}
97 K.ans=inf;K.query(K.rt);
98 printf("%d\n",K.ans);
99 }
100 }
原文地址：https://www.cnblogs.com/yyc-jack-0920/p/9903406.html
时间： 2024-10-01 10:49:48