# Patent

## Title
ISP-friendly rate allocation for P2P applications

## Abstract

An ISP-friendly rate allocation system and method that reduces network traffic across ISP boundaries in a peer-to-peer (P2P) network, Embodiments of the system and method continuously solve a global optimization problem and dictate accordingly how much bandwidth is allocated on each connection. Embodiments of the system and method minimize load on a server in communication with the P2P network, minimize ISP-unfriendly traffic while keeping the minimum server load unaffected, and maximize peer prefetching. Two different techniques are used to compute rate allocation, including a utility function optimization technique and a minimum cost flow formulation technique. The utility function optimization technique constructs a utility function and optimizes that utility function. The minimum cost flow formulation technique generates a minimum cost flow formulation using a bipartite graph have a vertices set and an edges set. A distributed minimum cost flow formulation is solved using Lagrangian multipliers.

## Background
BACKGROUND
Content distribution over the Internet is rapidly increasing. One area of considerable growth is video-on-demand over the Internet. However, providing high-quality Internet video-on-demand with a traditional client-server model is very costly. More importantly, the ever-mounting demand is adding significant pressure onto existing server-based infrastructures (such as data centers, content distribution networks (CDNs)), which are already under heavy burden to live up to their current load. As a result, high-profile failures are not uncommon.
Fortunately, on the heel of such crisis, peer-to-peer (P2P) networks come to rescue. Indeed, Internet video streaming (both on-demand and live broadcast) using various peer-to-peer or peer-assisted frameworks has been shown to greatly reduce the dependence on infrastructure servers, as well as bypass bottlenecks between content providers and consumers. Content distribution using a P2P network greatly reduces dependence on infrastructure servers and scales up to the demand of the Internet video era. However, it has also fundamentally altered the relationship among content owners, Internet service providers (ISPs), and consumers. In particular, ISPs on one hand are spending billions to maintain and upgrade their networks in order to support the ever increasing traffic due largely to P2P. On the other hand, ISPs are also being marginalized by content owners' direct reaching to consumers. As a result, unhappy ISPs started to put up various hurdles for P2P applications (such as throttling P2P traffic or even taking active measures to deter P2P traffic). These practices, however, created a huge backlash once they were discovered and made public. ISPs now realize that it is in their best interest to work collaboratively with content providers and consumers. However, to provide to ISPs to embrace P2P networks, any solution has to include a way for P2P applications to become ISP-friendly at the protocol level.
In general, ISP-unfriendly traffic refers to all levels of traffic unfriendly to ISPS. It includes traffic crossing ISP boundaries, between post-office protocols (PoPs), or even among different neighborhoods within the same ISP and PoP. Crossing ISP boundary traffic is a particularly significant concern of all the ISP-unfriendly traffic. ISP-unfriendly traffic is detrimental to ISPs because it increases an ISP's cost of doing business. This increase in the cost of doing business for the ISP is passed on to the user in the form of higher prices to access the Internet. This makes users unhappy. Thus, it is a double-edged sword for the ISP, who wants to keep prices low to attract customers but faces users wanting to exchange content using P2P content distribution.
At least two broad types of techniques have been used in an attempt to make P2P application more ISP-friendly. The aim of these techniques is to ensure quality of service (QoS) to an end user, to try and use as much of the P2P resources as possible, and to minimize adverse impact to the ISP.
One type of technique addresses the network topology of the P2P network. These topology-based techniques include biased neighbor selection and can be extended to take into account multi-tier network topology, from subnets, POPs, to ISPs. The topology building mechanism can also be implicit. Another type of techniques uses rate allocation to improve ISP-friendliness. Current rate allocation techniques apply a distributed optimization framework to P2P networks. However, each of these rate allocation techniques focuses on either application layer multicast or P2P live streaming.
SUMMARY
This Summary is provided to introduce a selection of concepts in a simplified form that are further described below in the Detailed Description. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter.
Embodiments of the ISP-friendly rate allocation system and method that reduces network traffic across ISP boundaries. This is called “ISP-friendly.” Instead of working on network topology, embodiments of the ISP-friendly rate allocation system and method use rate allocation to reduce traffic across ISP boundaries. Embodiments of the ISP-friendly rate allocation system and method assume the topology of the P2P network is fixed. Embodiments of the ISP-friendly rate allocation system and method operate at the packet level. Moreover, embodiments of the ISP-friendly rate allocation system and method reduce ISP-unfriendly traffic in a P2P network. In this manner, embodiments of the ISP-friendly rate allocation system and method guarantee quality of service to an end user, use as much of the P2P network as possible, and minimize adverse impact to an ISP.
In general, embodiments of the ISP-friendly rate allocation system and method continuously solve a global optimization problem and dictate accordingly how much bandwidth is allocated on each connection. Embodiments of the ISP-friendly rate allocation system and method operate at the protocol level. Embodiments of the system and method formulate an optimization problem for rate allocation that guarantees quality of service (QoS) for consumers, reduces server, load for content owners, and reduces ISP-unfriendly traffic. Moreover, embodiments of the system and method derive a distributed solution that can be executed by each P2P client independently, while collectively achieving the desired global optimal rate allocation.
Embodiments of the ISP-friendly rate allocation system and method include three hierarchical aims. The first aim is to minimize load on a server in communication with the P2P network. This is a first-order aim. The second aim is to minimize ISP-unfriendly traffic. There are many possible rate allocations that can result in the same minimum server load. Thus, the second-order aim is to minimize ISP-unfriendly traffic, while keeping the minimum server load unaffected. A third aim is to maximize peer prefetching. When both of the above aims are met, it is possible that peers still have spare upload bandwidth. This upload bandwidth is utilized to allow peers to download faster than real-time and cache future content. Prefetching can reduce server load as well as ISP-unfriendly traffic. Hence, the third-order aim is to maximize peer prefetching while leaving the first two aims unaffected.
Embodiments of the ISP-friendly rate allocation system and method can use two different techniques to compute rate allocation. A first technique is a utility function optimization technique. This technique constructs a utility function and then optimizes that utility function using a steepest descent technique. At each step the rate allocation is updated by using an update equation that guarantees convergence. Moreover, a modified utility function is generated from an original utility function by connecting two linear components of the original utility function with a concave smooth curve. The results of the utility function optimization technique yield an optimal rate allocation.
A second technique is a minimum cost flow formulation technique. This technique formulates a linear program formulation and then reformulates it as a minimum cost flow formulation. This reformulation can be done using a bipartite graph have a vertices set and an edges set. The minimum cost flow formulation is rewritten as a distributed minimum cost flow formulation. This distributed minimum cost flow formulation is updated using a rate allocation update equation. Moreover, the distributed minimum cost flow formulation is solved using Lagrangian multipliers that have the supply node Lagrangian multipliers and the demand node Lagrangian multipliers updated regularly. In addition, any surplus upload bandwidth is allocated using a convex program formulation.
Embodiments of the ISP-friendly rate allocation system and method also include a packet-level scheduling technique implement one of the two above rate allocation techniques in a P2P environment. The packet-level scheduling technique uses a token scheme to determine which neighbor peer's request to satisfy first. Specifically, the packet-level scheduling technique defines a rate allocation interval and a token allocation interval. The optimal rate allocation is updated at each rate allocation interval using an update equation. Moreover, at each token allocation interval, the peer is supplied with a number of tokens to give out. The number of tokens is proportional to the peer's upload bandwidth. At least some of the supplied tokens are allocated to a neighbor peer using the update equation.
It should be noted that alternative embodiments are possible, and that steps and elements discussed herein may be changed, added, or eliminated, depending on the particular embodiment. These alternative embodiments include alternative steps and alternative elements that may be used, and structural changes that may be made, without departing from the scope of the invention.



DRAWINGS DESCRIPTION
Referring now to the drawings in which like reference numbers represent corresponding parts throughout:
FIG. 1 is a block diagram illustrating a general overview of embodiments of the ISP-friendly rate allocation system and method disclosed herein.
FIG. 2 is a block diagram illustrating details of an exemplary embodiment of the ISP-friendly rate allocation system and method shown in FIG. 1 implemented in a peer-to-peer networking environment.
FIG. 3 is a flow diagram illustrating the operation of embodiments of the ISP-friendly rate allocation method of the ISP-friendly rate allocation system shown in FIG. 1.
FIG. 4 is a flow diagram illustrating the detailed operation of embodiments of the utility function optimization module shown in FIG. 1.
FIG. 5 is a block diagram illustrating an exemplary embodiment of a single utility function used by embodiments of the ISP-friendly rate allocation system and method shown in FIG. 1.
FIG. 6 is a block diagram illustrating an exemplary embodiment of a modified utility function used by embodiments of the ISP-friendly rate allocation system and method shown in FIG. 1.
FIG. 7 is a flow diagram illustrating the detailed operation of embodiments of the minimum cost flow formulation module shown in FIG. 1.
FIG. 8 is a block diagram illustrating an exemplary embodiment of a bipartite graph used by embodiments of the ISP-friendly rate allocation system and method shown in FIG. 1.
FIG. 9 is a flow diagram illustrating the detailed operation of embodiments of the packet-level scheduling module shown in FIG. 1.
FIG. 10 illustrates an example of a suitable computing system environment in which embodiments of the ISP-friendly rate allocation system and method shown in FIGS. 1-9 may be implemented.



DETAILED DESCRIPTION
In the following description of embodiments of the ISP-friendly rate allocation system and method reference is made to the accompanying drawings, which form a part thereof, and in which is shown by way of illustration a specific example whereby embodiments of the ISP-friendly rate allocation system and method may be practiced. It is to be understood that other embodiments may be utilized and structural changes may be made without departing from the scope of the claimed subject matter.
I. System Overview
FIG. 1 is a block diagram illustrating a general overview of embodiments of the ISP-friendly rate allocation system and method disclosed herein. It should be noted that the implementation shown in FIG. 1 is only one of many implementations that are possible. Referring to FIG. 1, an ISP-friendly rate allocation system 100 is shown implemented on a computing device 110. It should be noted that the computing device 110 may include a single processor (such as a desktop or laptop computer) or several processors and computers connected to each other.
In general, ISP-friendly rate allocation system 100 receives a request for content from its peers in the P2P network 120, processes the request and determines an optimal rate allocation, and then outputs the requested content that is sent to the requesting peers based on the optimal rate allocation 130. It should be noted that this optimal rate may change at any given time. Therefore, the ISP-friendly rate allocation system 100 is constantly updating the optimal rate allocation.
The ISP-friendly rate allocation system 100 can determine the optimal rate allocation using at least two different techniques. These techniques include a utility optimization technique and a minimum cost flow formulation technique. The ISP-friendly rate allocation system 100 includes a utility function optimization module 140 that computes the optimal rate allocation using the utility function optimization technique. The ISP-friendly rate allocation system 100 also includes a minimum cost flow formulation module 150 that computes the optimal rate allocation using the minimum cost flow formulation technique. The ISP-friendly rate allocation system 100 also includes a packet-level scheduling module 160 for implementing the optimal rate allocation at a packet level to fill requests for content from neighboring peers.
FIG. 2 is a block diagram illustrating details of an exemplary embodiment of the ISP-friendly rate allocation system 100 and method shown in FIG. 1 implemented in a peer-to-peer networking environment. In FIG. 2, a peer-assisted distribution platform 200 is shown that supports both video-on-demand (VoD) as well as bulk data dissemination. FIG. 2 illustrates the flow of the peer-assisted VoD scenario.
Referring to FIG. 2, a client 210 discovers videos via a web interface to the content library (step 1). Special URLs with embedded information redirect the client to a directory service (or tracker 220). The client 210 downloads video metadata from the tracker 220 and also obtains a list of peers who have started watching the same video earlier and are available for sharing (step 2). The client 210 then establishes direct connections with these peers whenever possible (step 3b) or gets help from signal servers if network address translation (NAT)-traversal is required (step 3a).
In the implementation shown in FIG. 2, the client 210 retrieves content from its peers in a best effort fashion. Often, however, this cannot provide sufficient quality of service (QoS), such as a continuous high bitrate for high definition and smooth video playback). Hence, the client 210 actively monitors its QoS and adaptively retrieves content from media servers 230 as needed (life line in step 4). Such media servers 230 are located in data centers, on edge networks, or hosted by content distribution networks.
In implementation shown in FIG. 2, the tracker 220 generates the peer list based on certain criteria. In general, peers within the same autonomous system (AS) with peering relationships with the client's origin AS are favored. Note that an AS is a subnetwork under separate administrative control having a common routing policy to the Internet. This criterion helps to build an ISP-friendly topology. Additionally, peers started in closer time stamps (compared to the client's 210 own start time) are favored, This is because matching peers in such a manner helps increase upload efficiency.
II. Operational Overview
In a peer-to-peer system, be it streaming or file sharing, each peer connects to multiple other peers (its neighbor peers) simultaneously. The peer downloads and uploads content data to its neighbor peers. Compared to download capacity, a peer's upload capacity is limited and frequently is the most prominent constraint of such systems. This is because most peers are connected via DSL and cable modem (even FiOS hosts have very asymmetric access). When multiple connections to neighbors contend for the limited upload capacity (say via TCP congestion control), each connection gets a fair share of the total capacity. Assuming the round-trip time (RTT) to neighbor peers is similar, the natural contention will then result in an implicit rate allocation. This means that upload capacity is evenly divided among connections (which could be called even bandwidth allocation).
Even bandwidth allocation, however, often is not sufficient. In order to overcome this limitation, embodiments of the ISP-friendly rate allocation system 100 and method continuously solve a global optimization problem and dictate accordingly how much bandwidth is allocated on each connection. It can be shown that explicit rate allocation does significantly outperform even bandwidth allocation across various scenarios.
There are at least three aims of embodiments of the ISP-friendly rate allocation system 100 and method. First, embodiments of the ISP-friendly rate allocation system 100 and method seek to minimize server load. This is a first-order and primary aim. Second, embodiments of the ISP-friendly rate allocation system 100 and method seek to minimize ISP-unfriendly traffic. There are many possible rate allocations, which could all result in the same minimum server load. Hence, the second-order aim is to minimize ISP-unfriendly traffic while keeping the minimum server load unaffected. Finally, embodiments of the ISP-friendly rate allocation system 100 and method seek to maximize peer prefetching. When both the first and second aims have been met, then it is possible that peers still have spare upload bandwidth. This is especially true in a surplus mode. Such upload bandwidth is utilized to allow peers download faster than real-time and cache future content. This feature, called “prefetching,” greatly reduces server load as well as ISP-unfriendly traffic. Hence, the third-order aim is to maximize peer prefetching while not affecting the first two aims.
FIG. 3 is a flow diagram illustrating the operation of embodiments of the ISP-friendly rate allocation method of the ISP-friendly rate allocation system 100 shown in FIG. 1. As outlined above, the method begins determining a rate allocation that will minimize a load on a server (box 300). Next, the method minimizes network traffic across ISP boundaries while keeping the minimum server load unchanged (box 310).
The method then determines that spare upload bandwidth is available for uploading (box 320). A selection then is made as to which rate allocation technique to use (box 330). If the utility function optimization technique is selected, then the method determines the rate allocation by optimizing a utility function (box 340). The details of this technique are discussed below. If the minimum cost flow formulation technique is selected, then the method determines the rate allocation using a minimum cost flow formulation (box 350). The method then outputs an optimal rate allocation (box 360).
III. Operational Details
The operational details of embodiments of the ISP-friendly rate allocation system 100 and method now will be discussed. These embodiments include embodiments of the program modules shown in FIG. 1. The operational details of each of these programs modules as well as the two rate allocation techniques mentioned above now will be discussed.
III.A. Utility Function Optimization Module
As mentioned above, the utility function optimization module 140 uses a utility function optimization technique. This technique is used to determine the rate allocation by optimizing a utility function. An overview of this utility function optimization technique first will be discussed, followed by a detailed explanation of the technique.
FIG. 4 is a flow diagram illustrating the detailed operation of embodiments of the utility function optimization module 140 shown in FIG. 1. The operation of the module 140 begins by computing a rate allocation that yields a minimum server load (box 400). Next, the module 140 computes a rate allocation that minimizes network traffic across ISP boundaries (box 410). These computations are described in detail below.
The module 140 then constructs a utility function to be optimized using a rate allocation that yields a minimum server load and the rate allocation that minimizes network traffic across ISP boundaries (box 420). A steepest descent technique then is applied to the utility function (box 430). This is performed by initializing the rate to zero (box 440) and then updating at each step using an update equation that guarantees convergence (box 450). This is explained in detail below.
As also explained below, the module 140 modifies the utility function. More specifically, the original utility function is modified by connecting two linear components of the utility function with a concave smooth curve (box 460). The original utility function is shown in FIG. 5 and the modified utility function is shown qualitatively in FIG. 6. The module 140 then outputs a solution to the optimized utility function that represents the optimal rate allocation (box 470). Each of these steps is described in detail as follows.
III.A.1. Brief Primer
By way of example, assume there are n peers in the system at any instant of time. Denote them as Peer k (k=1, 2, . . . , n), ordered by their arrival time. When Peer j arrives, it connects to a subset of peers already in the system (say including Peer i) and requests data from them. The rate xi,j is the rate Peer i allocates from its upload capacity Ui to serve Peer j. Each peer keeps track of its total upload capacity, which can initially be estimated based on historical values and then measured and updated once data starts to flow to its neighbors. Denote Sj as the set of all peers uploading to Peer j (Peer j's upstream neighbors). Denote Di as the set of all peers downloading from Peer i (Peer i's downstream neighbors). The aggregate rate Peer j receives from all of its upstream neighbors is denoted as:

xj=ΣiεSjxi,j.

Denote Rj as the desirable streaming rate of Peer j in order to maintain smooth video playback. In general, Rj could vary based on the amount of content in its local cache. For simplicity, however, let all Rj's be the same for now, equal to the video bitrate R. In this case, smooth video playback requires xj≧R. Hence, if xj<R, then Peer i will request data from the server at rate R−xj to make up for the deficit. If xj>R, on the other hand, Peer j will download data faster than real-time and cache for future. This is called “prefetching.” Prefetching combats peer chum, bandwidth jitter, and so forth. The ability to prefetch is one of the key differentiations between video on demand (VoD) and live streaming. The aggregate rate at which Peer i uploads to all its downstream neighbors is:

ΣjεDixi,j.

Clearly, this cannot exceed Peer is upload bandwidth, represented as,

ΣjεDixi,j≦Ui.

Furthermore, for unified representation, Peer 0 is denoted as a special peer, a peer that represents the server. Then, x0,j is the rate Peer j obtained from the server, which satisfies x0,j=max(0,R−xj).

In this document, the following assumptions are made for purposes of the analysis. First it is assumed that peers cache all the content they watch and keep them until they leave the system. Second, it is assumed that peers only upload to their downstream neighbors, and not vice versa. Finally, it is assumed that peers have more content than their downstream neighbors and thus are able to upload as fast as needed.
III.A.2. The 3-Stage Optimization Process
The details of embodiments of the utility function optimization module 140 now will be described. This system and method achieve each of the aforementioned objectives through a 3-stage optimization process.
III.A.2.i. First Stage: Minimize Server Load
The first stage of embodiments of the utility function optimization module 140 minimizes the server load. This can be realized as follows:




  
    
      
        min
        ⁢
        
          
            ∑
            j
            
              
            
          
          ⁢
          
            
          
          ⁢
          
            x
            
              0
              ,
              j
            
          
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          
            s
            .
            t
            .
            
              
                ∑
                
                  i
                  ∈
                  
                    {
                    
                      0
                      ,
                      
                        D
                        j
                      
                    
                    }
                  
                
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                x
                
                  i
                  ,
                  j
                
              
            
          
          =
          R
        
        ,
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        
          
            
              ∑
              
                j
                ∈
                
                  D
                  i
                
              
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              x
              
                i
                ,
                j
              
            
          
          ≤
          
            U
            i
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        
          
            x
            
              i
              ,
              j
            
          
          ≥
          
            0
            ⁢
            
              ∀
              i
            
          
        
        ,
        
          j
          ⁢
          
            
          
          ⁢
          
            
              (
              
                i
                ,
                
                  j
                  ≠
                  0
                
              
              )
            
            .
          
        
      
    
    
      
        (
        4
        )
      
    
  




After the optimization, the minimum server load can then be computed as,

U0min=Σjx*0,j,

where

{x*0,j}j=1n 

is an optimal solution.

III.A.2.ii Second Stage: Minimize ISP-Unfriendly Traffic

The second stage of embodiments of the utility function optimization module 140 seeks to limit the upload capacity of the server to U0min from the first stage. Next, a link cost is associated with ISP-unfriendly traffic. Denote gi,j(xi,j) as the cost for Peer i to upload to Peer j at rate xi,j. The link cost can take forms as simple as the sign function, where the cost is 0 for intra-ISP and 1 for inter-ISP traffic. However, the link cost could also be generalized to incorporate various levels of cost. By way of example, these levels could include using different costs to differentiate intra-POP and inter-POP traffic within the same ISP.
The rate allocation can be represented as follows:




  
    
      
        
          min
          ⁢
          
            
              ∑
              
                i
                ,
                j
              
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              
                g
                
                  i
                  ,
                  j
                
              
              ⁡
              
                (
                
                  x
                  
                    i
                    ,
                    j
                  
                
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            s
            .
            t
            .
            
              
            
            ⁢
            
              (
              2
              )
            
          
          ,
          
            (
            3
            )
          
          ,
          
            
              (
              4
              )
            
            ⁢
            
              
            
            ⁢
            and
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          
            ∑
            j
            
              
            
          
          ⁢
          
            
          
          ⁢
          
            x
            
              0
              ,
              j
            
          
        
        ≤
        
          
            
              U
              0
              min
            
            ⁡
            
              (
              
                from
                ⁢
                
                  
                
                ⁢
                1
                ⁢
                st
                ⁢
                
                  -
                
                ⁢
                stage
              
              )
            
          
          .
        
      
    
    
      
        (
        6
        )
      
    
  




After the optimization the minimum ISP-unfriendly traffic can be computed as,

GminΣi,jgi,j(x*i,j).

III.A.2.iii. Third Stage: Maximize Peer Prefetching

The third stage of embodiments of the utility function optimization module 140 seeks to allow peers download more than their demand. Peer prefetching then can be maximized when all peers fully utilize their upload capacity. Hence, this objective is translated into minimizing the remaining of peers' upload capacity. The first two objectives will not be affected, as long the minimum server load and ISP-unfriendly traffic is limited to U0min and Gmin, respectively. In this case, the rate allocation is as follows:




  
    
      
        
          min
          ⁢
          
            
              ∑
              i
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              (
              
                
                  U
                  i
                
                -
                
                  
                    ∑
                    j
                    
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    x
                    
                      i
                      ,
                      j
                    
                  
                
              
              )
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            s
            .
            t
            .
            
              
            
            ⁢
            
              (
              3
              )
            
          
          ,
          
            (
            4
            )
          
          ,
          
            
              (
              6
              )
            
            ⁢
            
              
            
            ⁢
            and
          
        
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          
            ∑
            
              i
              ∈
              
                {
                
                  0
                  ,
                  
                    D
                    j
                  
                
                }
              
            
            
              
            
          
          ⁢
          
            
          
          ⁢
          
            x
            
              i
              ,
              j
            
          
        
        ≥
        
          R
          ⁢
          
            
          
          ⁢
          
            
              (
              
                download
                ⁢
                
                  
                
                ⁢
                could
                ⁢
                
                  
                
                ⁢
                exceed
                ⁢
                
                  
                
                ⁢
                demand
              
              )
            
            .
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          
            ∑
            
              i
              ,
              j
            
            
              
            
          
          ⁢
          
            
          
          ⁢
          
            
              g
              
                i
                ,
                j
              
            
            ⁡
            
              (
              
                x
                
                  i
                  ,
                  j
                
              
              )
            
          
        
        ≤
        
          
            
              G
              min
            
            ⁡
            
              (
              
                from
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                nd
                ⁢
                
                  -
                
                ⁢
                stage
              
              )
            
          
          .
        
      
    
    
      
        (
        9
        )
      
    
  




III.A.2.iv. Unifying the Three Stages

The 3-stage optimization then is unified into a single optimization, as follows:




  
    
      
        
          
            min
            ⁢
            
              
            
            ⁢
            α
            ⁢
            
              
                ∑
                j
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                x
                
                  0
                  ,
                  j
                
              
            
          
          +
          
            β
            ⁢
            
              
                ∑
                
                  i
                  ,
                  j
                
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                
                  g
                  
                    i
                    ,
                    j
                  
                
                ⁡
                
                  (
                  
                    x
                    
                      i
                      ,
                      j
                    
                  
                  )
                
              
            
          
          +
          
            γ
            ⁢
            
              
                ∑
                j
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  
                    U
                    x
                  
                  -
                  
                    
                      ∑
                      j
                      
                        
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      x
                      
                        i
                        ,
                        j
                      
                    
                  
                
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    i
                    ∈
                    
                      {
                      
                        0
                        ,
                        
                          D
                          i
                        
                      
                      }
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            ≥
            R
          
          ,
          
            
              
                ∑
                
                  j
                  ∈
                  
                    D
                    i
                  
                
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                x
                
                  i
                  ,
                  j
                
              
            
            ≤
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          
            j
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  i
                  ,
                  
                    j
                    ≠
                    0
                  
                
                )
              
              .
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




It can be shows that, if α>>β>>γ, then the solution to Equation (10) will be equivalent to the 3-stage optimization. Both optimizations can be carried out with the existence of a central oracle, but becomes more difficult if all peers strike distributedly and independently.
III.A.3.Utility-Based Optimization
Embodiments of the utility function optimization module 140 next introduce an artificial utility function. This function not only bears an intuitive meaning, but can also directly translate the centralized optimization into a utility-based optimization, from which a completely distributed solution is merited.
III.A.3.i. Single Utility Function
Embodiments of the utility function optimization module 140 use a utility function for each peer in terms of the aggregate rate received from all its upstream neighbors (excluding the server or Peer 0). This utility function can be represented as:




  
    
      
        
          f
          ⁡
          
            (
            
              x
              j
            
            )
          
        
        =
        
          {
          
            
              
                
                  
                    a
                    ·
                    
                      x
                      j
                    
                  
                
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        x
                        j
                      
                    
                    ≤
                    R
                  
                
              
              
                
                  
                    
                      a
                      ·
                      R
                    
                    +
                    
                      b
                      ·
                      
                        (
                        
                          
                            x
                            j
                          
                          -
                          R
                        
                        )
                      
                    
                  
                
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        x
                        j
                      
                    
                    >
                    R
                  
                
              
            
            .
          
        
      
    
    
      
        (
        11
        )
      
    
  




where a, b are positive constants and a>b.

FIG. 5 is a block diagram illustrating an exemplary embodiment of a single utility function, f(x), used by embodiments of the ISP-friendly rate allocation system and method shown in FIG. 1. Note that a represents the first positive constant and b represents the second positive constant. At the intersection where the first positive constant a and the second positive constant b meet is the streaming rate R.
Here, a represents the value of getting rate from other peers up to streaming rate R. This means that for every unit of rate from other peers increases the utility by a. Also, a can be interpreted as the cost of server bandwidth. For example, in a peer-assisted VoD session, the server has to supplement all the rate peers cannot obtain from each other. At the same time, each peer has to maintain receiving data at least at streaming rate R. Hence, up to streaming rate R, every unit of bandwidth that a peer is able to obtain from other peers is a unit of bandwidth, the server can save.
A strictly positive b represents the value of prefetching. This is a major difference between live streaming and VoD streaming. In live streaming, peers have roughly synchronous playback times. This means that peers cannot buffer up much ahead of time and there can be very little prefetching (in other words, b=0). In VoD, however, a peer's playback time can be sufficiently far apart and there is great value for a peer to buffer up when possible to save for future use.
As a matter of fact, the piece-wise linear utility function ensures that, when the total utility is maximized, the server load is minimized. In addition, peer prefetching is maximized simultaneously.
Proposition 1: The solutions to,




  
    
      
        
          
            max
            
              x
              
                i
                ,
                j
              
            
          
          ⁢
          
            
              ∑
              j
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              f
              ⁡
              
                (
                
                  x
                  j
                
                )
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    j
                    ∈
                    
                      D
                      i
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            ≤
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          j
          ,
        
      
    
    
      
        (
        12
        )
      
    
  




minimize server load and maximize peer prefetching.

Proof sketch: Minimizing server bandwidth can be described as the following optimization problem:




  
    
      
        
          
            min
            
              x
              
                i
                ,
                j
              
            
          
          ⁢
          
            
              ∑
              j
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              a
              ·
              
                max
                ⁡
                
                  (
                  
                    0
                    ,
                    
                      R
                      -
                      
                        x
                        j
                      
                    
                  
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    j
                    ∈
                    
                      D
                      i
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            ≤
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          j
          ,
        
      
    
    
      
        (
        13
        )
      
    
  




since a is a positive constant.

This is equivalent to,




  
    
      
        
          
            max
            
              x
              
                i
                ,
                j
              
            
          
          ⁢
          
            
              ∑
              j
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              a
              ·
              
                min
                ⁡
                
                  (
                  
                    
                      x
                      j
                    
                    ,
                    R
                  
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    j
                    ∈
                    
                      D
                      i
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            ≤
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          
            j
            .
          
        
      
    
    
      
        (
        14
        )
      
    
  




The utility function given in Equation (11) can be rewritten as,

f(xj)=a·min(R, xj)+b·max(0, xj−R)=b·xj+(a−b)·min(R, xj)  (15)

Each of the solutions to Equation (12) using the proposed utility function given in Equation (11) must satisfy,

ΣjεDixi,j=Ui 

in order to maximize peer prefetching. This can be shown through contradiction. Hence, Equation (12) is equivalent to,





  
    
      
        
          
            
              max
              
                x
                
                  i
                  ,
                  j
                
              
            
            ⁢
            
              
                ∑
                j
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                b
                ·
                
                  U
                  j
                
              
            
          
          +
          
            
              ∑
              j
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  a
                  -
                  b
                
                )
              
              ·
              
                min
                ⁡
                
                  (
                  
                    R
                    ,
                    
                      x
                      j
                    
                  
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    j
                    ∈
                    
                      D
                      i
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            =
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          
            j
            .
          
        
      
    
    
      
        (
        16
        )
      
    
  




The term Σjb·Uj is a constant and a−b>0. This means that Equation (16) is equivalent to,





  
    
      
        
          
            max
            
              x
              
                i
                ,
                j
              
            
          
          ⁢
          
            
              ∑
              j
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              a
              ·
              
                min
                ⁡
                
                  (
                  
                    
                      x
                      j
                    
                    ,
                    R
                  
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    j
                    ∈
                    
                      D
                      i
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            =
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          
            j
            .
          
        
      
    
    
      
        (
        17
        )
      
    
  




The only difference between Equations (14) and (17) is the constraint on,

ΣjεDixi,j, min(xj,R)=min(ΣiεSixi,j, R)

is nondecreasing in,

xi,j, ∀i,j.

Thus, the solutions to Equation (17) are a subset of those to Equation (14). In summary, the solutions to Equation (12) are a subset of those to of Equation (14) and minimize server load and maximize peer prefetching. This concludes the sketch of the proof of the above proposition.

The connection between the utility optimization given in Equation (12) and the unified 3-stage optimization given in Equation (10) now will be explained. In fact, ISP-unfriendly link cost aside (or letting gi,j(xi,j)=0), it is not difficult to show that these two optimizations are equivalent and, in addition, shown that a=α and b=γ. Compared to the unified 3-stage optimization given in Equation (10), the utility function introduced here bears a much more intuitive meaning. In fact, it is straightforward to see that with the piece-wise utility function, a peer would not help a neighbor prefetch if it could instead use the upload bandwidth to satisfy other neighbors' basic streaming demand first. Hence, minimizing server load will have a higher priority than prefetching.
III.A.3.ii. Costs on ISP-Unfriendly Traffic
Recall that in order to reduce ISP-unfriendly traffic, embodiments of the utility function optimization module 140 associate a link cost with ISP-unfriendly traffic. In other words, the cost for Peer i to upload to Peer j at rate xi,j is gi,j(xi,j). The overall utility optimization problem then becomes,




  
    
      
        
          
            max
            
              x
              
                i
                ,
                j
              
            
          
          ⁢
          
            
              ∑
              
                1
                ≤
                j
                ≤
                n
              
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              (
              
                
                  f
                  ⁡
                  
                    (
                    
                      x
                      j
                    
                    )
                  
                
                -
                
                  
                    ∑
                    
                      i
                      ∈
                      
                        S
                        j
                      
                    
                    
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      g
                      
                        i
                        ,
                        j
                      
                    
                    ⁡
                    
                      (
                      
                        x
                        
                          i
                          ,
                          j
                        
                      
                      )
                    
                  
                
              
              )
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              s
              .
              t
              .
              
                
                  ∑
                  
                    j
                    ∈
                    
                      D
                      i
                    
                  
                  
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  x
                  
                    i
                    ,
                    j
                  
                
              
            
            ≤
            
              U
              i
            
          
          ,
          
            
              x
              
                i
                ,
                j
              
            
            ≥
            
              0
              ⁢
              
                ∀
                i
              
            
          
          ,
          
            j
            .
          
        
      
    
    
      
        (
        18
        )
      
    
  




A natural choice for

gi,j(•) is gi,j(xi,j)=ci,j·xi,j,

where ci,j is a positive constant that represents the cost of getting each unit of rate from Peer i to Peer j. When positive costs are used for various ISP-unfriendly traffic, it is intuitive that the maximization in Equation (18) will reduce undesirable traffic. However, this is potentially at the cost of increasing server load. The relationship among ci,j, a and b controls the tradeoff between the server load and the ISP-friendliness. To better understand this, consider a simple case where ci,j=0 if Peer i and Peer j are within the same ISP, and ci,j=c>0 otherwise.

Proposition 2: Compared to solutions for Equation (12), solutions to Equation (18) may have higher server rates. For each additional unit of server rate used due to incorporating ISP friendliness, there is at least a reduction of (a−b)/c units of ISP-unfriendly traffic.
Proof sketch: Among all the ISP-unaware solutions that minimize server load, a subset of them requires the least ISP-unfriendly traffic. Let { xi,j} be one of these solutions. Let {{tilde over (x)}i,j} be a solution for Equation (18). Then by definition,



  
    
      
        ∑
        
          1
          ≤
          j
          ≤
          n
        
        
          
        
      
      ⁢
      
        (
        
          
            f
            ⁡
            
              (
              
                
                  x
                  _
                
                j
              
              )
            
          
          -
          
            
              ∑
              
                i
                ∈
                
                  S
                  j
                
              
              
                
              
            
            ⁢
            
              
            
            ⁢
            
              
                g
                
                  i
                  ,
                  j
                
              
              ⁡
              
                (
                
                  
                    x
                    _
                  
                  
                    i
                    ,
                    j
                  
                
                )
              
            
          
        
        )
      
    
    ≤
    
      
        ∑
        
          1
          ≤
          j
          ≤
          n
        
        
          
        
      
      ⁢
      
        
          (
          
            
              f
              ⁡
              
                (
                
                  
                    x
                    ~
                  
                  j
                
                )
              
            
            -
            
              
                ∑
                
                  i
                  ∈
                  
                    S
                    j
                  
                
                
                  
                
              
              ⁢
              
                
              
              ⁢
              
                
                  g
                  
                    i
                    ,
                    j
                  
                
                ⁡
                
                  (
                  
                    
                      x
                      ~
                    
                    
                      i
                      ,
                      j
                    
                  
                  )
                
              
            
          
          )
        
        ⁢
              
    
  




  
    
      
        ∑
        
          1
          ≤
          j
          ≤
          n
        
        
          
        
      
      ⁢
      
        (
        
          
            f
            ⁢
            
              (
              
                
                  x
                  _
                
                j
              
              )
            
          
          -
          
            f
            ⁡
            
              (
              
                
                  x
                  ~
                
                j
              
              )
            
          
        
        )
      
    
    ≤
    
      
        ∑
        
          1
          ≤
          j
          ≤
          n
        
        
          
        
      
      ⁢
      
        (
        
          
            ∑
            
              i
              ∈
              
                S
                j
              
            
            
              
            
          
          ⁢
          
            
          
          ⁢
          
            (
            
              
                
                  g
                  
                    i
                    ,
                    j
                  
                
                ⁡
                
                  (
                  
                    
                      x
                      _
                    
                    
                      i
                      ,
                      j
                    
                  
                  )
                
              
              -
              
                
                  g
                  
                    i
                    ,
                    j
                  
                
                ⁡
                
                  (
                  
                    
                      x
                      ~
                    
                    
                      i
                      ,
                      j
                    
                  
                  )
                
              
            
            )
          
        
        )
      
    
  



Using Equation (15),




  
    
      
        LHS
        =
        
        ⁢
        
          
            (
            
              
                b
                ⁢
                
                  
                    ∑
                    j
                    
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      x
                      _
                    
                    j
                  
                
              
              +
              
                
                  (
                  
                    a
                    -
                    b
                  
                  )
                
                ·
                
                  min
                  ⁡
                  
                    (
                    
                      R
                      ,
                      
                        
                          x
                          _
                        
                        j
                      
                    
                    )
                  
                
              
            
            )
          
          -
          
            (
            
              
                b
                ⁢
                
                  
                    ∑
                    j
                    
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      x
                      ~
                    
                    j
                  
                
              
              +
              
                
                  (
                  
                    a
                    -
                    b
                  
                  )
                
                ·
                
                  min
                  ⁡
                  
                    (
                    
                      R
                      ,
                      
                        
                          x
                          ~
                        
                        j
                      
                    
                    )
                  
                
              
            
            )
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            b
            ·
            
              (
              
                
                  
                    ∑
                    j
                    
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      x
                      _
                    
                    j
                  
                
                -
                
                  
                    ∑
                    j
                    
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      x
                      ~
                    
                    j
                  
                
              
              )
            
          
          +
          
            
              (
              
                a
                -
                b
              
              )
            
            ·
            
              
                (
                
                  
                    min
                    ⁡
                    
                      (
                      
                        R
                        ,
                        
                          
                            x
                            _
                          
                          j
                        
                      
                      )
                    
                  
                  -
                  
                    min
                    ⁡
                    
                      (
                      
                        R
                        ,
                        
                          
                            x
                            ~
                          
                          j
                        
                      
                      )
                    
                  
                
                )
              
              .
            
          
        
      
    
  




As stated earlier,

Σj xj=ΣjUj,

and thus,

Σj xj≧Σj{tilde over (x)}j.

The term,

min(R, xj)−min(R,{tilde over (x)}j)

is the increase in server load due to incorporating ISP-friendliness. Thus,

LHS≧(a−b)·(additional server load)

due to incorporating ISP-friendliness. In addition,





  
    
      
        
          RHS
          =
          
            
              c
              ·
              
                
                  ∑
                  
                    1
                    ≤
                    j
                    ≤
                    n
                  
                  
                    
                  
                
                ⁢
                
                  (
                  
                    
                      ∑
                      
                        
                          i
                          ∈
                          
                            S
                            j
                          
                        
                        
                          i
                          ∉
                          
                            I
                            j
                          
                        
                      
                    
                    ⁢
                    
                      
                        x
                        _
                      
                      ij
                    
                  
                  )
                
              
            
            -
            c
          
        
        ⁣
        
          ·
          
            
              ∑
              
                1
                ≤
                j
                ≤
                n
              
              
                
              
            
            ⁢
            
              (
              
                
                  ∑
                  
                    
                      i
                      ∈
                      
                        S
                        j
                      
                    
                    
                      i
                      ∉
                      
                        I
                        j
                      
                    
                  
                
                ⁢
                
                  
                    x
                    ~
                  
                  ij
                
              
              )
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




where Ij denotes the subset of Peer j's neighbors in the same ISP as Peer j.

Here,




  
    ∑
    
      1
      ≤
      j
      ≤
      n
    
  
  ⁢
  
    (
    
      
        ∑
        
          
            i
            ∈
            
              S
              j
            
          
          
            i
            ∉
            
              I
              j
            
          
        
      
      ⁢
      
        x
        
          i
          ,
          j
        
      
    
    )
  




is nothing but the total ISP-unfriendly rate using allocation {xi,j}. Hence,

RHS=c·(saving in ISP-unfriendly rate).  (21)

Combining Equation (19) and Equation (21) yields,





  
    
      
        
          saving
          ⁢
          
            
          
          ⁢
          in
          ⁢
          
            
          
          ⁢
          ISP
          ⁢
          
            -
          
          ⁢
          unfriendly
          ⁢
          
            
          
          ⁢
          rate
        
        ≥
        
          
            
              a
              -
              b
            
            c
          
          ⁢
          
            
              (
              
                additional
                ⁢
                
                  
                
                ⁢
                server
                ⁢
                
                  
                
                ⁢
                load
              
              )
            
            .
          
        
      
    
    
      
        (
        22
        )
      
    
  




Intuitively, setting a<c will result in a pure ISP-friendly solution that eliminates ISP-unfriendly traffic completely. On the other hand, setting b>c will result in full utilization of a peer's upload bandwidth at all times.
III.A.4. Distributed Solution
Now that the rate allocation problem has been formulated as a convex optimization problem with linear constraints, it is straightforward to apply classical distributed solutions. In particular, some embodiments of the utility function optimization module 140 adopt the feasible steepest descent algorithm. This means that xi,j is first initialized to 0 and then updated at each step using the update equation given as:




  
    
      
        
          
            
              x
              .
            
            
              i
              ,
              j
            
          
          =
          
            Δ
            ·
            
              (
              
                
                  
                    
                      ∂
                      
                        
                      
                    
                    
                      ∂
                      
                        x
                        
                          i
                          ,
                          j
                        
                      
                    
                  
                  ⁢
                  
                    f
                    ⁡
                    
                      (
                      
                        x
                        j
                      
                      )
                    
                  
                
                -
                
                  
                    
                      ∂
                      
                        
                      
                    
                    
                      ∂
                      
                        x
                        
                          i
                          ,
                          j
                        
                      
                    
                  
                  ⁢
                  
                    
                      g
                      
                        i
                        ,
                        j
                      
                    
                    ⁡
                    
                      (
                      
                        x
                        
                          i
                          ,
                          j
                        
                      
                      )
                    
                  
                
              
              )
            
          
        
        ,
        
          

        
        ⁢
        
          
            x
            
              i
              ,
              j
            
          
          =
          
            
              [
              
                
                  x
                  
                    i
                    ,
                    j
                  
                
                +
                
                  
                    x
                    .
                  
                  
                    i
                    ,
                    j
                  
                
              
              ]
            
            +
          
        
      
    
    
      
        (
        23
        )
      
    
  




where [•]+ means an I2 projection onto a feasible set, which guarantees convergence for this problem given in Equation [10].

By way of example, due to the simple symmetric linear constraints, this step can be easily implemented as given by the following pseudo-code:



  
    while
    ⁢
    
      
    
    ⁢
    
      (
      
        
          ∑
          
            x
            
              i
              ,
              j
            
          
        
        >
        
          U
          i
        
      
      )
    
  




  
    
      {
      
        
          N
          i
        
        =
        
          
            
              number
              ⁢
              
                
              
              ⁢
              of
              ⁢
              
                
              
              ⁢
              neighbors
              ⁢
              
                
              
              ⁢
              with
              ⁢
              
                
              
              ⁢
              
                x
                
                  i
                  ,
                  j
                
              
            
            >
            
              0
              ⁢
              
                :
              
              ⁢
              
                

              
              ⁢
              
                foreach
                ⁡
                
                  (
                  j
                  )
                
              
              ⁢
              
                

              
              ⁢
              
                if
                ⁡
                
                  (
                  
                    
                      x
                      
                        i
                        ,
                        j
                      
                    
                    >
                    0
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                x
                
                  i
                  ,
                  j
                
              
            
          
          =
          
            
              x
              
                i
                ,
                j
              
            
            -
            
              
                min
                (
                
                  
                    x
                    
                      i
                      ,
                      j
                    
                  
                  ·
                  
                    
                      
                        U
                        i
                      
                      -
                      
                        ∑
                        
                          x
                          
                            i
                            ,
                            j
                          
                        
                      
                    
                    
                      N
                      i
                    
                  
                
                )
              
              ⁢
              
                :
              
            
          
        
      
      }
    
    .
  



Note that in order for Peer i to carry out the update step, the only external information required is the aggregate received bandwidth at Peer j. This information can be easily piggybacked in Peer j's packet requests.
III.A.5. Practical Considerations
The piece-wise linear utility function f(•) and linear link cost functions gi,j(•) have shortcomings in practice. First, f(•) has a sudden slope change at exactly video rate R. This causes unsteady convergence behavior and, more importantly, eliminates ISP-unfriendly prefetching completely. As a matter of fact, moderate prefetching is beneficial in reducing the server load and improving a peer's uplink utilization in the long run, especially in the presence of peer churning, bandwidth jitter or flash crowd (when a neighboring peer's playback points are very close to each other).
Thus, in some embodiments of the utility function optimization module 140, f(•) is modified slightly by connecting the two linear components with a concave smooth curve of width ΔR. FIG. 6 is a block diagram illustrating an exemplary embodiment of a modified utility function used by embodiments of the ISP-friendly rate allocation system and method shown in FIG. 1.
The implication is that moderate ISP-unfriendly prefetching is potentially allowed, but bound the amount by ΔR per peer in the worst case. Some embodiments use ΔR=0.1 R and a Deg. 3 polynomial for the curve. Next, some embodiments change the link cost function g(•) from linear to a flat quadratic one such that the optimization is strictly concave. This allows faster and steadier convergence. The link cost function g(•) is set such that,




  
    
      
        d
        
          
        
      
      
        d
        x
      
    
    ⁢
    
      g
      ⁡
      
        (
        x
        )
      
    
  
  =
  
    
      c
      +
      
        ɛ
        ⁢
        
          
        
        ⁢
        x
        ⁢
        
          
        
        ⁢
        with
        ⁢
        
          
        
        ⁢
        ɛ
      
    
    ⪡
    
      
        c
        R
      
      .
    
  




Note that in these embodiments, by making above modifications, the final solution will deviate from the optimal solution directly derived from Equation (10).

Additionally, to take a peer's evolving buffer level into consideration, peers also convert their buffer level into an equivalent received rate. Some embodiments of the P2P rate allocation system and method use a simple linear conversion where the aggregate received rate is incremented by,




  
    
      buffer
      ⁢
      
        
      
      ⁢
      level
    
    
      buffer
      ⁢
      
        
      
      ⁢
      length
    
  
  ·
  
    R
    .
  




For instance, suppose Peer j's buffer can hold 160 packets. Then if Peer j has 40 packets in the buffer, it will consider its aggregate received rate to be,





  
    
      1
      4
    
    ⁢
    R
  
  +
  
    
      ∑
      
        i
        ∈
        
          S
          j
        
      
    
    ⁢
    
      
        x
        
          i
          ,
          j
        
      
      .
    
  




III.B Minimum Cost Flow Formulation Module

As mentioned above, the minimum cost flow formulation module 150 uses a minimum cost flow formulation technique. This technique is used to determine the rate allocation by minimizing a cost flow formulation. An overview of this minimum cost flow formulation technique first will be discussed, followed by a detailed explanation of the technique.
FIG. 7 is a flow diagram illustrating the detailed operation of embodiments of the minimum cost flow formulation module 150 shown in FIG. 1. The operation of the module 150 begins by formulating a linear program formulation for determining rates (box 700). Next, the module 150 reformulates the linear program formulation as a minimum cost flow formulation (box 705).
In some embodiments of the minimum cost flow formulation module 150 a bipartite graph is used to reformulate the linear program formulation as a minimum cost flow formulation (box 710). As explained in detail below, the bipartite graph is denoted as G=(V, E), where V is a vertices set and E is an edges set (box 715). Next, the module 150 defines the vertices set V (box 720) and also defines the edges set E (box 725).
The module 150 then rewrites the minimum cost flow formulation in distributed form as a distributed minimum cost flow formulation (box 730). The distributed minimum cost flow formulation is updated using a rate allocation update equation (box 735). The distributed minimum cost flow formulation then is solved using Lagrangian multipliers (box 740). Solving involves updating the supply node Lagrangian multipliers using a supply node update equation (box 745) and updating the demand node Lagrangian multipliers using a demand node update equation (box 750). The module 150 also allocates surplus upload bandwidth using a convex program formulation (box 755). The module 150 then outputs an optimal rate allocation that minimizes ISP-unfriendly traffic (box 760). Each of these steps is described in detail as follows.
III.B.1 Notation
Peers are labeled in increasing order of arrival time as {1; 2; : : : ; n}. Server has label of 0. Let xi,j denote the rate from Peer i to Peer j (0≦i≦j<n), Let Ui be the uplink constraint of Peer i. The available server uplink capacity is given by U0. Let a be the associated cost of inter-ISP traffic and let β be the cost of serving traffic from the origin server. Let Ri denote the desired download rate of Peer i. This could vary with time as a function of the buffered amount of video.
IIII.B.2. Linear Programming Formulation
The linear programming formulation for determining the rates xi,j is given by:



  
    
      minimize
      ⁢
      
        
      
      ⁢
      
        
          ∑
          i
        
        ⁢
        
          β
          ⁢
          
            
          
          ⁢
          
            x
            
              0
              ⁢
              
                
              
              ⁢
              i
            
          
        
      
    
    +
    
      
        ∑
        
          i
          ,
          
            j
            ⁢
            
              
            
            ⁢
            cross
            ⁢
            
              
            
            ⁢
            ISP
          
        
      
      ⁢
      
        α
        ⁢
        
          
        
        ⁢
        
          x
          
            i
            ,
            j
          
        
      
    
  




  
    subject
    ⁢
    
      
    
    ⁢
    to
  




  
    
      
        ∑
        
          j
          >
          i
        
      
      ⁢
      
        x
        ij
      
    
    ≤
    
      
        U
        i
      
      ⁢
      
        
      
      ⁢
      
        ∀
        
          0
          ≤
          i
          ≤
          
            n
            .
          
        
      
    
  



III.B.3. Minimum Cost Flow Formulation

The above linear program can be formulated as a minimum cost flow problem as follows. First, construct a bipartite graph G=(V,E) where vertices,

V=V1∪V2 

are given by:

V1={v0, v1 . . . vn}

V2={v′1, v′2 . . . v′n, v′n+1}.

Each vi node has a supply of Ui that denotes the uplink capacity of Peer i (i>0) or server i (i=0). Each v′j node has a demand of Rj that denotes the desired download rate of Peer i. The super-sink node v′n+1 in V2 is added to model the unused uplink capacity of all peers-node v′n+1 has a demand of




  
    (
    
      
        
          ∑
          
            i
            =
            0
          
          n
        
        ⁢
        
          U
          i
        
      
      -
      
        
          ∑
          
            j
            =
            1
          
          n
        
        ⁢
        
          R
          j
        
      
    
    )
  
  .




Imagine that unused uplink capacity of all peers is used to send flow to the super-sink node v′n+1. With this, the total supply of nodes in V1 equals the total demand of nodes in V2.

The edge set E is given by the following:

    
    
        1. Edges (vi,v′j) for all 0≦i≦j≦n, denoting data transfer from Peer i (or, Server for i=0) to Peer j.
        2. Edges (vi,v′n+1) for all 0≦i≦n, denoting unused uplink capacity of Peer i (or, Server for i=0).
    
    


FIG. 8 is a block diagram illustrating an exemplary embodiment of a bipartite graph used by embodiments of the ISP-friendly rate allocation system and method shown in FIG. 1. As shown in FIG. 8, edges (vi,v′j) (0≦i<j≦n) have cost=−α if i, j are within same ISP. Otherwise, the cost=−1. Edges (v1, v′n+1) have cost=0. Note that any linear objective function can be modeled here. The term ci,j is used to denote the cost of the link (vi, v′j) and the corresponding rate (flow) by xi,j.
This is a variation of the transportation problem that is a special case of the minimum cost network flow problem. As set forth below, a distributed algorithm is developed to address this problem.
III.B.4. Distributed Algorithm
The linear program can be rewritten using the minimum cost flow formulation as:



  
    minimize
    ⁢
    
      
    
    ⁢
    
      
        ∑
        
          i
          <
          j
        
      
      ⁢
      
        
          c
          ij
        
        ⁢
        
          x
          ij
        
      
    
  




  
    
      subject
      ⁢
      
        
      
      ⁢
      to
    
    ,
    
      

    
    ⁢
    
      
        
          ∑
          
            j
            >
            i
          
        
        ⁢
        
          x
          ij
        
      
      =
      
        
          U
          i
        
        ⁢
        
          ∀
          
            0
            ≤
            i
            ≤
            n
          
        
      
    
  




  
    
      
        ∑
        
          i
          <
          j
        
      
      ⁢
      
        x
        ij
      
    
    =
    
      {
      
        
          
            
              R
              i
            
          
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                1
              
              ≤
              i
              ≤
              n
            
          
        
        
          
            
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  n
                
                ⁢
                
                  U
                  i
                
              
              -
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  n
                
                ⁢
                
                  R
                  i
                
              
            
          
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                i
              
              =
              
                n
                +
                1
              
            
          
        
      
    
  



Let m denote the number of variables xi,j, In other words, m denotes the number of edges in the bipartite graph G. The above constraints can be written in compact notation as Ax=b, where A is the (2n+2)×m matrix containing 0's and 1's corresponding to the LHS of the constraints, column vector x={xi,j}, and column vector b=(U0,U1 . . . Un,D1,D2, . . . , Dn+1) where




  
    D
    j
  
  =
  
    {
    
      
        
          
            R
            j
          
        
        
          
            
              if
              ⁢
              
                
              
              ⁢
              1
            
            ≤
            j
            ≤
            n
          
        
      
      
        
          
            
              
                ∑
                
                  i
                  =
                  0
                
                n
              
              ⁢
              
                U
                i
              
            
            -
            
              
                ∑
                
                  i
                  =
                  1
                
                n
              
              ⁢
              
                R
                i
              
            
          
        
        
          
            
              if
              ⁢
              
                
              
              ⁢
              j
            
            =
            
              n
              +
              1
            
          
        
      
    
  




The linear program now becomes

minimizeΣi<jci,jxi,j 

subject to

Ax=b,

III.B.4.i Augmented Lagrangian Technique

Embodiments of the minimum cost flow formulation module 150 use an augmented Lagrangian technique to solve the equivalent problem below for some parameter α>0. This problem can be represented as,

minimizeΣi<jci,jxi,j+ 1/2a∥Ax−b∥2 

subject to

Ax=b

Lagrangian multipliers are introduced and given by the column vector p for the constraints Ax=b. The augmented Lagrangian function is given by,




  
    
      L
      a
    
    ⁡
    
      (
      
        x
        ,
        p
      
      )
    
  
  =
  
    
      
        ∑
        
          i
          <
          j
        
      
      ⁢
      
        
          c
          ij
        
        ⁢
        
          x
          ij
        
      
    
    +
    
      
        p
        T
      
      ⁡
      
        (
        
          Ax
          -
          b
        
        )
      
    
    +
    
      
        1
        2
      
      ⁢
      a
      ⁢
      
        
          
          
            Ax
            -
            b
          
          
        
        2
      
    
  




The Lagrangian sub-problem is given by,

minimizexLa(x,p).

Embodiments of the minimum cost flow formulation module 150 next solve the Lagrangian sub-problem including the Gauss-Siedel method for quadratic programming. The solution will define the values of the xi,j variables in next time step and will be denoted by x(t+1). The values of the Lagrangian multipliers p then are given by,

p(t+1)=p(t)+a(t)[Ax(t+1)−b]

Lagrangian multipliers p next are divided into two groups: (i) si, 0≦i≦n corresponding to supply nodes, and (ii) dj, 1≦j≦n+1 corresponding to demand nodes in the bipartite graph. Then, the update step for the Lagrangian multipliers becomes,




  
    
      
        
          
            
              s
              i
            
            ⁡
            
              (
              
                t
                +
                1
              
              )
            
          
          =
          
            
              
                s
                i
              
              ⁡
              
                (
                t
                )
              
            
            -
            
              
                a
                ⁡
                
                  (
                  t
                  )
                
              
              [
              
                
                  U
                  i
                
                -
                
                  
                    ∑
                    
                      k
                      >
                      i
                    
                  
                  ⁢
                  
                    
                      x
                      ik
                    
                    ⁡
                    
                      (
                      
                        t
                        +
                        1
                      
                      )
                    
                  
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        24
        )
      
    
  




which is the supply node update equation, and





  
    
      
        
          
            
              d
              j
            
            ⁡
            
              (
              
                t
                +
                1
              
              )
            
          
          =
          
            
              
                d
                j
              
              ⁡
              
                (
                t
                )
              
            
            -
            
              
                a
                ⁡
                
                  (
                  t
                  )
                
              
              [
              
                
                  D
                  j
                
                -
                
                  
                    ∑
                    
                      k
                      <
                      j
                    
                  
                  ⁢
                  
                    
                      x
                      kj
                    
                    ⁡
                    
                      (
                      
                        t
                        +
                        1
                      
                      )
                    
                  
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        25
        )
      
    
  




which is the demand node update equation. Here, any non-decreasing sequence of positive numbers can be used for a(t).

III.B.4.ii. Solving the Lagrangian Sub-Problem

After algebraic simplification, the augmented Lagrangian function can be written as,




  
    
      L
      a
    
    ⁡
    
      (
      
        x
        ,
        p
      
      )
    
  
  =
  
    
      
        1
        2
      
      ⁢
      
        ax
        T
      
      ⁢
      
        A
        T
      
      ⁢
      Ax
    
    +
    
      
        (
        
          c
          +
          
            
              A
              T
            
            ⁢
            p
          
          -
          
            
              aA
              T
            
            ⁢
            b
          
        
        )
      
      ⁢
      
        x
        .
      
    
  




The standard form of the quadratic program for applying Gauss-Siedel method is,




  
    
      
        minimize
        n
      
      ⁢
      
        1
        2
      
      ⁢
      
        u
        T
      
      ⁢
      Pu
    
    +
    
      
        r
        T
      
      ⁢
      u
    
  




  where




  
    
      P
      =
      
        
          A
          ′
        
        ⁢
        
          Q
          
            -
            1
          
        
        ⁢
        
          A
          
            ′
            ⁢
            
              
            
            ⁢
            T
          
        
      
    
    ,
    
      r
      =
      
        
          c
          ′
        
        -
        
          
            A
            ′
          
          ⁢
          
            Q
            
              -
              1
            
          
          ⁢
          
            
              b
              ′
            
            .
          
        
      
    
  



Casting this to the form of the Lagrangian sub-problem that is to be solved yields,

A′=AT, Q=aT, r=c+AT(p−ab)

Let column vector a′j be the j-th column of A′T, and let column vector wj be the j-th row of A′Q−1. Then, the solution of the above (standard form) quadratic program is given by,





  
    
      u
      j
    
    =
    
      
        [
        
          
            u
            j
          
          -
          
            
              1
              
                
                  u
                  j
                  
                    ′
                    ⁢
                    
                      
                    
                    ⁢
                    T
                  
                
                ⁢
                
                  a
                  j
                  ′
                
              
            
            ⁢
            
              (
              
                
                  r
                  j
                
                -
                
                  
                    w
                    j
                    T
                  
                  ⁢
                  y
                
              
              )
            
          
        
        ]
      
      +
    
  
  ,




where y=-A-Tu.

Note that the matrix A contains exactly two 1's in each column. Using this structure of matrix A, the above equation can be simplified to obtain the update step for the original xi,j variables as follows,




  
    
      x
      ij
    
    ⁡
    
      (
      
        t
        +
        1
      
      )
    
  
  =
  
    
      
        x
        ij
      
      ⁡
      
        (
        t
        )
      
    
    -
    
      
        1
        
          2
          ⁢
          
            
          
          ⁢
          
            a
            ⁡
            
              (
              t
              )
            
          
        
      
      ⁡
      
        [
        
          
            c
            ij
          
          +
          
            
              s
              i
            
            ⁡
            
              (
              t
              )
            
          
          +
          
            
              d
              j
            
            ⁡
            
              (
              t
              )
            
          
        
        ]
      
    
    +
    
      
        1
        2
      
      [
      
        
          U
          i
        
        -
        
          
            ∑
            
              k
              >
              i
            
          
          ⁢
          
            
              x
              ik
            
            ⁡
            
              (
              t
              )
            
          
        
      
      ]
    
    +
    
      
        
          1
          2
        
        [
        
          
            D
            j
          
          -
          
            
              ∑
              
                k
                <
                j
              
            
            ⁢
            
              
                x
                kj
              
              ⁡
              
                (
                t
                )
              
            
          
        
        ]
      
      .
    
  




Since 0≦xi,j≦min(Ui, Dj), the xi,j value obtained at every step using the above equation needs to be projected to the interval [0; min(Ui, Dj)].

III.B.4.iii. Putting it All Together

Equation (26) can be used to update the rates xi,j(t) locally, since it requires communication between nodes i and j only. The term,

(Ui=Σk>ixik),

is the remaining unused uplink capacity of Peer i (or Server if i=0). The term,

(Dj−Σk<jxkj),

is the additional rate that Peer j desires.

Equations (24) and (25) can be used to update the Lagrangian multipliers. Each peer node owns two such variables and the update steps can be performed locally by using the value of the remaining unused uplink capacity of the peer or the value of the additional rate that the peer desires.
III.B.5. Allocating Surplus Upload Capacities
As mentioned earlier, the desired streaming rates Ri for each user could vary over time as a function of the buffered size of video. If excess uplink capacity is available in the system after allocating uplink bandwidth as per the earlier formulation, then embodiments of the minimum cost flow formulation module 150 can use this to provide additional download rates top users.
There could be many ways of allocating surplus bandwidth among competing users in the system. One possible objective might be to minimize the total (remaining) time of all peers to download the remaining (respective) portions of the video or to build some desired (large) buffered video size as an insurance against variability of network conditions in the future, Let zi be the additional download rate assigned to user i and let the size of the video remaining it needs to download be gi. Let wi be the surplus uplink capacity at node i (this could vary over time but is known at node i). Then, the instantaneous remaining download time of user i is given by,




  
    
      
        f
        i
      
      ⁡
      
        (
        
          z
          i
        
        )
      
    
    =
    
      
        
          
            g
            i
          
          
            
              z
              i
            
            +
            
              R
              i
            
          
        
        ⁢
        
          ∀
          i
        
      
      =
      1
    
  
  ,
  ...
  ⁢
  
    
  
  ,
  
    n
    .
  




Note that the functions fi(zi) are convex. In order to avoid increasing the objective function value defined in the earlier formulation, surplus upload capacity is allowed to be allocated only to peers within the same ISP (and thus use only intra-ISP links). This leads to the convex program formulation for allocating surplus uplink bandwidth:



  
    minimize
    ⁢
    
      
    
    ⁢
    
      
        ∑
        
          i
          =
          1
        
        n
      
      ⁢
      
        
          f
          i
        
        ⁡
        
          (
          
            z
            i
          
          )
        
      
    
  




  
    subject
    ⁢
    
      
    
    ⁢
    to
  




  
    
      
        ∑
        
          i
          <
          j
        
        
          
        
      
      ⁢
      
        y
        ij
      
    
    =
    
      
        z
        j
      
      ⁢
      
        ∀
        
          1
          ≤
          j
          ≤
          n
        
      
    
  




  
    
      
        ∑
        
          j
          >
          i
        
      
      ⁢
      
        y
        ij
      
    
    ≤
    
      
        w
        i
      
      ⁢
      
        ∀
        
          1
          ≤
          i
          ≤
          
            n
            .
          
        
      
    
  



This convex program formulation can be solved in a distributed manner using primal, primal-dual, or dual (subgradient) methods. Note that this and the previous optimization formulations will run simultaneously in the system. As a result of being allocated surplus bandwidth over some period of time, a user's buffered video size could increase over time, thus allowing it to decrease its desired streaming rate Ri. This, in turn may increase the surplus upload capacity in the system available for use by the convex program formulation. Thus, there is close interaction over time between execution of the two optimization formulations through the quantities Ri and wi.
III.C Packet-Level Scheduling Module
Thus far, it has been assumed that the rate allocation stays at an ideal network flow level. In other words, the upload bandwidth is treated as fluid and can be divided and utilized as finely as desired. The packet-level scheduling module 160 uses a practical packet-level algorithm that conforms nicely to the flow-level one.
FIG. 9 is a flow diagram illustrating the detailed operation of embodiments of the packet-level scheduling module 160 shown in FIG. 1. The operation of the module 160 begins by defining a rate allocation interval (box 900). Next, the optimal rate allocation is updated at each rate allocation interval using the update equation given above (box 910). A token allocation interval then is defined (box 920). At each token allocation interval, the peer is supplied with a number of tokens to give out (box 930). The number of tokens is proportional to the upload capacity of the peer. The module 160 then allocates at least some of the supplied tokens to a neighbor peer using the update equation (box 940).
III.C.1. Packet-Level Scheduling Technique
Embodiments of the packet-level scheduling module 160 use a sliding window. In particular, a sliding buffer is used for packet management. Peers maintain a buffer of interest and only download packets within the range. They slide the buffer forward every fixed interval. In some embodiments of the packet-level scheduling module 160, the fixed interval is 1 second. Peers cache all the content they have watched and make it available to their neighbor peers. In addition, peers advertise about their packet availability approximately once per second.
Embodiments of the packet-level scheduling module 160 also use packet request. More specifically, a peer ranks all the missing packets in its buffer by their importance. From each neighbor peer, the peer requests the most important packet possessed by the neighbor peer. This occurs as long as the number of pending requests at that neighbor peer is below a certain threshold. Embodiments of the packet-level scheduling module 160 determine the importance of a packet as follows.
The module 160 treats one second's worth of packets as one segment, and a segment is more important if it is closer to the playback deadline. Within a segment, the rarer a packet is in the peer's one-hop neighborhood the more important that packet. The module 160 also defines an urgency buffer, which is the first u second(s) of the buffer. The peers will simply request all the missing packets from the server. In some embodiments of the packet-level scheduling module 160, u=1. Unlike in the flow-level framework, peers are not prohibited from requesting packets from their downstream neighbors. Removing this restriction facilitates sharing among peers when their playback points are close.
Embodiments of the packet-level scheduling module 160 also use rate allocation and request response. In addition, the packet-level scheduling technique addresses how serving peers determine which neighbor peer's request to satisfy first. The module 160 uses a token system to implement the rate allocation algorithm at a packet level. A rate-allocation interval TRA and a token-allocation interval TTA are defined. Moreover, every TRA second, Peer i carries out rate update as in using the update equation given by Equation (23) above. The module 160 then computes the appropriate rate to allocate to each of its neighbor peers. This computation need not be synchronized among peers.
Every TTA second, Peer i gets a certain number of tokens to give out. It should be noted that peers need not synchronize their intervals. The number of tokens the peer obtains is proportional to its upload bandwidth. For instance, if Peer i has Ui kbps upload bandwidth, it gets Ui tokens each round. The peer then allocates xi,j tokens to neighbor Peer j as computed in the update equation given by Equation (23) above.
Note that Ui and xi,j need not be integers. Any partial tokens will simply be left with peers for use in future rounds. In deciding which neighbor peer's request to satisfy, Peer i chooses the neighbor peer with the highest level of tokens and deducts




  
    pocket
    ⁢
    
      
    
    ⁢
    size
  
  
    T
    TA
  




tokens from the neighbor peer.

Embodiments of the packet-level scheduling module 160 also perform request rejection. Specifically, Peer i should turn down requests that it is unlikely to fulfill, so that the requester does not wait unnecessarily. For this purpose, peers include time to playback deadline in the request of each packet (Dm). For example, suppose a neighbor Peer j has accumulated token level Lj (Peer i has this information). Peer i will turn down the request for Piece m from Peer j if by the time the packet expires, Peer j still would not have accumulated enough token. This can be represented by the equation,




  
    
      L
      j
    
    +
    
      
        x
        
          i
          ,
          j
        
      
      ·
      
        ⌊
        
          
            D
            m
          
          
            T
            TA
          
        
        ⌋
      
    
  
  <
  
    
      
        packet
        ⁢
        
          
        
        ⁢
        size
      
      
        T
        TA
      
    
    .
  




III.C.2. Computation Complexity and Overhead

Each iteration of the rate allocation at Peer i involves computing the update step for each of its neighbor peer j. Since the derivatives of the utility and cost functions are closed form expressions and pre-stored, this is a fairly straight forward computation with a few multiplications and additions. Each peer typically has no more than 20 neighbor peers including both upstream and downstream neighbor peers. Typically, the rate allocation is carried out once per second (such that TRA=1 second). For personal computers having GHz processors the computation overhead is marginal.
There is almost no communication overhead involved either. This is because the only additional information that peers have to transmit to each other is: (1) their aggregate received bandwidth; and, (2) the playback deadline for each packet. Peers can piggyback this information in their packet requests with almost no overhead.
IV. Exemplary Operating Environment
Embodiments of the ISP-friendly rate allocation system 100 and method are designed to operate in a computing environment. The following discussion is intended to provide a brief, general description of a suitable computing environment in which embodiments of the ISP-friendly rate allocation system 100 and method may be implemented.
FIG. 10 illustrates an example of a suitable computing system environment in which embodiments of the ISP-friendly rate allocation system and method shown in FIGS. 1-9 may be implemented. The computing system environment 1000 is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality of the invention. Neither should the computing environment 1000 be interpreted as having any dependency or requirement relating to any one or combination of components illustrated in the exemplary operating environment.
Embodiments of the ISP-friendly rate allocation system 100 and method are operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well known computing systems, environments, and/or configurations that may be suitable for use with embodiments of the ISP-friendly rate allocation system 100 and method include, but are not limited to, personal computers, server computers, hand-held (including smartphones), laptop or mobile computer or communications devices such as cell phones and PDA's, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
Embodiments of the ISP-friendly rate allocation system 100 and method may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform particular tasks or implement particular abstract data types. Embodiments of the ISP-friendly rate allocation system 100 and method may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices. With reference to FIG. 10, an exemplary system for embodiments of the ISP-friendly rate allocation system 100 and method includes a general-purpose computing device in the form of a computer 1010. It should be noted that the computing device 110 shown in FIG. 1 is one implementation of the computer 1010.
Components of the computer 1010 may include, but are not limited to, a processing unit 1020 (such as a central processing unit, CPU), a system memory 1030, and a system bus 1021 that couples various system components including the system memory to the processing unit 1020. The system bus 1021 may be any of several types of bus structures including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus also known as Mezzanine bus.
The computer 1010 typically includes a variety of computer readable media. Computer readable media can be any available media that can be accessed by the computer 1010 and includes both volatile and nonvolatile media, removable and non-removable media. By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes volatile and nonvolatile removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data.
Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the computer 1010. By way of example, and not limitation, communication media includes wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above should also be included within the scope of computer readable media.
The system memory 1030 includes computer storage media in the form of volatile and/or nonvolatile memory such as read only memory (ROM) 1031 and random access memory (RAM) 1032. A basic input/output system 1033 (BIOS), containing the basic routines that help to transfer information between elements within the computer 1010, such as during start-up, is typically stored in ROM 1031. RAM 1032 typically contains data and/or program modules that are immediately accessible to and/or presently being operated on by processing unit 1020. By way of example, and not limitation, FIG. 10 illustrates operating system 1034, application programs 1035, other program modules 1036, and program data 1037.
The computer 1010 may also include other removable/non-removable, volatile/nonvolatile computer storage media. By way of example only, FIG. 10 illustrates a hard disk drive 1041 that reads from or writes to non-removable, nonvolatile magnetic media, a magnetic disk drive 1051 that reads from or writes to a removable, nonvolatile magnetic disk 1052, and an optical disk drive 1055 that reads from or writes to a removable, nonvolatile optical disk 1056 such as a CD ROM or other optical media.
Other removable/non-removable, volatile/nonvolatile computer storage media that can be used in the exemplary operating environment include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like. The hard disk drive 1041 is typically connected to the system bus 1021 through a non-removable memory interface such as interface 1040, and magnetic disk drive 1051 and optical disk drive 1055 are typically connected to the system bus 1021 by a removable memory interface, such as interface 1050.
The drives and their associated computer storage media discussed above and illustrated in FIG. 10, provide storage of computer readable instructions, data structures, program modules and other data for the computer 1010. In FIG. 10, for example, hard disk drive 1041 is illustrated as storing operating system 1044, application programs 1045, other program modules 1046, and program data 1047. Note that these components can either be the same as or different from operating system 1034, application programs 1035, other program modules 1036, and program data 1037. Operating system 1044, application programs 1045, other program modules 1046, and program data 1047 are given different numbers here to illustrate that, at a minimum, they are different copies. A user may enter commands and information (or data) into the computer 1010 through input devices such as a keyboard 1062, pointing device 1061, commonly referred to as a mouse, trackball or touch pad, and a touch panel or touch screen (not shown).
Other input devices (not shown) may include a microphone, joystick, game pad, satellite dish, scanner, radio receiver, or a television or broadcast video receiver, or the like. These and other input devices are often connected to the processing unit 1020 through a user input interface 1060 that is coupled to the system bus 1021, but may be connected by other interface and bus structures, such as, for example, a parallel port, game port or a universal serial bus (USB). A monitor 1091 or other type of display device is also connected to the system bus 1021 via an interface, such as a video interface 1090. In addition to the monitor, computers may also include other peripheral output devices such as speakers 1097 and printer 1096, which may be connected through an output peripheral interface 1095.
The computer 1010 may operate in a networked environment using logical connections to one or more remote computers, such as a remote computer 1080. The remote computer 1080 may be a personal computer, a server, a router, a network PC, a peer device or other common network node, and typically includes many or all of the elements described above relative to the computer 1010, although only a memory storage device 1081 has been illustrated in FIG. 8. The logical connections depicted in FIG. 10 include a local area network (LAN) 1071 and a wide area network (WAN) 1073, but may also include other networks. Such networking environments are commonplace in offices, enterprise-wide computer networks, intranets and the Internet.
When used in a LAN networking environment, the computer 1010 is connected to the LAN 1071 through a network interface or adapter 1070. When used in a WAN networking environment, the computer 1010 typically includes a modem 1072 or other means for establishing communications over the WAN 1073, such as the Internet. The modem 1072, which may be internal or external, may be connected to the system bus 1021 via the user input interface 1060, or other appropriate mechanism. In a networked environment, program modules depicted relative to the computer 1010, or portions thereof, may be stored in the remote memory storage device. By way of example, and not limitation, FIG. 10 illustrates remote application programs 1085 as residing on memory device 1081. It will be appreciated that the network connections shown are exemplary and other means of establishing a communications link between the computers may be used.
The foregoing Detailed Description has been presented for the purposes of illustration and description. Many modifications and variations are possible in light of the above teaching. It is not intended to be exhaustive or to limit the subject matter described herein to the precise form disclosed. Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims appended hereto.

## Claims

1. A method for rate allocation in a peer-to-peer (P2P) network, comprising:
determining a first rate allocation that minimizes a load on a server in communication with the P2P network to generate a minimum server load;
determining a second rate allocation that minimizes network traffic across Internet service provider (ISP) boundaries;
constructing a utility function using the first rate allocation and the second rate allocation; and
optimizing the utility function to obtain an optimal rate allocation that minimizes network traffic across ISP boundaries while keeping the minimum server load unchanged.

2. The method of claim 1, further comprising:
determining that spare upload bandwidth is available for uploading; and
maximizing peer prefetching using the spare upload bandwidth while still minimizing network traffic across ISP boundaries and keeping the minimum server load unchanged.

3. The method of claim 1, further comprising representing the utility function to be optimized as:




  
    
      
        max
        
          x
          
            i
            ,
            j
          
        
      
      ⁢
      
        
          ∑
          
            1
            ≤
            j
            ≤
            n
          
        
        ⁢
        
          
            (
            
              
                f
                ⁡
                
                  (
                  
                    x
                    j
                  
                  )
                
              
              -
              
                
                  ∑
                  
                    
                      i
                      ∈
                      
                        S
                        j
                      
                    
                    ⁢
                    
                      
                    
                  
                
                ⁢
                
                  
                    g
                    
                      i
                      ,
                      j
                    
                  
                  ⁡
                  
                    (
                    
                      
                        x
                        i
                      
                      ,
                      j
                    
                    )
                  
                
              
            
            )
          
          ⁢
          
            
          
          ⁢
          
            s
            .
            t
            .
            
              
                ∑
                
                  j
                  ∈
                  
                    D
                    i
                  
                
              
              ⁢
              
                x
                
                  i
                  ,
                  j
                
              
            
          
        
      
    
    ≤
    
      U
      i
    
  
  ,
  
    
      x
      
        i
        ,
        j
      
    
    ≥
    
      0
      ⁢
      
        
      
      ⁢
      
        ∀
        i
      
    
  
  ,
  
    j
    .
  
  ,




where gi,j(xi,j) is a cost for Peer i to upload to Peer j at a rate xi,j, Ui is an upload capacity for Peer i, and Di denotes a set of peers downloading from Peer i.

4. The method of claim 1, further comprising computing the first rate allocation using the equation:


U0min=Σjx*0,j,


where,


{x*0,j}jn=1


is an optimal solution for the first rate allocation that yields the minimum server load.

5. The method of claim 1, further comprising computing the second rate allocation using the equation:


Gmin=Σi,jgi,j(x*i,j),


where gi,j(xi,j) is a cost for Peer i to upload to Peer j at a rate xi,j.

6. The method of claim 3, further comprising:
optimizing the utility function by connecting two linear components of the utility function with a concave smooth curve;
determining the optimal rate allocation by solving the optimized utility function.

7. The method of claim 2, further comprising determining the rate allocation using a minimum cost flow formulation.

8. The method of claim 7, further comprising:
rewriting the minimum cost flow formulation as a distributed minimum cost flow formulation that is a distributed solution executed by each peer in the peer-to-peer network independently, while collectively achieving the optimal rate allocation; and
using a bipartite graph in the distributed minimum cost flow formulation.

9. A method for reducing network traffic across Internet service provider (ISP) boundaries in a peer-to-peer (P2P) network, comprising:
generating a utility function that represents a minimum load on a server in communication with the P2P network, a minimum network traffic across ISP boundaries while keeping the minimum load unchanged, and a maximum peer prefetching using spare upload bandwidth while still minimizing network traffic across ISP boundaries and keeping the minimum server load unchanged, the utility function represented as:





  
    
      
        max
        
          x
          
            i
            ,
            j
          
        
      
      ⁢
      
        
          ∑
          
            1
            ≤
            j
            ≤
            n
          
        
        ⁢
        
          
            (
            
              
                f
                ⁡
                
                  (
                  
                    x
                    j
                  
                  )
                
              
              -
              
                
                  ∑
                  
                    
                      i
                      ∈
                      
                        S
                        j
                      
                    
                    ⁢
                    
                      
                    
                  
                
                ⁢
                
                  
                    g
                    
                      i
                      ,
                      j
                    
                  
                  ⁡
                  
                    (
                    
                      
                        x
                        i
                      
                      ,
                      j
                    
                    )
                  
                
              
            
            )
          
          ⁢
          
            
          
          ⁢
          
            s
            .
            t
            .
            
              
                ∑
                
                  j
                  ∈
                  
                    D
                    i
                  
                
              
              ⁢
              
                x
                
                  i
                  ,
                  j
                
              
            
          
        
      
    
    ≤
    
      U
      i
    
  
  ,
  
    
      x
      
        i
        ,
        j
      
    
    ≥
    
      0
      ⁢
      
        
      
      ⁢
      
        ∀
        i
      
    
  
  ,
  
    j
    .
  
  ,




where gi,j(xi,j) is a cost for Peer i to upload to Peer j at a rate xi,j, Ui is an upload capacity for Peer i, and Di denotes a set of peers downloading from Peer i; and
optimizing the utility function to obtain an optimal rate allocation for the P2P network that minimizes the network traffic across ISP boundaries.

10. The method of claim 9, further comprising:
applying a steepest descent technique to the utility function;
initialize xi,j to zero; and
update at each step using an update equation given by,





  
    
      x
      .
    
    
      i
      ,
      j
    
  
  =
  
    
      Δ
      ·
      
        
          (
          
            
              
                
                  ∂
                  
                    
                  
                
                
                  ∂
                  
                    x
                    
                      i
                      ,
                      j
                    
                  
                
              
              ⁢
              
                f
                ⁡
                
                  (
                  
                    x
                    j
                  
                  )
                
              
            
            -
            
              
                
                  ∂
                  
                    
                  
                
                
                  ∂
                  
                    x
                    
                      i
                      ,
                      j
                    
                  
                
              
              ⁢
              
                
                  g
                  
                    i
                    ,
                    j
                  
                
                ⁡
                
                  (
                  
                    x
                    
                      i
                      ,
                      j
                    
                  
                  )
                
              
            
          
          )
        
        .
        
          

        
        ⁢
        
          x
          
            i
            ,
            j
          
        
      
    
    =
    
      
        [
        
          
            x
            
              i
              ,
              j
            
          
          +
          
            
              x
              .
            
            
              i
              ,
              j
            
          
        
        ]
      
      +
    
  




where [xi,j+{dot over (x)}i,j]+ denotes an I2 projection on to a feasible set in order to guarantee convergence.

11. The method of claim 10, further comprising modifying f(xj) by connecting two linear components of f(xj) with a concave smooth curve of width ΔR, where R is a desirable streaming rate of a peer in the P2P network.

12. The method of claim 10, further comprising:
defining a rate-allocation interval as TRA; and
using Peer i to update the optimal rate allocation at every TRA second using the update equation to determine which neighboring peer's request to satisfy first.

13. The method of claim 12, further comprising:
defining a token-allocation interval as TTA;
every TTA second, supplying Peer i with a number of tokens to give out such that the number of tokens is proportional to Peer i's upload bandwidth; and
allocating xi,j tokens to neighbor Peer j using the update equation.

14. A method for reducing Internet service provider (ISP)-unfriendly traffic in a peer-to-peer (P2P) network, comprising:
formulating a linear program formulation for determining rates denoted by xi,j, where xi,j are rates to upload from Peer i to Peer j;
reformulating the linear program formulation as a minimum cost flow formulation;
rewriting the minimum cost flow formulation as a distributed minimum cost flow formulation that is a distributed solution executed by each peer in the peer-to-peer network independently, while collectively achieving an optimal rate allocation;
updating the distributed minimum cost flow formulation using a rate allocation update equation having supply node Lagrangian multipliers and demand node Lagrangian multipliers; and
solving the distributed minimum cost flow formulation using the supply node Lagrangian multipliers, the demand node Lagrangian multipliers, and a computing device having a processor to obtain the optimal rate allocation that minimizes the ISP-unfriendly traffic across ISP boundaries while keeping a minimum server load that minimizes a load on a server in communication with the P2P network unchanged to obtain the optimal rate allocation.

15. The method of claim 14, further comprising using a bipartite graph to reformulate the linear program formulation as the distributed minimum cost flow formulation.

16. The method of claim 15, further comprising:
denoting the bipartite graph as G=(V,E), where V is a set of vertices and E is a set of edges for the bipartite graph;
defining the vertices set V as V=V1∪V2 given by:


V1={v0, v1, . . . , vn} and V2={v′1,v′2, . . . , v′n,v′n+1}; and


defining the edges set E as:
(a) edges (vi, v′j) for all 0≦i≦j≦n, which denotes a data transfer from Peer I to Peer j;
(b) edges (vi, v′n+1) for all 0≦i≦n, which denotes an unused upload capacity of Peer i.

17. The method of claim 16, further comprising updating the distributed minimum cost flow formulation using the rate allocation update equation as given by,




  
    
      
        x
        ij
      
      ⁡
      
        (
        
          t
          +
          1
        
        )
      
    
    =
    
      
        
          x
          ij
        
        ⁡
        
          (
          t
          )
        
      
      -
      
        
          1
          
            2
            ⁢
            
              
            
            ⁢
            
              a
              ⁡
              
                (
                t
                )
              
            
          
        
        ⁡
        
          [
          
            
              c
              ij
            
            +
            
              
                s
                i
              
              ⁡
              
                (
                t
                )
              
            
            +
            
              
                d
                j
              
              ⁡
              
                (
                t
                )
              
            
          
          ]
        
      
      +
      
        
          1
          2
        
        [
        
          
            U
            
              i
              ⁢
              
                
              
            
          
          -
          
            
              ∑
              
                k
                >
                i
              
            
            ⁢
            
              
                x
                ik
              
              ⁡
              
                (
                t
                )
              
            
          
        
        ]
      
      +
      
        
          1
          2
        
        [
        
          
            D
            j
          
          -
          
            
              ∑
              
                k
                <
                j
              
            
            ⁢
            
              
                x
                kj
              
              ⁡
              
                (
                t
                )
              
            
          
        
        ]
      
    
  
  ,




where cij is a cost of a link (vi,v′j), si(t) are supply node Lagrangian multipliers, and dj(t) are demand node Lagrangian multipliers.

18. The method of claim 17, further comprising:
updating the supply node Lagrangian multipliers using the update equation:





  
    
      
        s
        i
      
      ⁡
      
        (
        
          t
          +
          1
        
        )
      
    
    =
    
      
        
          s
          i
        
        ⁡
        
          (
          t
          )
        
      
      -
      
        
          a
          ⁡
          
            (
            t
            )
          
        
        [
        
          
            U
            i
          
          -
          
            
              ∑
              
                k
                >
                i
              
            
            ⁢
            
              
                x
                ik
              
              ⁡
              
                (
                
                  t
                  +
                  1
                
                )
              
            
          
        
        ]
      
    
  
  ;
  and



updating the demand node Lagrangian multipliers using the update equation:





  
    
      
        d
        j
      
      ⁡
      
        (
        
          t
          +
          1
        
        )
      
    
    =
    
      
        
          d
          j
        
        ⁡
        
          (
          t
          )
        
      
      -
      
        
          a
          ⁡
          
            (
            t
            )
          
        
        [
        
          
            D
            j
          
          -
          
            
              ∑
              
                k
                <
                j
              
              
                
              
            
            ⁢
            
              
                x
                kj
              
              ⁡
              
                (
                
                  t
                  +
                  1
                
                )
              
            
          
        
        ]
      
    
  
  ,




where a(t) is any non-decreasing sequence of positive numbers.

19. The method of claim 17, further comprising allocating surplus upload bandwidth using a convex program formulation given by,




  
    minimize
    ⁢
    
      
    
    ⁢
    
      
        ∑
        
          i
          =
          1
        
        n
      
      ⁢
      
        
          f
          i
        
        ⁡
        
          (
          
            z
            i
          
          )
        
      
    
  
  ,




subject to,



  
    
      
        ∑
        
          i
          <
          j
        
      
      ⁢
      
        y
        ij
      
    
    =
    
      
        z
        j
      
      ⁢
      
        
      
      ⁢
      
        ∀
        
          1
          ≤
          j
          ≤
          n
        
      
    
  




  
    
      
        
          ∑
          
            j
            >
            i
          
        
        ⁢
        
          y
          ij
        
      
      ≤
      
        
          w
          i
        
        ⁢
        
          
        
        ⁢
        
          ∀
          
            1
            ≤
            i
            ≤
            n
          
        
      
    
    ,
  



where zi is an additional download rate assigned to Peer i, and wi is a surplus upload capacity at Peer i.