版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/sizeof_you/article/details/82941709
猫出的题···果然毒瘤
T1
猫说是送分题，好多人满分，多判了-1挂了10分···
改后代码：
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define LL long long
using namespace std;
LL A,B,t,ans;
LL gcd(LL a,LL b){
return b==0?a:gcd(b,a%b);
}
void solve(LL a,LL b){
if(ans==-1) return;
if(a==0 || b==0){
return;
}
if(b==1) {
ans+=a; return;
}
if(a>b){
ans+=a/b;
solve(b,a%b);
}
else{
ans+=b/a;
solve(a,b%a);
}
}
int main(){
freopen("dragon.in","r",stdin);
freopen("dragon.out","w",stdout);
scanf("%lld%lld%lld",&A,&B,&t);
if(A==0) return cout<<-1,0;
solve(A,B);
if(ans==-1) puts("-1");
else printf("%lld\n",ans*t);
return 0;
}T2
模型抽象题？
可以把各种状态之间的关系看成一棵树，树根就是b-a=c-b
从根往子树走就相当于扩展，因为扩展两种一定都可以，而往父亲走就相当于收缩，因为一定只有一种方案
那么就可以看从一开始到最后恰好经过k个点的方案数
因为这是一颗无限满二叉树，所以从每个节点往下其实都是等价的
所以可以先处理出来起点状态和终点状态向上k步的情况
然后dfs往上往下跳再转移就好了
（以及注释是我自己的理解我也看不太懂std写的什么玩意 如果有错请指出）
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#define maxn 105
#define LL long long
using namespace std;
int n,m,root=-1,str,f[maxn][maxn<<1][maxn<<1],k;
const int mod=1e9+7;
typedef vector<LL> Vec;
Vec c[maxn],d[maxn],now,pre,a,b;
Vec gether(Vec x){
if(x[1]-x[0]==x[2]-x[1]) return x;
Vec ret(3);
if(x[1]-x[0]>x[2]-x[1])
ret[0]=x[0],ret[1]=2*x[1]-x[2],ret[2]=x[1];
else
ret[0]=x[1],ret[1]=2*x[1]-x[0],ret[2]=x[2];
return ret;
}
int dfs(int cur,int p,int q){//还有k步，到最高点还有p步，从lca到b还要q步
if(p<root || cur<0 || q<0) return 0;
int &res=f[cur][p][q];
if(res>=0) return res;
res=0;
if(p>q){
(res+=dfs(cur-1,p-1,q))%=mod;
(res+=dfs(cur-1,p+1,q)*2%mod)%=mod;//往下走，两个子树是等价的
}
if(p==q){//走到了lca
(res+=dfs(cur-1,p-1,q-1))%=mod;
(res+=dfs(cur-1,p+1,min(str-n+m,q+1)))%=mod;
(res+=dfs(cur-1,p+1,q))%=mod;
} return res;
}
int main(){
freopen("rabbit.in","r",stdin);
freopen("rabbit.out","w",stdout);
LL A,B,C,D,E,F;
cin>>A>>B>>C>>D>>E>>F>>k;
a.push_back(A),a.push_back(B),a.push_back(C);
b.push_back(D),b.push_back(E),b.push_back(F);
sort(a.begin(),a.end()); sort(b.begin(),b.end());
c[0]=a; d[0]=b;
for(int i=1;i<=k;i++)
c[i]=gether(c[i-1]),d[i]=gether(d[i-1]);
for(n=0;n<=k;n++){
for(m=0;m<=k;m++)
if(c[n]==d[m]) break;
if(c[n]==d[m]) break;
}
if(n+m>k) return cout<<0,0;//
now=c[n],pre=now;
for(int i=n+1;i<=k;i++){
now=gether(now);
if(now==pre) {
root=i-n-1; break;//lca到根的距离
}
pre=now;
}
if((n+m)%2!=k%2) return cout<<0,0;//
memset(f,0xcf,sizeof f);
str=(k-n-m)/2+n;//从a往上走到的最高处
if(root>=0) root=str-n-root;//从根往下走的最低处
else root=0;
f[0][str-n+m][str-n+m]=1;//初始状态
cout<<dfs(k,str,str-n)<<endl;
return 0;
}T3
一道暴力能拿92的题qwq但是我暴力好像挂了不少分
正解是分治，因为根据题的性质可以把点分成一排一排的
然后发现l~r中，可以去掉某些点让两部分变成独立的，这就是类似分治的本质qwq
虽然说spfa会死，但是这道题你可以只考虑l~r之中的，所以不会t掉
那个分治写起来还是挺复杂的···
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define N 50005
#define M 200005
#define LL long long
#define inf 0x3f3f3f3f
using namespace std;
int n,m,k,p,q,cnt,head[M];
int bel[M*10],tot;
int sum[M],c[M];
LL dis[M*10],ans[M];
bool vis[M*10];
struct EDGE{
int to,w,nxt;
}edge[M<<1];
struct Query{
int s,t,id;
}ask[M],t[M];
inline void add(int x,int y,int z){
edge[++cnt].to=y; edge[cnt].nxt=head[x];
edge[cnt].w=z; head[x]=cnt;
}
struct data{
int v;
bool operator <(const data &a) const{
return dis[v]>dis[a.v];
}
};
priority_queue <data> Q;
inline void spfa(int l,int r,int s){
for(int i=sum[l-1]+1;i<=sum[r];i++)
dis[i]=1LL<<60;
vis[s]=1,dis[s]=0,Q.push((data){s});
while(!Q.empty()){
int u=Q.top().v; vis[u]=0; Q.pop();
for(int i=head[u];i;i=edge[i].nxt){
int v=edge[i].to;
if(l>bel[v] || bel[v]>r) continue;
if(dis[v]>dis[u]+edge[i].w){
dis[v]=dis[u]+edge[i].w;
if(!vis[v]) vis[v]=1,Q.push((data){v});
}
}
}
}
void divide(int l,int r,int ql,int qr){
if(ql>qr) return;
if(r-l+1<=k+5){
for(int i=ql;i<=qr;i++){
spfa(l,r,ask[i].s);
ans[ask[i].id]=min(ans[ask[i].id],dis[ask[i].t]);
}return;
}
int ll=r+l-k>>1,rr=ll+k+1;
int qll=ql,q1,q2;
for(int i=ql;i<=qr;i++) t[i]=ask[i];
for(int i=ql;i<=qr;i++)
if(bel[t[i].s]<=ll && bel[t[i].t]<=ll) ask[qll++]=t[i];
q1=qll-1;
for(int i=ql;i<=qr;i++)
if(bel[t[i].s]>=rr && bel[t[i].t]>=rr) ask[qll++]=t[i];
q2=qll-1;
for(int i=ql;i<=qr;i++)
if(!(bel[t[i].s]<=ll&&bel[t[i].t]<=ll || bel[t[i].s]>=rr&&bel[t[i].t]>=rr))
ask[qll++]=t[i];
for(int i=sum[ll]+1;i<=sum[rr-1];i++){
spfa(l,r,i);
for(int j=ql;j<=qr;j++)
ans[ask[j].id]=min(ans[ask[j].id],dis[ask[j].s]+dis[ask[j].t]);
}
divide(l,ll,ql,q1);
divide(rr,r,q1+1,q2);
}
inline int rd(){
int x=0,f=1; char c=' ';
while(c<'0' || c>'9') f=c=='-'?-1:1,c=getchar();
while(c>='0' && c<='9') x=x*10+c-'0',c=getchar();
return x*f;
}
int main(){
freopen("transport.in","r",stdin);
freopen("transport.out","w",stdout);
n=rd(),m=rd(),k=rd(),p=rd(),q=rd();
for(int i=1;i<=n;i++){
c[i]=rd();
for(int j=1;j<=c[i];j++) bel[++tot]=i;
sum[i]=sum[i-1]+c[i];
}
memset(ans,0x3f,sizeof ans);
for(int i=1;i<=m;i++){
int x=rd(),y=rd(),z=rd();
add(x,y,z); add(y,x,z);
}
for(int i=1;i<=q;i++){
ask[i].s=rd(),ask[i].t=rd(); ask[i].id=i;
}
divide(1,n,1,q);
for(int i=1;i<=q;i++){
if(ans[i]>inf) puts("-1");
else printf("%lld\n",ans[i]);
}
return 0;
}
//5 11 1 2 5
//2 2 2 2 2
//1 2 1
//1 3 2
//2 3 1
//3 4 1
//3 5 2
//4 5 1
//5 6 1
//5 7 2
//6 7 1
//7 8 1
//9 10 1
//1 4
//2 8
//5 7
//9 4
//3 2