2020.12.16 模擬賽x+1
A. 接力比賽
跑兩遍揹包,再進行一些玄學的剪枝
程式碼
#include<cstdio> #include<algorithm> #define rg register inline int read(){ rg int x=0,fh=1; rg char ch=getchar(); while(ch<'0' || ch>'9'){ if(ch=='-') fh=-1; ch=getchar(); } while(ch>='0' && ch<='9'){ x=(x<<1)+(x<<3)+(ch^48); ch=getchar(); } return x*fh; } const int maxn=1e3+5,maxm=1e6+5; int n,m,maxa,maxb,mmax; long long fa[maxm],fb[maxm],ans=0; struct asd{ int val,cost; }b1[maxn],b2[maxn]; bool cmp(asd aa,asd bb){ return aa.val>bb.val; } int main(){ n=read(),m=read(); for(rg int i=1;i<=n;i++){ b1[i].cost=read(),b1[i].val=read(); } for(rg int i=1;i<=m;i++){ b2[i].cost=read(),b2[i].val=read(); } mmax=700; std::sort(b1+1,b1+1+n,cmp); std::sort(b2+1,b2+1+m,cmp); for(rg int i=1;i<=mmax;i++){ maxa+=b1[i].cost; } for(rg int i=1;i<=mmax;i++){ maxb+=b2[i].cost; } for(rg int i=1;i<=maxa;i++){ fa[i]=-0x3f3f3f3f3f3f3f3f; } for(rg int i=1;i<=maxb;i++){ fb[i]=-0x3f3f3f3f3f3f3f3f; } for(rg int i=1;i<=mmax;i++){ for(rg int j=maxa;j>=b1[i].cost;j--){ fa[j]=std::max(fa[j],fa[j-b1[i].cost]+b1[i].val); } } for(rg int i=1;i<=mmax;i++){ for(rg int j=maxb;j>=b2[i].cost;j--){ fb[j]=std::max(fb[j],fb[j-b2[i].cost]+b2[i].val); } } for(rg int i=1;i<=maxa;i++){ ans=std::max(ans,fa[i]+fb[i]); } printf("%lld\n",ans); return 0; }
B. 樹上競技
單獨考慮每一條邊的貢獻
設該邊一側關鍵點的個數為 \(x\) ,則另一側關鍵點的個數為 \(m-x\)
那麼最終的集合點一定會選在點比較多的那一側,因為我們要讓儘量少的點經過這條邊
設該邊一側點的個數為 \(s\) ,則另一側點的個數為 \(n-s\)
該邊的答案就是 \(\sum_{i=1}^{m-1}C_s^i \times C_{n-s}^{m-i} \times min(i,m-i)\)
帶著 \(min\) 不好處理,我們把它展開
\[\sum_{i=1}^{\frac{m-1}{2}} C_{S}^{i} C_{n-S}^{m-i} \times i + \sum_{i=1}^{\frac{m-1}{2}} C_{S}^{m-i} C_{n-S}^{i} \times i + [m\%2=0] \times C _{S}^{\frac{m}{2}} \times C_{n-S}^{\frac{m}{2}} \times \frac{m}{2}\]前兩項本質是相同的,最後一項加的時候特判一下即可
變一下形式\(\sum_{i=1}^{\frac{m-1}{2}} C_{S}^{i} C_{n-S}^{m-i} \times i=\sum_{i=1}^{\frac{m-1}{2}} C_{S-1}^{i-1} C_{n-S}^{m-i} \times s\)
令 \(k=\frac{m-1}{2}\)
則 \(G(S) = \sum_{i=1}^{k} C_{S-1}^{i-1} C_{n-S}^{m-i}\)
所以我們只需要求出 \(G(s)\) 即可
我們現在考慮 \(G(s)\) 的組合意義,可以理解為 \(n-1\) 個物品裡一共要選 \(m-1\)
要從 \(G(i)\) 變成 \(G(i+1)\),發現答案變少的部分就是前面 \(i-1\) 個選了 \(k-1\) 個,而 \(i\) 也被選中了
只有這種情況會被 \(G(i)\) 計算而不會被 \(G(i+1)\) 計算
\(O(n)\) 遞推即可計算,最後列舉每一條邊算貢獻再累加即可
程式碼
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define rg register
inline int read(){
rg int x=0,fh=1;
rg char ch=getchar();
while(ch<'0' || ch>'9'){
if(ch=='-') fh=-1;
ch=getchar();
}
while(ch>='0' && ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);
ch=getchar();
}
return x*fh;
}
const int mod=1e9+7,maxn=1e6+5;
int h[maxn],tot=1,n,m;
struct asd{
int to,nxt;
}b[maxn<<1];
void ad(rg int aa,rg int bb){
b[tot].to=bb;
b[tot].nxt=h[aa];
h[aa]=tot++;
}
int siz[maxn],ans;
void dfs(rg int now,rg int lat){
siz[now]=1;
for(rg int i=h[now];i!=-1;i=b[i].nxt){
rg int u=b[i].to;
if(u==lat) continue;
dfs(u,now);
siz[now]+=siz[u];
}
}
int ny[maxn],jc[maxn],jcc[maxn],fa[maxn];
int getC(rg int nn,rg int mm){
return 1LL*jc[nn]*jcc[mm]%mod*jcc[nn-mm]%mod;
}
void pre(){
ny[1]=1;
for(rg int i=2;i<maxn;i++){
ny[i]=1LL*(mod-mod/i)*ny[mod%i]%mod;
}
jc[0]=jcc[0]=1;
for(rg int i=1;i<maxn;i++){
jc[i]=1LL*jc[i-1]*i%mod;
jcc[i]=1LL*jcc[i-1]*ny[i]%mod;
}
}
int anss[maxn];
int js(rg int aa,rg int bb){
if(aa>bb) std::swap(aa,bb);
if(anss[aa]!=-1) return anss[aa];
rg int nans=0;
for(rg int i=1;i<=m-1;i++){
rg int j=m-i;
if(aa<i) break;
if(bb<j) continue;
nans+=1LL*getC(aa,i)*getC(bb,j)%mod*std::min(i,j)%mod;
if(nans>=mod) nans-=mod;
}
return anss[aa]=nans;
}
void dfs2(rg int now,rg int lat){
for(rg int i=h[now];i!=-1;i=b[i].nxt){
rg int u=b[i].to;
if(u==lat) continue;
ans+=js(siz[u],n-siz[u]);
if(ans>=mod) ans-=mod;
dfs2(u,now);
}
}
int getsum(int num){
return 1LL*num*(num+1)/2LL%mod;
}
bool haha=1;
int main(){
memset(h,-1,sizeof(h));
memset(anss,-1,sizeof(anss));
n=read(),m=read();
for(rg int i=2;i<=n;i++){
fa[i]=read();
if(fa[i]!=i-1) haha=0;
ad(fa[i],i);
ad(i,fa[i]);
}
pre();
if(n<1000000 && m>=200000) {
rg int cs;
if(m!=1){
pre();
for(rg int i=1;i<=n;i++){
if(i-1<m/2 || n-i<m/2) continue;
cs=1LL*getC(i-1,m/2)*getC(n-i,m/2)%mod;
ans+=1LL*cs*(m/2)%mod*getsum(i-1)%mod*ny[i-1]%mod;
if(ans>=mod) ans-=mod;
ans+=1LL*cs*(m/2)%mod*getsum(n-i)%mod*ny[n-i]%mod;
if(ans>=mod) ans-=mod;
}
}
} else {
dfs(1,0);
dfs2(1,0);
}
printf("%d\n",ans);
return 0;
}
C. 虛構推理
可以列舉每一個時刻,然後取最小值
但是秒不一定是整數,所以直接列舉精度不大好處理
於是可以模擬退火
把時針和分針提前排好序可以做到 \(log(n)\) 計算答案,可以多做不少遍模擬退火
程式碼比較玄學,多交幾遍可以 \(AC\)
程式碼
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<ctime>
#define rg register
const int maxn=5e4+5;
typedef double db;
const db eps=1e-6,xs=0.999;
struct asd{
int a,b,c;
}jl[maxn];
struct jie{
db jd;
int id;
}orza[maxn],orzb[maxn];
bool cmp(jie aa,jie bb){
return aa.jd<bb.jd;
}
int n,jla,jlb;
db jlc,ans=1000;
db jsfz(rg db a1,rg db b1,rg db c1,rg db a2,rg db b2,rg db c2){
rg db ans1=(db)b1*6.0+(db)c1/10.0;
rg db ans2=(db)b2*6.0+(db)c2/10.0;
return fabs(ans1-ans2);
}
db jssz(rg db a1,rg db b1,rg db c1,rg db a2,rg db b2,rg db c2){
rg db ans1=a1*30+(db)b1/2.0+(db)c1/120.0;
rg db ans2=a2*30+(db)b2/2.0+(db)c2/120.0;
return fabs(ans1-ans2);
}
int ef1(db aa){
rg int l=1,r=n,mids;
if(aa>=360) aa-=360;
if(aa<0) aa+=360;
while(l<=r){
mids=(l+r)>>1;
if(orza[mids].jd>=aa) r=mids-1;
else l=mids+1;
}
return l;
}
int ef2(db aa){
rg int l=1,r=n,mids;
if(aa>=360) aa-=360;
if(aa<0) aa+=360;
while(l<=r){
mids=(l+r)>>1;
if(orzb[mids].jd>=aa) r=mids-1;
else l=mids+1;
}
return l;
}
db jsa(rg int now,rg db aa,rg db bb,rg db cc){
rg int id;
if(now<1 || now>n) return 0;
rg db ans1,ans2;
id=orza[now].id;
ans1=jsfz(aa,bb,cc,jl[id].a,jl[id].b,jl[id].c);
ans2=jssz(aa,bb,cc,jl[id].a,jl[id].b,jl[id].c);
ans1=std::min(ans1,360.0-ans1);
ans2=std::min(ans2,360.0-ans2);
return std::max(ans1,ans2);
}
db jsb(rg int now,rg db aa,rg db bb,rg db cc){
rg int id;
if(now<1 || now>n) return 0;
rg db ans1,ans2;
id=orzb[now].id;
ans1=jsfz(aa,bb,cc,jl[id].a,jl[id].b,jl[id].c);
ans2=jssz(aa,bb,cc,jl[id].a,jl[id].b,jl[id].c);
ans1=std::min(ans1,360.0-ans1);
ans2=std::min(ans2,360.0-ans2);
return std::max(ans1,ans2);
}
db js(rg db aa,rg db bb,rg db cc){
rg db nans=0,jd1=jssz(aa,bb,cc,0,0,0),jd2=jsfz(aa,bb,cc,0,0,0);
rg int wz=ef1(jd1+180.0);
nans=std::max(nans,jsa(wz-1,aa,bb,cc));
nans=std::max(nans,jsa(wz,aa,bb,cc));
nans=std::max(nans,jsa(1,aa,bb,cc));
nans=std::max(nans,jsa(n,aa,bb,cc));
wz=ef2(jd2+180.0);
nans=std::max(nans,jsb(wz-1,aa,bb,cc));
nans=std::max(nans,jsb(wz,aa,bb,cc));
nans=std::max(nans,jsb(1,aa,bb,cc));
nans=std::max(nans,jsb(n,aa,bb,cc));
return nans;
}
void mnth(){
db t=1.0;
rg int aa=12,bb=60;
while(t>eps){
if(aa) aa--;
if(bb) bb--;
rg int na=(jla+aa*(rand()%2-1))%12;
rg int nb=(jlb+bb*(rand()%2-1))%60;
db nc=jlc+(rand()*2-RAND_MAX)%60*t;
if(na<0) na+=12;
if(nb<0) nb+=60;
if(nc<0) nc+=60;
if(nc>=60) nc-=60;
db nw=js(na,nb,nc);
db cha=nw-ans;
if(cha<0){
jla=na,jlb=nb,jlc=nc,ans=nw;
} else if(exp(-cha/t)*RAND_MAX>rand()){
jla=na,jlb=nb,jlc=nc;
}
t*=xs;
}
}
int main(){
srand(time(0));
scanf("%d",&n);
rg char ch1,ch2;
for(rg int i=1;i<=n;i++){
scanf("%02d%c%02d%c%02d",&jl[i].a,&ch1,&jl[i].b,&ch2,&jl[i].c);
jl[i].a%=12;
orza[i].id=orzb[i].id=i;
orza[i].jd=jssz(jl[i].a,jl[i].b,jl[i].c,0,0,0);
orzb[i].jd=jsfz(jl[i].a,jl[i].b,jl[i].c,0,0,0);
}
std::sort(orza+1,orza+1+n,cmp);
std::sort(orzb+1,orzb+1+n,cmp);
ans=js(jla,jlb,jlc);
for(rg int i=1;i<=350;i++){
mnth();
}
printf("%.8f\n",ans);
return 0;
}
D. 記憶碎片
考慮 \(DP\)
設 \(f[i][s]\) 表示當前考慮了前 \(i\) 短的邊,聯通狀態為 \(s\) 的方案數,\(s\) 的表示方法可以使用最小表示法(其實就是\(sort\)之後拿\(vector\)存一下再用 \(map\) 對映)
狀態數不是很多,最多隻有 \(37338\)
轉移時只需要分兩種情況
1、當前的邊不是最小生成樹中的邊,此時只能從某一個聯通塊中選兩個沒有邊的點去連
2、當前的邊是最小生成樹中的邊,考慮當前邊連線了哪兩個之前沒有相連的聯通塊即可
常數巨大
程式碼
#include<cstdio>
#include<algorithm>
#include<vector>
#include<map>
#define rg register
const int maxn=45,maxm=40005,mod=1e9+7;
int n,f[maxn*maxn][maxm],tot,sum[maxm],m;
bool vis[maxn*maxn];
std::vector<int> g[maxm],cs;
std::map<std::vector<int>,int> mp;
void dfs(rg int now,rg int lat){
if(lat==0){
if(now==0){
g[++tot]=cs;
for(rg int i=0;i<g[tot].size();i++){
rg int u=g[tot][i];
sum[tot]+=1LL*u*(u-1)/2%mod;
if(sum[tot]>=mod) sum[tot]-=mod;
}
}
return;
}
dfs(now,lat-1);
if(now>=lat){
cs.push_back(lat);
dfs(now-lat,lat);
cs.pop_back();
}
}
int main(){
scanf("%d",&n);
m=(n-1)*n/2;
rg int aa;
for(rg int i=1;i<n;i++){
scanf("%d",&aa);
vis[aa]=1;
}
dfs(n,n);
for(rg int i=1;i<=tot;i++){
std::reverse(g[i].begin(),g[i].end());
mp[g[i]]=i;
}
f[0][1]=1;
for(rg int i=0;i<m;i++){
for(rg int j=1;j<=tot;j++){
if(f[i][j]){
if(vis[i+1]){
for(rg int k=0;k<g[j].size();k++){
for(rg int o=k+1;o<g[j].size();o++){
cs.clear();
cs.push_back(g[j][k]+g[j][o]);
for(rg int p=0;p<g[j].size();p++){
if(p!=k && p!=o){
cs.push_back(g[j][p]);
}
}
std::sort(cs.begin(),cs.end());
f[i+1][mp[cs]]+=1LL*f[i][j]*g[j][k]%mod*g[j][o]%mod;
if(f[i+1][mp[cs]]>=mod) f[i+1][mp[cs]]-=mod;
}
}
} else {
f[i+1][j]+=1LL*f[i][j]*(sum[j]-i)%mod;
if(f[i+1][j]>=mod) f[i+1][j]-=mod;
}
}
}
}
printf("%d\n",f[m][tot]);
return 0;
}