大佬：這些題都太水啦，隨隨便便就AK啦
我：BJ考的都是些什么神仙題啊？？？
聽說是結論題
我推了一個結論然后YY了一個線段樹做法打了200+棄療了打了個暴力，最后發現結論有點點問題還好只打了暴力？
怎么寫啊，我不會．
聽說又是結論題
我不會呀
到底有什么結論呀．．
最后來一道lca裸題送溫暖嗎，可惜我沒領出題人這個情，預處理從2開始,然后GG了．
1 //Achen
2 #include<algorithm>
3 #include<iostream>
4 #include<cstring>
5 #include<cstdlib>
6 #include<vector>
7 #include<cstdio>
8 #include<queue>
9 #include<cmath>
10 #include<set>
11 #define For(i,a,b) for(int i=(a);i<=(b);i++)
12 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
13 const int mod=998244353,N=300007; 14 typedef long long LL; 15 typedef double db; 16 using namespace std; 17 int n,m; 18 LL kf[N][51],pr[N][51]; 19
20 template<typename T> void read(T &x) { 21 char ch=getchar(); x=0; T f=1; 22 while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar(); 23 if(ch=='-') f=-1,ch=getchar(); 24 for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f; 25 } 26
27 int ecnt,fir[N],nxt[N<<1],to[N<<1]; 28 void add(int u,int v) { 29 nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v; 30 nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u; 31 } 32
33 int R[N],f[N][20]; 34 void dfs(int x,int fa) { 35 f[x][0]=fa; 36 For(i,1,18) f[x][i]=f[f[x][i-1]][i-1]; 37 for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) { 38 R[to[i]]=R[x]+1; 39 dfs(to[i],x); 40 } 41 } 42
43 int lca(int x,int y) { 44 if(R[x]<R[y]) swap(x,y); 45 Rep(i,18,0) if(R[f[x][i]]>=R[y]) 46 x=f[x][i]; 47 if(x==y) return x; 48 Rep(i,18,0) if(f[x][i]!=f[y][i]) 49 x=f[x][i],y=f[y][i]; 50 return f[x][0]; 51 } 52
53 //#define DEBUG
54 int main() { 55 #ifdef DEBUG 56 freopen("sum.in","r",stdin); 57 freopen("sum.out","w",stdout); 58 #endif
59 read(n); 60 For(i,1,n-1) { 61 int x,y; 62 read(x); read(y); 63 add(x,y); 64 } 65 dfs(1,0); 66 For(i,1,n) kf[i][0]=1; 67 For(i,1,n) For(j,1,50) { 68 kf[i][j]=kf[i][j-1]*i%mod; 69 pr[i][j]=(pr[i-1][j]+kf[i][j])%mod; 70 } 71 read(m); 72 For(i,1,m) { 73 int x,y,k; 74 read(x); read(y); read(k); 75 int z=lca(x,y),w=z; z=f[z][0]; 76 LL ans=((pr[R[x]][k]+pr[R[y]][k])%mod-2LL*pr[R[z]][k]%mod+mod)%mod; 77 ans=(ans-kf[R[w]][k]+mod)%mod; 78 printf("%lld\n",ans); 79 } 80 return 0; 81 } 82 /*
83 5 84 1 2 85 1 3 86 2 4 87 2 5 88 2 89 1 4 5 90 5 4 45 91 */
真有意思
可以手玩兩個小數據
我不會呀
這道題難得找得到題解
但是我不想看
等sxyA了給我講
考前一天剛好被llj講到
然而我和sxy都因為沒有判-1爆成10分
列出柿子高斯消元
發現矩陣很有特點，先O(n)地消成１～n-1行的i,i+1有元素，n行的n有元素，再O(n)地消成對角線即可．
我常數太大，被洛谷卡了
llj也被洛谷卡了，但是本地過了，老張的機子真快
sxy用更優秀的寫法成功bzoj 和luogu都是rank1，但是感覺本質是一樣的，這也算常數優化嗎？總之Orz orz
1 //Achen
2 #include<algorithm>
3 #include<iostream>
4 #include<cstring>
5 #include<cstdlib>
6 #include<vector>
7 #include<cstdio>
8 #include<queue>
9 #include<cmath>
10 #include<set>
11 #define For(i,a,b) for(register int i=(a);i<=(b);i++)
12 #define Rep(i,a,b) for(register int i=(a);i>=(b);i--)
13 const int mod=1000000007; 14 typedef long long LL; 15 typedef double db; 16 using namespace std; 17 LL T,n,p,m,k; 18 LL p1,p2,C[1555],inv[1555],power1[1555],power2[1555]; 19
20 template<typename T> void read(T &x) { 21 char ch=getchar(); x=0; T f=1; 22 while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar(); 23 if(ch=='-') f=-1,ch=getchar(); 24 for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f; 25 } 26
27 inline LL ksm(LL a,LL b) { 28 LL rs=1,bs=a; 29 while(b) { 30 if(b&1) rs=rs*bs%mod; 31 bs=bs*bs%mod; 32 b>>=1; 33 } 34 return rs; 35 } 36
37 void mo(LL &x) { if(x>=mod) x-=mod; } 38
39 LL g[1555][1555]; 40 inline int solve(int n) { 41 For(i,1,n) { 42 LL tp=ksm(g[i][i],mod-2); 43 if(g[i][i]==0) return -1; 44 LL A=g[i][n+1]*tp%mod; 45 LL B=(mod-g[i][i+1])*tp%mod; 46 For(j,i+1,n) { 47 mo(g[j][i+1]+=B*g[j][i]%mod); 48 mo(g[j][n+1]+=(mod-A*g[j][i]%mod)); 49 g[j][i]=0; 50 } 51 } 52 Rep(i,n,1) { 53 For(j,i+1,n) { 54 mo(g[i][n+1]+=(mod-g[j][n+1]*g[i][j]%mod)); 55 } 56 (g[i][n+1]*=ksm(g[i][i],mod-2))%=mod; 57 } 58 return 1; 59 } 60
61 //#define DEBUG
62 int main() { 63 #ifdef DEBUG 64 freopen("heal1.in","r",stdin); 65 //freopen("heal.out","w",stdout);
66 #endif
67 read(T); 68 while(T--) { 69 read(n); read(p); read(m); read(k); 70 C[0]=1; C[1]=k; inv[0]=inv[1]=1; 71 p1=ksm(m+1,mod-2); p2=p1*m%mod; 72 power1[0]=power2[0]=1; 73 power1[1]=p1; power2[1]=p2; 74 For(i,2,n+1) { 75 inv[i]=mod-mod/i*inv[mod%i]%mod; 76 power1[i]=power1[i-1]*p1%mod; 77 power2[i]=power2[i-1]*p2%mod; 78 } 79 For(i,2,n+1) { 80 if(i<=k) C[i]=C[i-1]*max(1LL,(k-i+1))%mod*inv[i]%mod; 81 else C[i]=0; 82 } 83 memset(g,0,sizeof(g)); 84 For(i,1,n) { 85 g[i][i]++; 86 LL tpp2; 87 if(i!=n) { 88 tpp2=ksm(p2,k-min((int)k,i)); 89 Rep(j,min((int)k,i),0) { 90 LL t=p2*C[j]%mod*power1[j]%mod*tpp2%mod; 91 mo(g[i][i-j]+=mod-t); tpp2=tpp2*p2%mod; 92 } 93 tpp2=ksm(p2,k-min((int)k-1,i)-1); 94 Rep(j,min((int)k-1,i),-1) { 95 LL t=p1*C[j+1]%mod*power1[j+1]%mod*tpp2%mod; 96 mo(g[i][i-j]+=mod-t); tpp2=tpp2*p2%mod; 97 } 98 } 99 else { 100 LL tpp2=ksm(p2,k-min((int)k,i)); 101 Rep(j,min((int)k,i),0) { 102 LL t=C[j]*power1[j]%mod*tpp2%mod; 103 mo(g[i][i-j]+=mod-t); tpp2=tpp2*p2%mod; 104 } 105 } 106 g[i][n+1]=1; 107 } 108 int tp=solve(n); 109 if(tp==-1) puts("-1"); 110 else printf("%lld\n",g[p][n+1]); 111 } 112 //cerr<<clock()<<endl;
113 return 0; 114 } 115 /*
116 2 117 2 1 1 1 118 2 2 1 1 119
120 2 121 20 15 54 100 122 512 112 255 55665 123 */