C20220806T1 暴力计算

2022-08-31 13:04:42 阅读：189 来源： 互联网

标签：std 暴力 rep C20220806T1 flag 计算 c11 dp define

给定一张图，按照边权走，求总边权达到 \(M\) 时用的最短长度。 \(n\leq 100,M\leq 10^{18}\) 。

首先可以用 \(dp[i][j][k]\) 表示从 \(i\) 出发通过 \(2^k\) 步走到 \(j\) 能产生的最大价值。那么有 \(dp[i][j][k]=\max\{dp[i][x][k-1]+dp[x][j][k-1]\}\) ，可以用floyd来转移，然后再反向跳 \(dp\) 数组即可得到答案。

#include<bits/stdc++.h>
#define ll long long
#define pb push_back
#define mp std::make_pair
#define pii std::pair<int,int>
#define chkmin(_A,_B) (_A=std::min(_A,_B))
#define chkmax(_A,_B) (_A=std::max(_A,_B))
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
#define copy(x,t) memcpy(x,t,sizeof(x))
//--------------------FastInOut--------------------//
class IO{
    public:
        inline char read(){
            static const int IN_LEN =1<<18|1;
            static char buf[IN_LEN],*s,*t;
            return (s==t)&&(t=(s=buf)+fread(buf, 1, IN_LEN, stdin)),(s==t)?-1:*s++;
        }
        template<typename _Tp>inline IO &operator >>(_Tp &x){
            static char c11, boo;
            for (c11=read(),boo=0;!isdigit(c11);c11=read()){
                if (c11==-1)
                    return *this;
                boo|=(c11=='-');
            }
            for(x=0;isdigit(c11);c11=read())
                x=x*10+(c11^'0');
            if(boo)
                x=-x;
            return *this;
        }
        inline void push(const char &c){
            putchar(c);
        }
        template<typename _Tp>inline IO &operator <<( _Tp x){
            if (x<0)
                x=-x,push('-');
            static _Tp sta[35];
            _Tp top=0;
            do{
                sta[top++]=x%10,x/=10;
            }while(x);
            while(top)
                push(sta[--top]+'0');
            return *this;
        }
        inline IO &operator <<(char lastChar){
            push(lastChar);
            return *this; 
        }
}FIO;
#define INF 1e18
ll f[65][205][205],rec[205],tmp[205];
ll T,n,m;
int main(){
    FIO>>T;
    while(T--){
        FIO>>n>>m;
        int r=0;
        ll ans=0;
        rep(i,1,n) rep(j,1,n){
            FIO>>f[0][i][j];
            f[0][i][j]=(!f[0][i][j])?-INF:f[0][i][j];
         }
        rep(w,1,60){
            rep(i,1,n) rep(j,1,n){
                f[w][i][j]=-INF;
                rep(k,1,n) f[w][i][j]=std:: max(f[w][i][j],f[w-1][i][k]+f[w-1][k][j]);
             }
            bool flag=false;
            rep(i,1,n) if (f[w][1][i]>=m){
                flag=true; break;
             }
            if (flag){
                r=w; break;
             }
         }
        rep(i,1,n) rec[i]=f[r-1][1][i];
        ans+=(ll)(1ll<<(r-1));
        drp(w,r-2,0){
            rep(i,1,n){ tmp[i]=-INF;
                rep(j,1,n){
                    tmp[i]=std:: max(tmp[i],rec[j]+f[w][j][i]);
                 }
             }
            bool flag=true;
            rep(i,1,n) if (tmp[i]>=m) flag=false;
            if (flag){
                copy(rec,tmp);
                ans+=(ll)(1ll<<w);
             }
         }
         FIO<<ans+1<<'\n';
     }
    return 0;
}

标签：std,暴力,rep,C20220806T1,flag,计算,c11,dp,define
来源： https://www.cnblogs.com/zhouzizhe/p/16642690.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

C20220806T1 暴力计算