杂题

HH去散步
矩阵乘法优化dp
对边dp
dp[i]=西格玛dp[j];
最后找是连接节点的边，计数

// luogu-judger-enable-o2
#include<cstdio>
#include<cstring>
using namespace std;
const int M=121,p=45989; int u[M],v[M],n,m,t,S,T,st[M],end[M],ans;
struct node {
    int a[M][M];node(){memset(a,0,sizeof a);}
    node operator *(node b) const{node ans;
        for(int i=0;i<2*m;i++) for(int j=0;j<2*m;j++)
        for(int k=0;k<2*m;k++) ans.a[i][j]=(ans.a[i][j]+a[i][k]*b.a[k][j]%p)%p;
        return ans;
    }
}a,tmp;
int main(){
    scanf("%d%d%d%d%d",&n,&m,&t,&S,&T);
    for(int i=0;i<m;i++) scanf("%d%d",&u[i*2],&v[i*2]),u[i*2+1]=v[i*2],v[i*2+1]=u[i*2];
    for(int i=0;i<2*m;i++) tmp.a[i][i]=1;
    for(int i=0;i<2*m;i++) if(u[i]==S) st[i]=1;
    for(int i=0;i<2*m;i++) for(int j=0;j<2*m;j++) 
    if(u[i]==v[j]&&(i^(j^1))) a.a[i][j]=1;t--;
    for(;t;a=a*a,t>>=1) if(t&1) tmp=tmp*a;
    for(int i=0;i<2*m;i++) for(int j=0;j<2*m;j++)end[i]=(end[i]+tmp.a[i][j]*st[j])%p;
    for(int i=0; i<2*m; i++)
        v[i]==T?ans=(ans+end[i])%p:0;
    printf("%d\n", ans);
}

super gcd
更像减损发+高精
。。。。
无语的题
还要优化有公因子2的情况

// luogu-judger-enable-o2
// luogu-judger-enable-o2
// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include<cstdio>
#include<cstring>
#include<iostream>
const int M=10001,base=100000000;
struct node{
    int a[M],len;node(){std::memset(a,0,sizeof a);len=0;}
    node operator -(node b)const{
        node ans;ans.len=std::max(len,b.len);
        for(int i=1;i<=ans.len;i++) {
            ans.a[i]+=(a[i]-b.a[i]);
            if(ans.a[i]<0) ans.a[i+1]--,ans.a[i]+=base;
        }
        while(!ans.a[ans.len]&&ans.len) ans.len--;
        return ans;
    }
    node operator /(int b)const{
        node ans;ans.len=len;
        for(int i=ans.len;i>=1;i--) {
            if(a[i]&1) ans.a[i-1]+=base;
            ans.a[i]=(ans.a[i]+a[i])/b;
        }
        while(!ans.a[ans.len]&&ans.len) ans.len--;
        return ans;
    }
     node operator *(int b)const{
        node ans;ans.len=len+1;
        for(int i=1;i<=ans.len;i++) {
            ans.a[i+1]+=(ans.a[i]+1ll*a[i]*b)/base;
            ans.a[i]=(ans.a[i]+1ll*a[i]*b)%base;
        }
        while(!ans.a[ans.len]) ans.len--;
        return ans;
    }
    bool operator <(node b){
    	if(len<b.len) return 1;
    	if(len>b.len) return 0;
    	for(int i=len;i>=1;i--) if(a[i]<b.a[i]) return 1;else if(a[i]>b.a[i]) return 0;
    }
    void print(){printf("%d",a[len]);
        for(int i=len-1;i>=1;i--) printf("%08d",a[i]);
    }
    bool check(){return (len&&a[1]%2==0);}
}a,b,ans;char s[M],t[M];int cnt;
node gcd(node x,node y){
    while((y.len!=0)){
        if(x<y) std::swap(x,y);//x.print();puts("");y.print();puts(";");
        node tmp=x;if(x.check()&&y.check()) cnt++,x=x/2,y=y/2;
        else x=(y.check()?(y/2):y),y=((tmp-y).check())?((tmp-y)/2):(tmp-y);
    }
    return x;
}
int main(){
    std::cin>>s>>t;int ls=strlen(s)-1,lt=strlen(t)-1;
    for(int i=ls,tmp=1;i>=0;i--){
    	if(tmp==base) tmp=1;
    	a.a[(ls-i)/8+1]+=tmp*(s[i]-'0');tmp*=10;
    }a.len=ls/8+1;
    for(int i=lt,tmp=1;i>=0;i--){
    	if(tmp==base) tmp=1;
    	b.a[(lt-i)/8+1]+=tmp*(t[i]-'0');tmp*=10;
    }b.len=lt/8+1;ans=gcd(a,b);
    for(;cnt--;)ans=ans*2;
    ans.print();
}

P1783 海滩防御
二分答案+bfs验证是否成功，1-n是否联通，也可以并查集

// luogu-judger-enable-o2
#include<cstdio>
#include<cmath>
#include<cstring>
#include<queue> 
#include<algorithm>
using namespace std;
const int M=640001;
int n,m,head[M],nex[M*2],to[M*2],tot,L,R,x[801],y[801];
bool vis[801];
double cost[M*2];
double dist(int i,int j){
    return sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}
void add(int x,int y,double z){
    nex[++tot]=head[x];cost[tot]=z;
    to[tot]=y;head[x]=tot;
}
queue<int>q;
bool check(double mid){
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=m;i++){
        if(vis[i]) continue;
        while(!q.empty())q.pop();
        q.push(i);L=R=x[i];
        while(!q.empty()){
            int w=q.front();q.pop();
            for(int i=head[w],tmp;i;i=nex[i]){
                if(cost[i]<=mid&&!vis[tmp=to[i]]){
                    q.push(tmp);vis[tmp]=1;
                    if(x[tmp]<L) L=x[tmp];if(x[tmp]>R) R=x[tmp];
                    if(L<=mid&&R+mid>=n) return 1;
                }
            }
        }
    }
    return 0;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;++i)scanf("%d%d",&x[i],&y[i]);
    for(int i=1;i<m;++i)
        for(int j=i+1;j<=m;++j)
            add(i,j,dist(i,j)),add(j,i,dist(i,j));
    double l=1,r=100001;
    while(l+0.01<r){double mid=(l+r)/2;
    	if(check(mid)) r=mid;
    	else l=mid;
    }
        printf("%.2lf",r/2);
        return 0;
    }